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ABSTRACT 


For  steady  motion  of  a  propeller  operating  in  an  inviscid  fluid  having  an  unbounded 
irrotational  flow  Field,  an  expression  for  the  velocity  potential  (in  excess  of  the  body 
motion)  is  derived  in  terms  of  the  boundary  values.  From  this  expression  perturbation  so¬ 
lutions  are  determined— one  for  small  thickness— or  camber-to-chord  ratio  and  one  for  small 
chord-to-diameter  ratio. 

The  first  problem  (lifting-surface  theory)  is  a  regular-perturbation  problem,  and  the 
second  (lifting-line  theory)  is  a  singular-perturbation  problem  which  requires  construction  of 
matched  asymptotic  expansions.  Two  terms  of  each  series  are  found.  Numerical  techniques 
are  not  discussed.  The  outer  solution  for  the  lifting  line  is  the  same  as  that  published  in 
the  literature.  The  formal  lifting-surface  analysis  differs  from  other  developments  in  several 
ways.  The  most  important  of  these  is  that  for  propellers  with  variable  pitch,  warp,  and 
rake,  the  normal  to  the  blade  has  a  radial  component  which  requires  consideration  of  the 
radial  velocity  in  determining  the  blade  shape.  For  the  case  considered  the  sign  of  the  con¬ 
tributions  in  the  inner  radii  differ  from  the  outer  radii  values;  thus,  this  additional  term 
might  cause  little  effect  on  thrust  but  could  be  important  for  cavitation  performance.  A 
design  procedure  is  discussed  which  involves  only  quantities  appropriate  for  the  lifting- 
surface  analysis. 
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INTRODUCTION 


A-OBJECTIVE 

Calculations  pertaining  to  lifting  bodies  are  based  on  linearized  versions  of  the  mathematical  model 
used  to  describe  the  flow  field  about  such  objects.  The  traditional  names  for  the  two  approximations 
used  are  lifting-surface  and  lifting-line.  Each  of  these  is  associated  with  a  different  small  parameter. 

These  parameters  are  discussed  in  Section  B. 

Until  recently,  developments  in  the  linearized  theory  have  usually  been  presented  on  an 
intuitive  basis,  rather  than  on  a  formal  mathematical  basis.  In  wing  theory,  the  formal  linearization  pro¬ 
cedures  have  not  altered  the  lower-order  theory  significantly;  however,  they  have  permitted  higher-order 
effects  to  be  calculated.  Perhaps  the  best  example  of  the  formalization  of  the  previous  intuitive  procedure 
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is  the  lifting-line  theory  presented  by  Friedrichs  and  Van  Dyke.  In  particular,  Van  Dyke  showed  that 
the  integral  equation  could  be  simplified  to  a  quadrature  and  that  higher-order  terms  could  be  found.  The 
contribution  of  Friedrichs  and  Van  Dyke  will  be  discussed  further  in  later  paragraphs  and  are  mentioned 
here  as  an  example  of  the  formalization. 

To  date,  a  similar  formal  development  appropriate  for  propeller  geometry  has  not  been  presented. 
Accordingly,  the  objective  of  the  present  work  is  to  develop  the  formal  linearized  solutions  for  propeller 
lifting-surface  theory  and  lifting-line  theory. 

B— REGULAR  AND  SINGULAR-PERTURBATION  PROBLEMS 

A  hubless  propeller  has  three  characteristic  lengths:  the  diameter,  a  representative  chordlength, 
and  a  representative  thickness  or  camber.  For  lifting-surface  theory,  the  ratio  of  thickness  and/or  camber 
to  either  chordlength  or  diameter  is  the  small  parameter,  and  the  exact  formulation  is  implicity  expanded 
in  a  perturbation  series.  Such  a  series  is  herein  called  a  regular-perturbation  series,  since  useful  information 
can  be  obtained  without  further  analysis.  However,  the  resulting  formulation  is  not  uniformly  valid  since 
it  fails  at  the  leading  edge.  The  lifting-surface  equations  are  derived  in  Chapter  3. 

Instead  of  the  thickness  and  camber  ratio  being  considered  small,  the  chord-to-diameter  ratio  can 
be  taken  as  the  perturbation  parameter.  This  leads  to  lifting-line  theory  which  is  the  subject  of  Chapter  4. 

Unfortunately  a  vanishing  chordlength  leads  to  a  singular-perturbation  problem  as  opposed  to  the 
regular-perturbation  problem  in  lifting-surface  theory.  It  is  often  sufficient  to  distinguish  between  them  by 
defining  the  regular  problem  as  that  for  which  the  order  of  the  governing  equations  and  the  number  of 
boundary  conditions  remain  fixed  as  the  parameter  goes  to  zero.  This  definition  is  adopted  here,  although 
an  additional  requirement  is  violated.  The  additional  requirement  is  that  the  resulting  formulation  must  be 
uniformly  valid  throughout  the  flow  field.  As  already  mentioned,  the  point  of  view  adopted  here  is  that 
useful  information  can  be  obtained  from  the  regular-perturbation  series  without  consideration  of  the 
singular  region. 


^Friedrichs,  K.O.,  “Special  Topics  in  Fluid  Dynamics,”  New  York  University  (1953);  also  published  by 
Gordon  and  Breach,  New  York  (1966).  (A  complete  listing  of  references  is  given  on  page(s)  115). 

^Van  Dyke,  M.,  “Lifting  Line  Theory  as  a  Singular  Perturbation  Problem,”  Archiwum  Mechaniki 
Stosowanej,  Vol.  3,  No.  16,  pp.  601-614  (1964). 

$ 

A  singular  region  can  be  anticipated  since  the  formulation  requires  a  derivative  of  the  shape  function. 
Although  the  shape  relative  to  the  reference  surface  is  uniformly  small,  the  derivative  of  the  thickness 
function  is  infinite  at  the  leading  edge  and  small  away  from  it. 
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For  a  singular-perturbation  problem,  the  order  of  the  governing  equation  is  decreased,  and/or  one 
or  more  of  the  boundary  conditions  has  to  be  discarded.  In  problems  involving  a  lifting  body  of  small 
chordlength,  the  limiting  process  as  the  chordlength  goes  to  zero  replaces  the  body  by  a  lifting  line,  and, 
hence,  the  body  boundary  conditions  cannot  be  applied.  (The  order  of  the  equations  remains  unchanged.) 
To  regain  the  details  of  the  flow  near  the  body,  this  region  must  be  magnified  by  stretching  the 
coordinates.  The  flow  field  near  the  body  is  calculated  in  these  stretched  variables  and  must  be  matched 
to  the  flow  determined  for  the  lifting  line.  The  mathematical  techniques  involved  in  the  analysis  are 
explained  by  Van  Dyke/*  Kaplun,^  and  Cole.^ 

The  first  application  of  singular-perturbation  concepts  to  wings  of  finite  span  was  made  by 
Prandtl.^  Later,  Friedrichs^  formalized  the  concept  of  inner  and  outer  regions  and  derived  Prandtl’s  inte- 
gral  equation  for  circulation.  Van  Dyke  obtained  higher-order  terms  in  the  perturbation  series  and  showed 
that  the  circulation  was  obtained  as  a  quadrature  rather  than  from  an  integral  equation.  A  singular- 

n 

perturbation  solution  has  also  been  obtained  for  a  swept  wing  by  Thurber  and  for  a  nonplanar  wing  by 

o  9 

Rotta.  The  formal  perturbation  series  assume  zero  thickness.  Van  Dyke  compared  his  calculations 
with  results  from  exact  theory  for  the  total  lift  of  elliptic  wings  and  found  good  agreement  for  aspect 
ratios  as  low  as  about  2.7.  He  was  able  to  identify  terms  in  the  solution  which  gave  streamline  curvature 
and  which  slightly  modified  the  angle  of  attack. 

For  marine  propellers  the  chord-to-diameter  ratio  is  about  unity,  and  the  thickness-to-chord  ratio 
is  about  one-tenth.  One  intuitively  expects  then  that  the  lifting-surface  formulation  will  be  more  accurate 
than  the  lifting-line  theory  for  the  same  number  of  terms. 

C— DESIGN  AND  PERFORMANCE  PROBLEMS 

In  the  present  work  a  distinction  is  made  between  the  design  problem  and  the  performance- 
evaluation  problem.  This  distinction  does  not  apply  to  the  mathematical  formulation  but  to  the  application 
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Van  Dyke,  M.,  “Perturbation  Methods  in  Fluid  Mechanics,”  Academic  Press,  Inc.,  New  York  (1964). 

^Kaplun,  S.,  “Fluid  Mechanics  and  Singular  Perturbations,”  Edited  by  P.A.  Lagerstrom,  et  al.,  Academic 
Press,  Inc.,  New  York  (1967). 

^Cole,  J.D.,  “Perturbation  Methods  in  Applied  Mathematics,”  Blaisdell  (1968). 

^Prandtl,  L.j  “Application  of  Modern  Hydrodynamics  to  Aeronautics,”  National  Advisory  Committee 
for  Aeronautics  Report  116  (1921). 

n 

Thurber,  J.K.,  “An  Asymptotic  Method  for  Determining  the  Lift  Distribution  of  a  Swept-Back  Wing  of 
Finite  Span,”  Communications  on  Pure  and  Applied  Mathematics,  Vol.  18,  pp.  733-756  (1965). 

'"’Rotta,  N.R.,  “The  Non-Planar,  Moderate  Aspect  Ratio,  Subsonic  Winds,”  Ph.D.  Thesis,  New  York 
University,  University  Microfilms  69-4576  (1968). 
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of  theoretical  results  since  certain  quantities  are  assumed  known  in  one  problem  while  they  would  be 
unknown  in  the  other. 

In  the  typical  design  problem,  one  specifies  the  blade  loading,  the  diameter,  number  of  blades, 
blade  outline  and  hub  details.  One  seeks  the  blade-section  thickness,  camber  and  pitch  angle  for  operation 
in  a  given  flow  field.  Variations  often  occur  in  the  specifications;  in  early  design  procedures,  the  camber¬ 
line  shape  was  specified,  and  the  radial  magnitude  needed  to  obtain  a  specified  thrust  was  sought.  In  some 
modern  design  techniques,  the  spanwise  and  chordwise  shape  of  the  loading  function  as  well  as  the  thrust 
is  specified,  and  the  camberline  shape  is  to  be  calculated.  In  both  of  these  design  problems,  the  thickness 
distribution  is  specified  from  strength  and  cavitation  considerations.  When  the  thickness  distribution  and 
loading  function  are  specified,  the  design  calculations  become  quadratures  in  the  lifting-surface  problem. 
This  point  is  elaborated  in  Chapter  3.  Recently  techniques  for  determining  the  thickness  distribution  from 

Q 

a  specified  pressure  distribution  have  been  developed  by  Hille/  This  would  allow  a  design  problem  for 
which  the  total  pressure  distribution  was  specified  and  for  which  the  thrust  and  geometry  would  be  calcu¬ 
lated.  However,  the  combination  of  strength  considerations  and  the  simplicity  of  the  design  based  on 
quadratures  will  probably  insure  continued  usage  of  thickness  specification. 

In  contrast  to  the  possible  variations  for  the  design  problem,  the  performance-evaluation  problem 
is  straightforward  to  describe.  The  geometry  is  specified  and  from  it  one  seeks  to  calculate  some  desired 
quantities  when  the  propeller  operates  in  a  given  upstream  flow  field.  A  formal  solution  of  this  problem 
is  obtained  as  a  quadrature  in  the  lifting-line  theory  of  Chapter  4.  Unfortunately,  an  insufficient  number 
of  terms  in  the  series  is  determined  to  permit  accurate  calculations,  so  investigators  in  this  area  currently 
attempt  a  solution  of  the  integral  equation  from  the  lifting-surface  analysis. 

D— HISTORICAL  DEVELOPMENT  OF  PROPELLER 
THEORY  AND  STATE  OF  THE  ART 

Propeller  theory  is  based  on  the  principles  applied  to  conventional  wings,  except  that  the  pro¬ 
peller  geometry  adds  a  significant  complication.  As  discussed  in  Section  B,  Prandtl^  laid  the  foundation 
for  lifting-line  analysis  of  wings.  In  the  same  paper  he  included  an  approximate  solution  for  propellers 
with  small  chord-to-diameter  ratios.  Several  years  later,  Goldstein^  obtained  a  more  accurate  solution  for 
a  propeller  lifting-line  with  a  specific  circulation  distribution;  however,  it  was  not  until  more  than  two 
decades  later  that  Moriya^  presented  a  derivation  of  the  induced  velocities  at  a  lifting-line  of  an  arbitrary 

^Hille,  R.,  “Bestimmung  der  Dickenlinie  von  Propellerfliigelprofilen  bel  vorgegebener  Druckverteilung,” 
Report  262  of  the  Institut  fur  Schiffbau  der  Universitat  Hamburg,  West  Germany  (1970). 

10Goldstein,  S.,  “On  the  Vortex  Theory  of  Screw  Propellers,”  Procedure  Royal  Society  of  London, 

Series  A,  Vol.  123,  pp.  440-465  (1929). 

^Moriya,  T.,  “On  the  Integration  of  Biot-Savart’s  Law  in  Propeller  Theory,”  (in  Japanese),  Journal  of 
the  Society  for  Aeronautical  Science,  Japan,  Vol.  9,  No.  89,  of  1015-1020  (1942);  English  translation  in 
Selected  Scientific  and  Technical  Papers  by  Tomijiro  Moriya,  Moriya  Memorial  Committee,  University  of 
Tokyo,  Japan,  pp.  74-80  (1959). 
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propeller.  This  derivation  was  written  in  Japanese  and  went  unnoticed  for  many  years.  Previous  to 
Moriya’s  work,  Kawada^  had  derived  the  potential  for  an  infinitely  long  vortex  filament  which  Lerbs^ 
used  to  develop  charts  useful  in  design.  Since  the  work  of  Lerbs,  only  minor  refinements  have  taken  place 
in  lifting-line  theory.  A  derivation  of  the  lifting-line  equations  and  explanation  of  the  numerical  techniques 
used  in  their  evaluation  are  given  by  Morgan  and  Wrench.^ 

Even  while  development  of  the  lifting-line  theory  was  in  progress,  it  was  recognized  that  marine 
propellers  were  too  broad  for  direct  application  of  the  lifting-line  results.  Efforts  were  devoted  to  finding 
correction  factors  which  would  allow  use  of  the  lifting-line  results  by  modifying  the  propeller  design  so  that 
they  would  actually  deliver  the  required  thrust.  The  calculation  of  such  correction  factors  was  quite 
tedious,  and  a  considerable  amount  of  simplification  was  incorporated.  Ludwieg  and  Ginzel^  derived  an 
approximate  correction  factor  for  wide-bladed  propellers  which  was  extensively  used  before  high-speed 
computers  made  possible  more  accurate  calculations  for  the  individual  design.  Strscheletzky  ^  presented  a 

1  n 

basis  and  some  numerical  results  for  a  general  lifting-surface  theory;  however,  it  was  not  until  the  use  of 
computers  became  popular  in  the  early  1960’s  that  lifting-surface  calculations  became  practical.  Typical  of 


12Kawada,  S.,  “Induced  Velocity  by  Helical  Vortices,”  Journal  of  Aeronautical  Sciences,  Vol.  3,  pp.  86- 
87  (1936). 

^Lerbs,  H.W.,  “Moderately  Loaded  Propellers  with  a  Finite  Number  of  Blades  and  an  Arbitrary  Distri¬ 
bution  of  Circulation,”  Society  of  Naval  Architects  and  Marine  Engineers  Transactions,  Vol.  60,  pp.  73-123 
(1952). 

^Morgan,  W.B.  and  J.W.  Wrench,  “Some  Computational  Aspects  of  Propeller  Design,”  Methods  of 
Computational  Physicis,  No.  4,  pp.  301-331,  Academic  Press,  Inc.,  New  York  (1965). 

15Ludwieg,  H.  and  I.  Ginzel,  “Zur  Theorie  der  Breitblattschraube ,”  Aerodynamische  Versuchsanstalt, 
Gottingen,  Report  44/A/08  (1944);  see  Ginzel,  G.I.,  “Theory  of  the  Broad-Bladed  Propeller  ”  Aeronautical 
Research  Council,  Current  Papers  208  (1955). 

16Strscheletzky,  M.,  “Hydrodynamische  Grundlagen  zur  Berechnung  der  Schiffsschrauben ,”  Verlag  G. 
Braun,  Karlsruhe,  West  Germany  (1950). 

^Strscheletzky,  M.,  “Berechnungskurven  fur  dreiflugelige  Schiffsschrauben,”  Verlag  G.  Braun,  Karlsruhe, 
West  Germany  (1955). 
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the  lifting-surface  investigations  are  the  analyses  presented  by  Yamazaki  Sparenberg,  Kerwin/  ’ 

no  nq  ^f)  o  t  on  on 

Kerwin  and  Leopold,  ’  y  Pien,  Pien  and  Strom-Tejsen,  Hanaoka,  and  Murray.  Reviews  of  many 


^Yamazaki,  R.,  “A  Study  on  Screw  Propellers,”  Memoirs  of  the  Faculty  of  Engineering,  Kyushu 
University,  Japan,  Vol.  19,  No.  1,  pp.  1-75  (1960). 

^Yamazaki,  R.?  “On  the  Theory  of  Screw  Propellers,”  Fourth  Symposium  on  Naval  Hydrodynamics 
(1962), 

^^Yamazaki,  R.,  “On  the  Theory  of  Screw  Propellers,”  Memoirs  of  the  Faculty  of  Engineering,  Kyushu 
University,  Japan,  Vol.  23,  No.  2,  pp.  97-112  (1963);  corrected  version  of  1962  paper. 
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1 Yamazaki,  R,,  “On  the  Theory  of  Screw  Propellers  in  Non-Uniform  Flows,”  Memoirs  of  the  Faculty  of 
Engineering,  Kyushu  University,  Japan,  Vol.  25,  No.  2,  pp.  107-174  (1966). 

77 

Yamazaki,  R.,  “On  the  Propulsion  of  Ships  in  Still  Water  (Introduction),”  Memoirs  of  the  Faculty  of 
Engineering,  Kyushu  University,  Japan,  Vol.  27,  No.  4,  pp.  187-220  (1968). 

77 

Yamazaki,  R.,  “Theory  of  Unsteady  Propeller  Forces,”  Seventh  Symposium  on  Naval  Hydrodynamics, 
Rome  (1968). 

^Yamazaki,  R.,  “On  the  Theory  of  Unsteady  Propeller  Forces,”  Memoirs  of  the  Faculty  of  Engineering, 
Kyushu  University,  Japan,  Vol.  28,  No.  3,  pp.  157-206  (1969). 

75 

Sparenberg,  J.A.,  “Application  of  Lifting  Surface  Theory  to  Ship  Screws,”  International  Shipbuilding 
Progress,  Vol.  7,  No.  67,  pp.  99-106  (1960). 

76 

^Kerwin,  J.E.,  “The  Solution  of  Propeller  Lifting  Surface  Problems  by  Vortex  Lattice  Methods,” 
Massachusetts  Institute  of  Technology,  Naval  Architecture  Department  Report  (June  1961). 

77 

* Kerwin,  J.E.,  “Linearized  Theory  for  Propellers  in  Steady  Flow,”  Massachusetts  Institute  of  Technology, 
Naval  Architecture  Department  Report  (1963). 

no 

^°Kerwin,  J.E.  and  R.  Leopold,  “Propeller  Incidence  Correction  Due  to  Blade  Thickness,”  Journal  of 
Ship  Research,  Vol.  7,  No.  2,  pp.  1-6  (1963). 

on 

^Kerwin,  J.E.  and  R.  Leopold,  “A  Design  Theory  for  Subcavitating  Propellers,”  Society  of  Naval 
Architects  and  Marine  Engineers  Transactions,  Vol.  72,  pp.  294-335  (1964). 

^^Pien,  P.C.,  “The  Calculation  of  Marine  Propellers  Based  on  Lifting  Surface  Theory,”  Journal  of  Ship 
Research,  Vol.  5,  No.  2,  pp.  1-14  (1961). 

7 1 

Pien,  P.C.  and  J.  Strom-Tejsen,  “A  General  Theory  for  Marine  Propellers,”  Seventh  Symposium  on 
Naval  Hydrodynamics,  Rome  (1968). 

72 

Hanaoka,  T.,  “Hydrodynamics  of  an  Oscillating  Screw  Propeller,”  Fourth  Symposium  on  Naval 
Hydrodynamics  (1962). 

77 

Murray,  M.T.,  “Propeller  Design  and  Analysis  by  Lifting  Surface  Theory,”  International  Shipbuilding 
Progress,  Vol.  14,  No.  160,.  pp.  433-451  (1967). 
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of  these  investigations  can  be  found  in  the  paper  by  Wu  and  the  books  of  Isay.  ’  Cox J ;  has  surveyed 
the  state  of  the  art  for  subcavitating  propeller  design  and  concludes  that  for  propellers  described  as  lightly 
to  moderately  loaded,  experimental  results  for  propeller  performance  agree  with  the  requirements  set  in  the 
design  procedure.  (Since  the  design  procedures  are  based  on  inviscid  flow,  viscous  corrections  must  be  made. 
Presently  this  is  done  on  an  empirical  basis.)  Cox  points  out  that  research  is  still  needed  in  steady-flow 
design  theory  to  extend  the  theory  to  heavily  loaded  propellers.  The  present  work  contributes  to  the  heavily 
loaded  design  problem  since  the  lifting-surface  theory  of  Chapter  3  is  carried  out  to  second  order. 

While  these  design  theories  were  being  developed,  the  problem  of  predicting  the  performance  was 
also  being  investigated.  Kerwin  calculated  the  performance  of  a  series  of  propellers  over  an  operating 
range,  using  approximate  curved-flow  corrections  derived  for  use  in  design.  Differences  of  10  percent,  com¬ 
pared  with  experimental  results,  were  found,  although  occasionally  excellent  agreement  was  obtained. 

Analyses  of  the  lifting-surface  formulations  to  predict  performance  from  the  geometry  have  been 
attempted  since  Kerwin’s  investigation.  To  date  most  of  these  investigations  have  not  compared  predictions 
with  experimental  results.  Yamazaki^^  made  some  comparisons  with  experiments,  and  he  concluded 
good  agreement  when  the  pitch-to-diameter  ratio  was  small.  Murray  J  calculated  the  performance  of  pro¬ 
pellers  he  had  previously  designed  and  found  discrepancies  between  the  calculations.  For  one  propeller  the 
predicted  and  design  value  of  the  thrust  differed  by  8  percent.  Isay  and  Armonat  calculated 


*^Wu,  T.Y.,  “Some  Recent  Developments  in  Propeller  Theory,”  Schiffstechnik,  Vol.  12,  No.  60,  pp.  1-11 
(1965). 

^Isay,  W.H.,  “Propellertheorie,  Hydrodynamische  Probleme,”  Springer-Verlay,  Berlin  (1964). 

*^Isay,  W.H.,  “Moderne  Probleme  der  Propellertheorie’’ Springer-Verlag,  Berlin  (1970). 

37 

J  Cox,  G.G.,  “State-of-the-Art  for  Subcavitating  Propeller  Design  Methods,”  Appendix  II  of  the  Report  of 
the  Propeller  Committee,  12th  International  Towing  Tank  Conference,  Rome  (1969). 

30 

Kerwin,  J.E.,  “Machine  Computation  of  Marine  Propeller  Characteristics,”  International  Shipbuilding 
Progress,  Vol.  6,  No.  60,  pp.  343-354  (1959). 

39 

Isay,  W.H.  and  R.  Armonat,  “Zur  Berechnung  der  Potentialtheoretischen  Druckverteilung  am 
Flugelblatt  eines  Propellers,”  Schiffstechnik,  Vol.  13,  No.  67,  pp.  75-89  (1966). 

Lerbs1,5  defines  loading  ranges  based  on  the  importance  of  the  trailing- vortex  position  in  the  calcu¬ 
lations  of  induced  velocities  at  the  lifting  line. 
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40 

bladc-prcssurc  distribution  but  made  no  comparisons  with  experiments.  Sugai  made  extensive  calculations 
of  blade  pressure  distributions,  thrust  coefficient,  and  comparative  camberlinc  studies  but  did  not  compare 
his  results  with  experiments. 

None  of  the  performance  evaluation  studies  mentioned  previously  included  thickness  effects.  All 

1001 

were  linearized  formulations,  although  Yamazaki  formulated  a  theory  which  did  not  explicitly 
require  linearization  but  did  require  zero  thickness.  Later  papers  by  Yamazaki  included  thickness 
effects  but  were  based  on  linearized  formulations.  He  also  gave  equations  governing  the  position  of  the  shed 
vortex  sheet.  In  each  case  described,  evaluation  of  the  numerical  formulation  usually  required  further 
approximation.  The  linearization  and  approximations  point  out  the  complexity  of  the  problem  resulting 
from  the  geometry  of  the  propeller  and  its  shed  vortex  sheet. 

Other  problems  in  propeller  theory  which  are  being  investigated  concern  unsteady  effects  and 
viscous  effects.  Study  of  these  two  aspects  of  propeller  theory  is  still  in  its  infancy,  and  relatively  few 
papers  have  been  presented.  Most  of  the  literature  on  unsteady  effects  can  be  found  in  Yamazaki, 

Pien  and  Strom-Tejsen,^  and  Tsakonas  and  Jacobs.^  Several  references  concerning  boundary  layers 
applicable  to  propeller  blades  can  be  found  in  Armonat.^  Since  these  two  areas  are  not  covered  in  the 
present  work,  no  further  details  of  these  investigations  will  be  given. 

E— STATEMENT  OF  THE  PROBLEM 

It  is  appropriate  at  this  point  to  define  the  problem  to  be  considered;  the  flow  field  of  a  hubless 
propeller  with  Z  identical,  equally  spaced,  blades  operating  in  an  unbounded,  inviscid,  incompressible  fluid. 
The  rectilinear  and  angular  velocity  are  coaxial  and  constant.  The  exact  lifting-surface  theory  is  first 
formulated  without  further  restrictions;  however,  the  singular-perturbation  problem  is  further  restricted  in 
geometry.  Experience  with  wings,  Van  Dyke,  Thurber,  and  Rotta,  indicates  that  more  terms  can  be 
calculated  for  the  simplest  geometry,  i.e.,  without  sweep  and  dihedral.  Accordingly,  the  blade  should  be 
neither  raked  nor  warped  (these  terms  arc  defined  in  Chapter  1).  The  reason  for  this  restriction  is  explained 
in  Chapter  4. 

1  o 

The  usual  assumptions  (Lerbs)  about  the  significant  interaction  effects  of  importance  to  the 
position  of  the  shed  vortex  sheet  are  not  made  in  the  lifting-surface  analysis.  The  position  of  the  vortex 
sheets  can  be  taken  into  account  as  the  perturbation  solutions  evolve.  An  approximation  is  made  in  the 
lifting-line  analysis  which  corresponds  to  moderately  loaded  theory. 


^Sugai,  K.,  “Hydrodynamics  of  Screw  Propellers  Based  on  a  New  Lifting  Surface  Theory,”  Selected 
Papers  from  the  Journal  of  the  Society  of  Naval  Architects  of  Japan,  Vol.  4,  Tokyo,  pp.  96-106  (1970). 

^Tsakonas,  S.  and  W.R.  Jacobs,  “Propeller  Loading  Distributions,”  Journal  of  Ship  Research,  Vol.  13, 
No.  4,  pp.  237-257  (1969). 

^Armonat,  R.,  “Untcrsuchung  der  Druckverteilung  eines  Propellers  unter  Beriicksichtigung 
grenzschiehtbedingter  Massstabseffektc,”  Schiffstechnik,  Vol.  16,  No.  81,  pp.  41-54  (1969). 
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F— DESCRIPTION  OF  CONTENTS 

Because  the  geometry  plays  such  an  essential  role  in  propeller  theory,  the  first  chapter  is  devoted  to 
an  examination  of  the  geometrical  specification  of  a  propeller  blade.  The  normal  to  a  blade  is  determined, 
and  certain  differentials  are  examined.  The  findings  are  used  in  later  chapters. 

Since  no  derivation  of  the  exact  potential  function  associated  with  the  finite  dimensions  of  the 
propeller  could  be  found  in  the  published  literature,  an  expression  has  been  formulated  in  terms  of  boundary 
values  in  Chapter  2.  By  a  straightforward  linearization,  the  regular-perturbation  problem  has  been  obtained 
in  Chapter  3.  The  applications  to  both  design  and  performance  calculations  is  discussed.  The 
perturbation  solution  is  carried  out  to  second  order.  Although  the  performance  calculations  result  in  an 
integral  equation  in  the  regular-perturbation  problem,  in  Chapter  4  the  performance  problem  is  considered 
from  the  viewpoint  of  singular  perturbations,  and  a  quadrature  results.  Two  terms  in  the  series  have  also 
been  calculated  in  this  analysis. 

In  Appendix  A,  an  integral  expression  for  the  potential  of  a  lifting-line  has  been  derived  from  the 
Biot-Savart  Law.  (Although  this  expression  can  be  integrated  in  terms  of  known  functions,  the  nested 
infinite  series  which  result  are  not  thought  to  be  practical  for  computation.)  In  Appendix  B  the  intermediate 
expansion  of  the  outer  potential  is  obtained.  This  expansion  is  used  in  Chapter  4. 
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CHAPTER  1 

PROPELLER  GEOMETRY 


The  details  of  propeller  theory  are  considerably  more  complicated  than  wing  theory  because  of  the 
geometry.  In  this  chapter,  the  analytical  specification  of  the  blade  shape  is  developed,  and  certain  ex¬ 
pressions  needed  later  in  the  analysis  are  obtained.  First  the  coordinate  systems  are  described,  and  the 
relations  between  them  are  derived.  In  addition,  the  normal  to  the  blade  is  found.  In  another  section,  the 
blade  sections  are  discussed,  and  the  small  parameters  (thickness-  and  camber-to-chord  ratio  or  chord- 
diameter  ratio)  are  explicitly  considered  in  the  expressions  for  the  derivatives  appearing  in  the  normal.  In 
the  last  section,  the  gradients  of  functions  expressed  in  these  coordinate  systems  are  obtained. 

A— COORDINATE  SYSTEMS  AND  BLADE 
GEOMETRY 

To  describe  the  flow  field  quantitatively,  two  right-handed  coordinate  systems  are  used,  both 
attached  to  the  moving  propeller.  The  first  of  these  is  a  Cartesian  reference  frame  (xty,z)  with  z  pointing 
away  from  the  axis  of  rotation  along  a  reference  line  in  the  blade  and  with  *  pointing  along  the  axis  of 
rotation  with  positive  displacement  measured  downstream  from  the  propeller.  The  second  system  is  cy¬ 
lindrical  polar  coordinates  (C5,6txJ.  Radial  distances  are  measured  from  the  axis  of  rotation  in  a  (y,z)  plane, 
and  0  is  measured  from  the  z  axis  in  the  clockwise  direction  looking  along  the  positive  x  axis;  see  Figure  1. 


Unit  base  vectors  in  the  (x,y>z)  system  are  the  usual  (i£k)  and  base  vectors  in  the  (co,d,x)  system  are 
denoted  by  (e&eQ,i),  respectively,  where 
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and 


(M) 


err=  -  sin  0  j  +  cos  0  k 


£ 0  ~  -  cos  Oj_-  sin  0  k_ 

In  general,  blade-section  offsets  for  propellers  arc  measured  from  a  reference  line  lying  on  a  cy¬ 
linder  of  constant  radius.  This  viewpoint  is  adopted  for  some  of  the  analyses  (Chapters  2  and  3)  but  for 
another  portion  of  the  work  (Chapter  4)  a  slightly  different  viewpoint  is  necessary.  This  second  viewpoint 
is  that  the  geometry  is  defined  for  sections  cut  by  a  plane  at  constant  z  on  the  reference  blade  (the  reference 
blade  is  the  one  which  lies  as  close  as  possible  to  the  z-axis).  Although  a  set  of  unit  orthogonal  reference 
vectors  is  assigned  for  each  viewpoint,  neither  system  is  an  orthogonal  curvilinear  coordinate  system,  except 
for  a  special  case  in  the  second  system.  (The  orthogonal  coordinate  system  with  the  cylindrical  surface 
would  be  cylindrical  polar  coordinates.) 

For  the  sections  defined  on  a  right  circular  cylinder,  the  blade-section  reference  line  on  the  surface  of 
the  cylinder  is  called  the  geometric  pitch  line.  It  makes  an  angle  \pp(co)  with  the  y,z  plane.  The  geometric 
pitch  line  usually  lies  along  the  nose-tail  line,  which  is  a  line  on  the  cylindrical  surface  joining  the  leading 
point  and  the  trailing  point  of  the  blade  section  cut  by  the  cylinder.  Another  important  line  is  the  blade- 
reference  line  (also  called  stacking  line)  which  is  defined  as  the  line  through  an  approximate  datum  point  at 
each  radius.  The  datum  point  is  usually  on  the  nose-tail  line  and  generally  will  be  the  midchord  point.  The 

blade  is  raked  when  the  blade-reference  line  is  given  by  x  -  tii(to);  0  =  6^  ~  — —  (b  *  0,  7,...  Z  -  7 ). 

When  the  reference  line  is  given  by  d  =  6^  +  W( c5),  x  =  0,  the  blade  is  said  to  be  warped.  A  combination 
of  warp  and  rake  is  called  skew;  usually  the  displacement  is  along  a  helix  which  passes  through  a  straight 
line  in  the  plane  x  -  0.  On  the  cylindrical  surface  co  =  constant ,  a  coordinate  system  (%j,  ^  £3 J  is  con¬ 
structed  with  on  the  surface  of  the  cylinder  and  measured  from  the  blade-reference  line  along  the 
constant  pitch  of  the  nose-tail  line.  Positive  values  point  in  the  downstream  direction.  The  variable  ^ 
is  normal  to  the  axis  and  points  in  the  upstream  direction.  Unit  orthogonal  reference  vectors  (ej^.e ~) 
are  in  the  direction,  respectively,  but  the  coordinate  system  is  only  locally  orthogonal,  i.e.,  these 

are  not  the  unitary  base  vectors  described  by  Wills.^  The  system  is  left  handed  to  permit  specification  of 
the  blade-section  geometry  in  the  conventional  two-dimensional  orientation  (positive  camber  in  the  positive 
^2  direction).  A  schematic  of  the  blade-coordinate  system  is  shown  in  Figure  2.  Although  the  blade-section 
reference  line  is  shown  as  the  nose-tail  line,  it  is  not  necessary  in  the  following  development  and  will  not 
always  be  used  that  way. 


^Wills,  A.P.,  “Vector  Analysis  with  an  Introduction  to  Tensor  Analysis,”  Dover  Publications,  Inc.,  New 
York  (1958). 
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x  =  M(u})  +  £/  sin  v?  (u>)  -  %2  cos  Vp  (“>) 


y  -  -  to  sin  \oh  +  W(u>)  + 


£/  cas  'Pn  +  %2  sin  % 


z  =  oj  cos  I  0  +  W(  00)  + 


0  cos  Vp  +  ?2  sin  Pp 


One  can  express  the  position  vector  in  these  variables  by  substituting  into  the  expression 

L  =  x  L  +  >'i  +  z  h. 

and  the  unit  base  vectors  in  the  £y  and  ^  direction  can  be  found: 


-1  =  TT~  =  s/,?  cos  ®  cos  VpJ  *  sin  9  cos  $  k_ 


sin  $pl  +  cos  <pp  eg 


£ 2  -  — =  -  cos  ^PpL'  cos  9  sin  \p^_  -  sin  6  sin  ifipk^ 

=  -  cos  v?pL  +  sin  {pp  e_  Q 

where  6  is  given  in  Equation  (1-2). 

The  unit  base  vector  in  the  Go  direction  is  the  unit  cylindrical  polar  vector  in  the  co  direction.  It  is 
given  by  Equation  (1-1): 

sin  9  +  k  cos  6 

(The  partial  of  r  with  respect  to  co  will  not  give  this  value  since  the  £j,  £2  ^ase  vect°rs  are  functions  of  co.) 
By  construction  Cy  x  -  ~JL£o 

The  inverse  of  Equations  (1-3)  is 


£y  =  cos  Kpp  yy^~z^ (tan  1  j  -  9^  -  W)  +  (x-di)  sin  \pp 

=  sin  <pp  ]/y2  +  z2  (tan1  (—V  8b-W)-  fx-SlJ  cos  yp 


co  =  Vy2  +  z 2 


Blade  sections  are  given  by 


=  E<ij,  u) 
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for 


c2  (<*>)<%]  <c}  (<Z) 
o  < Zo<R 

where  c^co)  is  the  distance  to  the  leading  edge, 

Cj(ZS)  is  the  distance  to  the  trailing  edge,  and 
R  is  the  propeller  radius. 

E  is  the  sum  of  a  single-valued  (camber)  and  a  double-valued  (thickness)  function  which  combine  to  describe 
the  offsets  from  the  axis  on  the  cylinder  of  radius  S3.  The  function  E  is  considered  in  detail  later  in  this 
chapter.  Figure  3  shows  the  profile  geometry. 


*1 


A  point  on  the  surface  of  a  blade  is  given  as  a  function  of  two  variables  (%j9  S3)  by 

£  =  f  <R  +  sin  iflp-  E  cos  Vp)± 

_  £7  cos  <£>  +  E  sin  <zl 

db  +  W+  1  p _ 

ST 

6b  +  w  + 

Two  tangent  vectors  on  the  blade  surface  are  given  by 
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ds 


dv 

~~  3dT 

sec.  e.g.,  Wills^  and  a  vector  normal  to  the  surface  (but  not  a  unit  normal)  is  given  by 


0-8) 


N  =  +i/Xr2 


(1-9) 


where  the  sign  is  positive  on  the  suction  or  back  side  and  negative  on  the  pressure  or  face  side. 

The  components  of  these  vectors  in  the  (ij,k)  system  are  lengthy  but  when  expressed  in  the 
(e sYstem  they  are  shorter.  In  particular,  the  normal  is 


N=± 


9ZT  J  /  9 E\  c% 

'15 


J6i  /  9 E  \ 

( _ dW  *1  cos  *p  +  f.  sin  /  b£  \  ) 

\  d<3*  z  )  \  *p  '  95; C0'  *p\ 


(1-10) 


For  E  =  0,  the  resulting  expression  for  the  normal  is 


£ 2 


(d<Pp  sin  ipp  cos  <f>p' 


doo  CO 


dR  _ dW 

+  ——cos  -  co —  sin 

p  dz>  p 


(i-n) 


which  is  the  normal  to  the  helical  surface  of  pitch  angle  \p  ( oj)  through  the  blade-reference  line.  Components 
in  the  (i,j,k)  directions  can  be  found  when  Equations  (1-4)  and  (1-1)  are  substituted  into  these  expressions. 

Previous  investigators  have  assumed  that  the  only  significant  contribution  to  was  e^-  Note  that 
even  for  blades  without  rake  or  warp,  an_£^  component  exists.  This  remaining  term  can  be  simplified  by 
expressing  in  terms  of  the  pitch,  Pfoo).  The  pitch  is  the  distance  of  axial  advancement  in  one  complete 

revolution  of  a  screw  of  pitch  angle  ip^ZS): 
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tan  <pp(  c o)  =  ■ 


Using  this  relationship,  one  finds 


%  sin  *p  cos  Vp  dP 


d^2tt  (of  +  P2/^) 


dP  cosz  i pp 
doo  2itoS 


(M2) 


Thus  for  helical  surfaces  of  constant  pitch,  this  component  of  the  normal  is  zero.  Since  few  modern  pro¬ 
pellers  have  constant  pitch,  this  term  must  be  retained. 

For  the  second  system  of  blade  specification  the  sections  are  cut  out  by  a  plane  at  a  given  z  value. 
As  will  be  explained  in  Chapter  4,  neither  rake  nor  warp  is  considered  with  this  viewpoint.  Hence  the 
blade-reference  lines  are  the  6 ^  values  and  for  the  reference  blade  along  the  z  axes,  the  nose  tail  line  is  on  a 
plane  and  forms  an  angle  ^z)  with  the  y-direction.  The  coordinate  system  is  similar  to  the  previous  one, 
except  that  ^  are  measured  in  an  x,y  plane  at  a  constant  z  value. 

The  Cartesian  components  are 


X  =  sin  Ifig  *  %2  cos 
y  =  -(%l  cos  <pB  +  \2  sin  <f>B) 


n-i3> 


Similar  to  the  previous  discussion,  the  unit  vectors  are 


ej=-~  =  sin  vBi  -  cos  <pB[ 

°Z] 


cos  ipgi  -  sin  i pBf 


(1-14) 


ur 

and  the  third  vector  is  the  k  vector  such  that  k  =  -  ej  X  e2  and  here  also  k  f 


The  blade  surface  is  here  assumed  to  be  given  by 


$2  =  e(Z 1>  zJ 


(1-15) 


where  E  is  a  sum  of  a  double-  and  a  single-valued  function  describing  the  thickness  and  camber  distributions, 
respectively. 
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As  before,  a  point  on  the  surface  is  found  by  substituting  this  expression  in  the  (’artesian  com¬ 
ponents  given  by  liquation  (1-13),  and  two  tangent  vectors  to  the  surface  are  found  by  differentiation.  The 
cross  product  of  these  two  tangent  vectors  gives  a  normal  vector, 


\  dz  dz  d%j  dz  J  ” 


This  is  similar  to  the  expression  for  the  normal  where  the  section  shape  was  measured  on  the  surface  of  a 
cylinder;  however,  the  terms  for  the  component  normal  to  the  template  plane  show  marked  differences. 


B— SECTION  GEOMETRY 


In  this  section,  the  specific  form  of  the  shape  function  is  examined.  In  general  this  shape  function 
can  be  written 


%2  ~  E($l  • 

where  is  either  the  cylindrical  radius  or  the  z  variable. 

The  blade  shape  is  made  up  of  a  single-valued  and  a  double-valued  function.  The  single-valued 
function  is  the  camberline,  and  the  double-valued  function  is  the  thickness.  These  functions  can  be  con¬ 
sidered  separately. 

First  the  thickness  can  be  represented  by 

h=  -  ft^1' 


7'T~V(cr  M“ 

V  1 


tj  \» 


Cj  *  and  C2  -  are  sh°wn  Figure  3,  and  in  general  afl  =  on(^).  At  the  leading  edge  =  c~>t 

and  at  the  trailing  edge  =  Cj. 

Often  the  same  form  is  maintained  spanwise  and  then  an>  where  t(^)  Is  the  thickness 

to  chord  ratio  and  the  an  are  now  constants. 

For  the  regular  perturbation  problem, 

rfy  =  eT($3)  (1-18) 

that  is,  the  thickness  is  small  relative  to  the  chord.  By  assumption  all  derivatives  remain  of  first  order 
although  by  inspection  this  breaks  down  at  the  leading  and  trailing  edge,  where  at  the  leading  edge 


o.q 


(1-19) 
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Even  though  this  singularity  is  not  compatible  with  the  0(e)  assumptions  in  the  derivation, 
solutions  obtained  by  ignoring  this  discrepancy  have  been  useful.  Several  techniques  for  improvement 
(Friedrichs*  and  Van  Dyke^)  are  known  but  they  usually  break  down  after  a  few  terms. 

In  the  singular-perturbation  problem,  the  chord  length  is  the  small  quantity,  and  in  this  case  it  is 
convenient  to  put 


ci  ~  *hi  (y 

c2  =  eh2^3> 


(1-20) 


where  e  can  be  considered  the  ratio  of  maximum  chordlength  to  diameter.  In  this  case,  the  equation  for 
thickness  becomes 


The  normal  vector  involves  derivatives  of  tliis  function  with  respect  to  both  and  First,  the  derivative 
with  respect  to  is 


dE< 


T  1/2  (h1+h2)-  j-  N 


~n 


‘  44)  (Kr 


=0  €(hl  h2> 


hj  -  h2 


N 

2  na~ 
n=l 


n  \*(*rh2). J 


n-1 


(1-22) 


In  Chapter  4,  the  form  of  this  expression  in  inner  or  stretched  variables  is  needed.  Let  the  variable  be  Oj 
=  £y/e,  then  converting  this  to  the  stretched  coordinate,  one  finds 


T  l/2(h]  +  h2)  -  0 


bE. 

(°i>  y 


1  N 

2 


^(hrOj)(Oj-h2)  n=0  n\hl~h2 
7  -  Oj)(arh2) 


N 

2 


hl~h2. 


(hj  -  h2)  n=l 

This  is  entirely  a  zero-order  term  with  respect  to  the  chordlength  parameter. 


(1-23) 
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Second,  l he  derivative  of  the  thickness  function  witli  respect  to  the  variable  is 


■:T  -(Cj  c2)'+$i  (Cj  +  c2)'  N 


Wmm~7=  =  2  M - 

3  y  (cr  tj) (%rc2>  11=0  \ci~c2 


%-c2t  2  ‘'At-tA" 

n-0  \C1  c2 


-  (°1  ~c2)'  N  /  £,  \ 

'(ertyrtrv  '  — 1  h  2  nan(—L-)n-1 


(C1  c2>  n=l  \C1  c2 


where  the  prime  denotes  derivative  with  respect  to  When  this  is  converted  to  the  inner  variable,  it  be- 


dEr  -  fhj  h^'+Ojfhj  +  h2f  N  /  °J 

(°ly  ~  2  / -  ^  an  f - - 

Y  (^1  “  ^ 2 )  n~®  \ll  ~  ^2 


+  EV<h  -o])(oI-h2)  2  a'n 


n=0  \hj  -  h2 


*V<hi  -  Oj)  (Oj  ■  h2) 


(hj-h/  N 


o,  2  na 


( hl~h2 ?  n=1  \hl~h2 


Hence,  this  expression  remains  of  order  e 

For  the  camberline  function,  the  appropriate  expression  is 

*2 = Ec  (b  y 


(crty(h-c2>  N  l  h  \ 

EJb  y = -  2  bn( - r  +y  (1'26) 

(Cj-c2)  n=0  \cj-c2J 

In  some  current  design  procedures,  the  shape  of  the  camberline  is  computed  at  each  radius,  and  in  general  the 
coefficients  in  the  summation  are  not  a  constant  times  the  same  function  of  ^  The  derivatives  of  this 
function  follow  the.  same  patterns  as  the  thickness  function.  For  the  regular-perturbation  problem,  they  are 
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both  the  same  order  and  for  the  singular-perturbation  problem,  the  ^  derivative  is  of  zero-order  and  the 
derivative  is  of  order  e.  For  most  camberline  shapes,  the  derivative  is  bounded. 

C— DIFFERENTIALS  EXPRESSED  IN  PROFILE  COORDINATES 

In  some  portions  of  later  chapters,  it  is  necessary  to  calculate  the  gradient  of  a  function  which  is 
given  in  the  profile  coordinate  system.  For  the  hypothetical  functions  Jj  =  J  j  (%j,  ^  £>)  and  ^2  ~  ^ 2  ^7 » 

^2*  z),  the  gradient  can  be  determined  by  straightforward  application  of  the  chain  rule.  Consider  the 
function  J  j  first. 

a ji  a j}  a j} 

vjj= - 1  + - 1  + - k  (1.27) 

dx  by  dz 

where 


20 


VJ . 


ay  2 
9*/ 


’■£/  + 


bh 

d%2 


—2  + 


bJ 2  dipfi  j  vj  2  vj  2 ' 


a  j?  a  j 

•[«/ - ?2“ 

a  z  dz  \  a  %2 


(1-30) 


This  expression  is  used  in  Chapter  4  as  well  as  one  for  the  Laplace  equation  expressed  in  the  pro¬ 
file  coordinates  (£y,  ^2<  z)-  From  repeated  application  of  the  chain  rule,  one  finds 


a 2j2  a 2j2  a 2j2  %2j2 

+  -  = -  + 


bx2  by2  d£,2  a^2 


a2y2  (d<pg\ 


dz “ 


l  dz 


a  'j. 


b2J, 


a  j- 


bJ-, 


o  5 7  0  "  s?  ‘  5/ 

3*/  9^  3^2  9|; 


2£;  $2! 


'<%\2  32/2  d2 

- - 


ay?  ay? 


c/2  /  a^a$2  <fc'  |_  9^2  9|;J 

/  32y 


a2/A  a2y 


+  2  (  %2 

dz  \  9  9£^  9zy 


2  \  u 
+ - 


9 


Obviously,  the  Laplace  equation  involves  the  sum  of  the  previous  two  equations.  It  has  been  broken  up  to 
call  attention  to  a  feature  of  the  expression  for  This  feature  is  that  any  scale  change  applied  to 

both  and  ^  cancels.  This  property  is  utilized  in  Chapter  4. 
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CHAPTER  2 

GENERAL  FORMULATION 


A—DESCRIPTION  OF  FLOW  FIELD* 

A  Z-bladed  propeller  of  finite  radius  R  is  assumed  to  be  rotating  with  constant  angular  velocity  12 
and  to  be  advancing  at  a  uniform  rate  U  into  an  unbounded,  inviscid,  incompressible  fluid;  however,  one 
for  which  the  Kutta  condition  is  satisfied  at  the  trailing  edge.  The  propeller  produces  a  thrust  which 
requires  that  the  average  pressure  on  one  side  of  the  blade,  called  either  the  back  or  suction  side,  is  lower 
than  on  the  other,  called  either  the  face  or  pressure  side.  The  Kutta  condition  requires  that  no  flow  cross 
the  trailing  edge.  The  general  flow  pattern  then  is  that  fluid  is  pushed  from  the  high-pressure  side  to  the 
low-pressure  side  around  only  the  leading  edge  and  the  tips.  The  tip  flow  deflects  the  streamlines  in  the 
vicinity  of  the  tips  toward  the  axis  of  rotation  on  the  low-pressure  side  and  away  from  the  axis  on  the  high- 
pressure  side.  This  flow  behavior,  together  with  the  requirement  of  continuity  in  the  pressure  field,  results 
in  a  discontinuity  in  direction,  but  not  magnitude,  of  the  fluid  velocity  leaving  the  blades  at  the  trailing 
edge.  This  layer  of  fluid  slip  is  called  a  vortex  sheet. 

The  flow  model  thus  consists  of  the  lifting  surface  and  the  vortex  sheet  advancing  and  rotating  into 
still  fluid. 

B— GOVERNING  EQUATIONS  AND  BOUNDARY 
CONDITIONS 

The  unsteady  flow  field  described  in  Section  A  can  be  made  steady  by  considering  the  flow 
relative  to  the  coordinate  systems  rotating  with  the  blade.  In  a  rotating-coordinate  system,  the  equations  of 
motion,  e.g,,  Kochin,  Kibel'  and  Roze^  Equation  2.7,8,  are 

q  .  grad  q  + — grad  p  «  -  212  X  £2  X  (12  X  r)  (2-1) 

“  P  ~ 


div  q  =  0 


(2-2) 


where  q(r) 

P 

P 

£2  -  -  £2/ 
r 


is  the  velocity  vector, 
is  the  pressure, 
is  the  density, 

is  the  angular  velocity  of  the  propeller  for  the  assumed  right-hand  rotation  shown  in  Figure  1, 
is  the  position  vector  of  a  point  in  the  flow  field. 


This  discussion  is  patterned  after  that  given  by  Prandtl  for  planar  wings. 

^Kochin,  N.E.,  et  al.,  “Theoretical  Hydromechanics,”  (Translation  of  Fifth  Russian  Edition),  Interscience 
Publications,  Inc.,  New  York  (1964). 
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The  flow  field  in  the  inertial  reference  frame  is  irrotational  everywhere  except  on  the  boundaries. 
In  the  rotating  coordinate  system  the  vorticity  is  given  by 

curl  q  =  -  2S7  =  2  (2-3) 

The  boundary  conditions  are  that  the  upstream  velocity  is  specified;  there  is  no  flow  through  the 
blade;  the  vortex  sheet  moves  with  the  fluid,  i.e.,  the  normal  velocity  at  the  sheet  is  zero;  and  the  velocity 
at  the  trailing  edge  is  finite  (Kutta  condition).  If  %  represents  the  pressure  side  of  the  blades,  Sg+  the 
suction  side,  Sy+  the  continuation  of  Sg+  onto  the  vortex  sheet,  Sy-  the  other  side  of  the  vortex  sheet, 
and  T(z)  the  trailing  edge,  then  the  boundary  conditions  become 

q+  q0  =  U±+  £2cj£g  as  *  *  -  00  (24) 

1  •£“  Oforr^eSb  (2-5) 

q  <+00  f°r a*  T  (2-6) 

where  U  is  the  advance  velocity  in  the  inertial  reference  frame 

n  is  the  normal  pointing  into  the  fluid 

Sfr  =  S  ^  S  \r  and 

b  denotes  one  of  the  Z  blades,  b  =  0,1 ,  ...  Z  -  1 

* 

If  the  velocity  is  taken  as 


q=  V<f>+  VxB  (2-7) 

where  B  is  the  vector  potential,  and  (p  is  the  scalar  potential,  then  Equation  (2-3)  gives 

VX(VX  B)=2 ft/ 


and  Equation  (2-2)  gives 

If  the  boundary  condition  at  upstream  infinity  is  associated  with  i?,  an  expression  for  B_  can  be  immediately 
stated: 

B  -  oS(  Six  -  0U)e_£  (2-8) 

%  $ 

where  U  must  be  a  constant  for  this  expression  to  satisfy  Equations  (2-3)  and  (2-4),  and  VXB  =  Qq- 
$ 

The  symbolic  vector  notation  is  used  for  the  differential  operators  in  curvilinear  coordinates. 

The  velocity  component  U  is  commonly  taken  as  a  function  of  radius  to  approximate  nonuniform  in¬ 
flow.  One  would  expect  such  shear  flows  to  be  axially  variable  also.  Shear  flows  imply  vorticity,  and 
careful  consideration  should  be  given  to  the  propagation  of  vorticity  from  upstream  infinity  to  the  propeller 
and  on  downstream. 
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The  remaining  unknown  is  now  only  the  scalar  function  0,  and  the  problem  can  be  reformu¬ 
lated  in  terms  of  it.  The  boundary  conditions  are: 


F0  +  0  as  x  t  £o  +  oo 

(2-9) 

V0  .  n  -  -  qQ  .  n  for  £€  Sb 

(2-10) 

|  v0|  <  •/*  oo  for  ref 

(2-11) 

In  addition  the  pressure  is  continuous  across  the  vortex  sheet. 

The  three  simultaneous  equations  of  motion  described  by  Equation  (2-1)  have  been  replaced  by  a 
single  equation  for  the  scalar  potential  0  which  depends  only  on  the  kinematics  of  the  flow.  The  dynamics 
of  the  flow  enter  only  in  the  determination  of  pressure  which  can  be  obtained  by  integrating  Equation  (2-1). 
The  integration  is  performed  by  using  the  expansion  formula 


q  .  grad  q  =  l/2V(q  .  q)  -  q  X  curl  q 
=  l/2V(q.q)-2tt  Xq 

and  by  noting  that  for  £2  ~  -  Sli 


Q.X(£lXr)  =  .to2  fie  — 

—  —  —  “CO 

-fi2  ~2 

= - For 

2 

Hence,  Equation  (2-1)  can  be  written  as 

1  Pi  O  _  9 

Vfcq  .q+---Sl2  c o2)  =  0 

2-  —  p  2 

Since  the  expression  is  zero  for  all  spatial  derivatives,  it  must  be  constant.  Thus 


P 

,q+— 
-  P 


where  K  is  a  constant  everywhere  in  the  fluid. 

Hence  the  pressure  is 

The  condition  that  the  pressure  be  continuous  across  the  votex  sheet  can  be  expressed  as 


p  =  j  |aT  +  Q?  oj2  -  q  .  q 


(l-2.)±esv+  (I'&reSy. 

or 


(2-12) 


(2-13) 


(2-14) 
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(2-15) 


(2q()+  H)  ■  H\~ e,<V+=  0 

LeSy 

In  the  next  section  an  expression  for  <p  is  found  which  satisfies  the  boundary  conditions  in 
Equations  (2-9)  through  (2-1 1)  and  Equation  (2-15). 

C— SOLUTION  IN  TERMS  OF  BOUNDARY  VALUES 

The  solution  of  Poisson’s  or  Laplace’s  equation  in  a  three-dimensional  volume  can  be  expressed 
using  Green’s  second  identity;  see,  e.g.,  Tychonov  and  Samarskii 

,11,  rr  11 
/jf/^77,  “  777  V  <t>)  •  v  — -  —  a  .  w  as  (2-16) 

v  Id  Id  Id  Id 

2 

where  v  is  the  total  fluid  volume, 

2  is  the  surface  which  bounds  the  volume,  including  the  surface  at  |r|  «♦ 00 , 

/  is  the  position  vector,  measured  from  a  fixed  point  in  the  fluid, 
n  is  the  unit  normal,  pointing  from  the  bounding  surface  into  the  fluid. 

If  0  satisfies  Laplace’s  equation,  the  volume  integral  is  zero  everywhere  in  the  fluid,  except  at  the  point  /  = 
0.  This  point  can  be  excluded  from  the  volume  by  enclosing  it  in  a  sphere.  Taking  the  limit  of  the  sur- 
face  integral  over  the  sphere  as  the  radius  goes  to  zero,  one  finds  the  following  standard  formula 

(2-,7) 

where  r  is  an  arbitrary  point  in  the  fluid,  and 
£  is  a  point  on  the  boundary. 


^Tychonov,  A.N.  and  A.A.  Samarski,  “Partial  Differential  Equations  of  Mathematical  Physics,”  Vol.  1, 
Holden-Day,  Inc.,  San  Francisco,  Calif  (1964). 

*Tychonov  and  Samarski, for  example,  show  that  other  shapes  excluding  the  singularity  are  also 
satisfactory  as  long  as  the  maximum  dimension  of  the  excluded  region  goes  to  zero. 
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It  is  necessary  to  know  the  value  of  0  on  the  boundary  as  well  as  in  the  flow  field.  That  is,  one 

seeks 

lim  <p(r)  =  lim  +  r\e)  =  <fr  (rQ) 
l+lo  V  +  o 

where  £  is  a  point  in  the  flow  field 
is  a  point  on  the  surface 

£  is  an  arbitrary  vector  pointing  into  the  fluid  fromj;0  and  the 
±  sign  is  used  when  rQ  is  on  the  ±  surface 

Two  situations  occur:  one  when^  is  a  point  on  the  boundary  of  a  solid  of  finite  thickness,  and  one  when 

j r  is  a  point  on  a  boundary,  both  sides  of  which  are  in  the  flow  field,  across  which  0  is  discontinuous.  In 

the  second  case,  the  limitj;  -+ •  produces  singularities  in  both  sides  of  the  surface. 

$ 

For  the  first  case,  for  which  rQ  is  on  a  single-sided  surface,  one  finds 


+  1  tf(  -rc4  7 

W 7* ■  'M-  * 


(2-1 8a) 


and  for  the  second  case,  for  whichj-  is  on  a  double-sided  surface,  one  finds 


^ + *  ffr  •  fir  ^ r-H 


(2-1 8b) 


where  the  improper  surface  integral  is  obtained  by  evaluating  the  integral  over  the  bounding  surface,  ex¬ 
cluding  the  region  surrounding  the  pointy,  and  then  by  taking  the  limit  as  the  maximum  dimension  of  the 
excluded  region  tends  to  zero. 

Since  the  present  problem  is  a  second  (or  Neumann)  boundary-value  problem,  it  is  important  to 
determine  the  values  of  the  gradient  of  0  on  the  boundary.  To  find  this  limit,  the  pointy  is  excluded  from 
the  surface  S  by  surrounding  it  with  different  shapes  s(£Qy  A),  where  X  characterizes  the  dimensions  of  the 
region.  With  the  same  notation  as  used  previously,  one  seeks 

lim  V<p(r)  =  lim  V0  (r  +r\e)=( 

L*Ln~  r\  -  o 


For  a  point  on  the  boundary  of  a  solid  of  finite  thickness 


-  lim\+0 


-j-(n  .  ftp)  ds 


2 '  s(ln>  V 


(2-19a) 


The  ±  superscripts  here  denote  different  sides  of  the  surface  as  well  as  the  limit  operation  when 
appropriate. 
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/ 

+  —  lim 
27T 

\+o 


gj  JL.3(!j^L± 

m  VJr^  I roSp 

^rS(ro,\) 


(fo-DY^s 


6  n 
Y()-o 

X 


(2- 1 9a) 
Cont’d 


where  n()  =  n(r(}),  and  <j>()  =  (pfr()). 

Wlicn  s(r()f\)  is  the  region  of  2  interior  to  a  right  circular  cylinder  of  radius  X  with  axis  then  K$  =  2tt; 
when  s(r(),\)  is  the  region  of  2  interior  to  a  square  of  half-side  X  with  axis  nQ,  then  K$  =  4  V^7~when 
s(r(r\)  is  the  region  of  2  interior  to  a  rectangle,  two  sides  of  which  are  distant  X  from_r^  and  the  other  two 
sides  of  which  are  an  arbitrary  finite  distance  from_r0,  then  Ks  =  4.  We  call  the  first  integral  in  Equation 
(2-19a)  a  Cauchy  principal-value  integral. 

For  a  point  on  a  boundary  in  the  flow  field 


(V<t>)+  +  (V(t>)  -  —  lim 

27r  \-o 


-Vi 


(n  .  V0y  (Is 


fin  3(lo'Z)  ■’! 

0  I  - - (l0 

-  A^\J^s(r0,X)  VM-*  Lv^ 


in  ds 


n0+  <P0+  +  Ho  <t>o 


(2-1 9b) 


The  limit  operation  indicated  in  Equations  (2-19)  is  considered  to  define  a  singularity  in  the 
mathematical  literature  (Tychonov  and  Samarski)^  while  in  the  aeronautical  literature  (Mangier,^  Ashley 
and  Landahl)^  such  forms  are  accepted.  In  Chapter  5,  a  form  of  the  expression  for  velocity  is  obtained 
for  which  only  Cauchy  principal  values  are  needed  so  that  numerical  analysis  can  be  performed  on  the  more 
convenient  form. 

The  values  of  the  gradient  on  the  boundary  are  generally  to  be  interpreted  as  VF(r) L  r  =  (VF)~' 

This  expression  is  not  equal  to  VF(rQ)  =  V(F±)  since  the  normal  components  are  not  equal.  Generally  the 
normal  component  is  given  by  a  separate  boundary  condition,  and  one  need  not  explicitly  consider  the 
distinction.  Throughout  the  text  the  notation  for  differentials  of  functions  known  on  the  surface  V(F~)  is 
indicated,  rather  than  the  field  value  evaluated  on  the  body  ( VF)~.  Although  this  notation  sometimes 
indicates  a  ridiculous  result  if  interpreted  literally,  it  is  used  since  it  permits  a  shorthand-like  notation  to  be 
developed. 


^Mangier,  K.W.,  “Improper  Integrals  in  Theoretical  Aerodynamics,”  Aeronautical  Research  Council, 
Current  Papers  94  (1952). 

^ Ashley,  H.  and  M.  Landahl,  “Aerodynamics  of  Wings  and  Bodies,”  Addison-Wesley  Publishing  Co., 
Inc.,  Reading,  Mass.  (1965). 
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The  surface  2  consists  of  the  Z  bounding  surfaces  of  sB^sv  (where  SB  is  the  blade  surface,  and 
the  shed  vortex  surface)  and  the  surface  for  |r|  ■+  At  this  point  it  is  necessary  to  suppose  that  the 
integral  over  the  bounding  surface  at  infinity  gives  at  most  a  constant.  (This  supposition  will  be  proven.) 
Hence,  to  within  this  additive  constant,  the  potential  for  the  propeller  can  be  written 


W-f  2  1 

4ir  b=0  JJ 


,  1  1 
<p(s)nfs)  .  V  — j-  — 

N  til 


n(s)  .  V  <Ms)  >  dS 


(2-20) 


where  S ^  is  one  of  the  M  bounding  surfaces 
Sb  =  SV+  w  SB+  ”  SB-  ^  SV- 

The  boundary  condition  in  Equation  (2-10)  can  be  used  to  simplify  Equation  (2-20)  since  the 
second  part  of  the  integral  becomes 


dS 


(2-21) 


The  reduction  in  area  in  Equation  (2-21)  follows  because  is  continuous  at  S,  and  n  (on  Sy+)  =  -n  (on  5^-). 
The  normal  on  is  directed  oppositely  to  that  on  Sy-,  and  hence  Equation  (2-20)  can  be  written 


(2-22) 
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The  circulation  in  a  fluid  is  defined  as  the  integral 


(2-23) 


where  is  the  vector  arc  length  along  the  closed  curve  which  completely  encircles  the  shed  vorticity. 
Since  c[q  is  a  continuous  point  function,  it  contributes  nothing  to  the  integral  in  Equation  (2-23).  For  a 
curve  which  lies  entirely  within  the  irrotational  flow  except  at  the  one  point  at  which  it  cuts  the  shed 
vortex  sheet 


r(%)~  k'^sv  (2_24) 

Further  simplification  is  possible  by  considering  the  continuous  vortex  sheet  to  be  composed  of  lines  of 
constant  circulation  which  leave  the  trailing  edge  of  the  blade  at  the  point  T(z). 

In  Figure  4,  from  curve  C2,  the  circulation  about  the  blade  at  T(z)  is  given  by 

T(T(z))  =  0+  (T(z))  -  <t>(T(z))  (2-25) 


z 


Figure  4  —  Schematic  of  Vortex  Sheet 
for  Determining  Circulation 


$ 

The  word  circulation  will  be  associated  with  the  scalar  jump  in  potential  across  a  surface  as  in 
Equation  (2-24),  and  the  word  vorticity  will  be  associated  with  the  vector  jump  in  velocity  across  a  surface. 

If  A  is  the  vorticity  vector  then  A.  =  n  X(<]+  *■  ff)  ~  H  X  V  (<^~  -  0‘ )  =  n  X  VT,  where  n  is  the  unit  normal; 
see  Reference  48. 

48 

Milne-Thomson,  L.M.,  “Theoretical  Hydrodynamics,”  Fifth  Edition,  The  Macmillan  Company,  New 
York  (1968). 
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Now  the  circulation  atr  e  Sv  is  determined  from  the  curve  Cj.  This  curve  is  connected  to  by  the  curves 
and  lying  along  a  line  of  constant  circulation*  Hence,  the  total  curve  is  simply  connected  and 


V(t>.d!  +  T(r  e  Sv>  -  T(T(z)) 


V(4*-<t>)  .d/+r(reSv)-r(T(z)) 


(2-26) 


By  definition,  Cj  lies  along  a  curve  of  constant  circulation  and  from  Equation  (2-24),  the  integrand  in 
Equation  (2-26)  is 

ar 

vr.t,=—=0  (2-27) 

dtl 

where  ij  is  a  unit  vector  along  the  lines  of  constant  circulation.  Hence 

T(reSv)  =  T(T(z))  (2-28) 

Conditions  are  sought  which  describe  the  position  vector  of  this  line  of  constant  circulation  which  leaves 
the  trailing  edge  from  the  point  Tj  e  T(z).  To  emphasize  that  this  point  is  located  at  the  blade,  let  it  be 
written  Zq.  The  position  vector  of  the  line  of  constant  circulation  can  be  described  as 

i=£(z0’  d;  ebf 

=  fj(z0,  6;  6b)erfd)+f2(z0,  6;  6b)L  (2-29) 


where  z 


0 


6 


h- 


2jrb 

Z 


is  the  parameter  describing  the  starting  position  of  the  line, 
is  the  independent  variable,  and 

is  the  parameter  describing  the  individual  vortex  sheets. 


The  functions  f](zg,  d;  6^)  and  0;  6^)  are  the  radius  and  axial  position  of  the  line,  respectively.  The 

parameter  Zq  ranges  from  zero  to  the  propeller  radius,  and  the  variable  6  ranges  from  the  trailing  edge 
®Te/z())  t0  °°*  ^ *0  ”  z0’  ^en  no  contracti°n  °f  the  shed  vorticity  takes  place. 

sjcsjc 

To  establish  differential  equations  governing  the  position  of  the  line,  Equations  (2-10)  and  (2-15) 
will  be  used  together  with  the  definition  that  the  strength  of  the  circulation  is  constant  along  the  line.  Since 
the  rate  of  change  of  F  is  zero  along  the  line, 


Alternatively,  one  could  take  f  -  $(zq,  x;  9 ^).  However,  the  usual  representation  given  in  Equation  (2-29) 
was  more  convenient. 

**  4.0  7 

The  derivation  here  parallels  that  given  by  Ciolkowski  and  Thurber. 

^Ciolkowski,  S.I.,  “The  Swept-Back  Wing,”  Ph.D.  Thesis,  New  York  University,  University  Microfilms 
24,694  (1955). 
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(2-30) 


vr k)  .-r=  o 

(16 

Since  the  vortex  sheet  is  actually  continuous,  the  function  6J  may  be  considered  a  function 
of  two  independent  variables  which  describe  the  surface.  Lines  of  constant  circulation  are  obtained  by 
holding  Zq  constant.  Tangent  vectors  in  the  surface  Sy  are  found  by  constructing; 

ar 


11  e 

ar 


—2 


dz 


0 


see,  e.g.,  Wills 

A  vector  normal  to  the  sheet  Sy  (but  not  a  unit  normal)  is  given  by 

N=  — 

*  bd  dzg 

which  is  assumed  to  be  non-zero. 

Hence,  Equation  (2-10)  can  be  expressed 


(q0+  V0)  .N 

Since  T  =  -  <£-|,  Equation  (2-34)  gives  two  equations: 


N=  0 


r  e S, 


vr  .  N  =  0 


N  .  U  +  vj<^  +  0-j 

Let  \p  =  /i 0+  +  <p~  \  .  Then  Equation  (2-36)  is 
LCSy 

N  .  (2q0  +Vijj)  =  0 

Now  Equation  (2-15)  can  be  written  as 


(2-31) 


(2-32) 


(2-33) 


(2-34) 


(2-35) 

(2-36) 


(2-37) 


Vr  .  (qQ  +  1/2 V\jj)  =  0  (2-38) 

It  would  appear  that,  since  only  two  unknowns, /y  and  are  sought,  the  scalar  Equations  (2-27)  or  (2-30) 
and  (2-37)  would  be  sufficient  to  establish  differentials  of  the  trajectory.  However,  derivatives  with  respect 
to  Zq  enter  into  the  equations  and  are  not  desirable.  For  a  more  elegant  form  of  the  differential  equation 
for  Equation  (2-35)  must  be  used. 
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liquations  (2-35)  and  (2-30)  show  that 


vrx(I]xN])  =  (vr -(vr  ■ii>n  =  o 

Since  Equation  (2-39)  is  zero,  any  multiple  of  it  is  also  zero.  Hence, 


(2-39) 


(2q()  +  Vip)  X  (  v  T  X  (il  X  N))  =  0 


(2-40) 


In  another  form  this  equation  is 


(2q()  +  Vi p)  .  (J.j  \  N)  J  Vr-  (2q()  +  Vp)  .  VfJ  ljXN=0 


Equation  (2-38)  simplifies  this  equation  to  be 


(2q0+vt)  -(i/xyv)}  vr 


This  equation  holds  everywhere  on  the  sheet.  Since 
Vr  y  everywhere,  it  requires  that 


(2q0+VM  .(tjxWa0 


(241) 


Thus  for  the  orthogonal  base  vector  system  JV,  N)}  the  average  velocity  at  the  vortex 
sheet,  (qg  +  l/2Vty),  has  no  N  component  from  Equation  (2-37)  and  no  t_j  XjV  component  from 
Equation  (241).  Hence  the  average  velocity  is  parallel  to  the  curves  of  constant  circulation,  and  the 
equation  for  f  can  be  constructed  from  two  of  the  three  scalar  equations  which  make  up 


The  vector_Ty  is 


and  on  the  shed  vortex  sheet 


(2q()  +  Vl//)XJy  =  0 
df}  bf2 


,2i0  +  _r  =i  S*  (2  at,  *2  „  *  ( 


3  n . 

2U+  dx  I1 


(242) 


Hence,  Equation  (2-42)  in  component  form  is 


f]  \2U 


2  ft/,  + 


1  3  A9^ 


fj  3  6  j  3  6 


(2-43) 


dip  (  dip)  9 
— r—  -12U  +  — - -  =0 


9  x  /  d  6 


(2-44) 
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Thus  from  Equation  (2-43) 


and  from  Equation  (2-45) 


(  /  3  0\  3  0 

»  rf'*7,T')-f>ra" 


3  0 

_U  +  1/2  bx 

it  }  d  4/ 

- -  _ r. 

2 fj  be 


1  9^ 

2  9a; 

be  "  _ 


(245) 


(246) 


1  3  \I ) 
~  30~ 


(2-47) 


Equations  (2-46)  and  (2-47)  are  nonlinear  coupled  integro-differential  equations  for  the  radial  and  axial 
position  of  the  lines  of  constant  circulation. 

The  integral  over  the  shed  vortex  sheet  in  the  expression  for  0,  Equation  (2-22),  can  be  simplified 
using  the  coordinates  Zq ,  6  used  to  describe  the  vortex  sheet;  see  Figure  4.  First,  in  the  expression 

j_  Lrl 

Vs\L-s\  | r-s? 

the  position  of  the  circulation  element  is  used  for  .s  so  the  integration  is  performed  for  constant  circulation 
strength.  Second,  the  value  of  the  vector  surface  area  element  n  dS  is  given  by 


ndS  =ljX±2  dz0  da 

b£  ar 

“  x  ~  dz  ^  da 

da  3  Zq  ” 


(2-48) 


where  a  is  a  dummy  variable  for  0;  see,  e.g.,  Wills.^  With  these  changes,  the  integral  over  a  shed  vortex 
sheet  becomes 


^  j 


where  b  denotes  an  individual  blade. 
Hence,  0  becomes: 


T(z0>dzc 


°o  /3| 

X-] 

Z-l 

\ba 

bz0 / 

'  Ir-il 

TE(zq) 

(249) 
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y  z-l  /y.  r-s  Z-/ 

~r  £  H  <t>(s)n  ~  dS  +  x  h,v(') 

4lT  h=0  S~  \r-sr  h=() 


1  t-J  rr  q0  -n 

+y~l  ft  - d* 

4nb=0JJSR  |r-s| 


(2-50) 


As  described  in  Chapter  1,  the  propeller  geometry  is  traditionally  defined  in  terms  of  blade 
sections  lying  on  cylinders  of  constant  radius.  The  nose-tail  line  lies  on  the  cylinder,  and,  hence,  blade  shape 
is  measured  from  a  helical  surface  of  radially  varying  pitch.  On  the  cylindrical  surface  oo-  constant,  a  left- 
handed  coordinate  system  (to  define  section  geometry  in  the  traditional  orientation)  is  constructed  with 
on  the  surface  of  the  cylinder  and  is  measured  positively  in  the  downstream  direction  along  the  constant 
pitch  of  the  nose-tail  line.  The  variable  is  normal  to  the  ^  axis  and  points  in  the  upstream  direction. 

Unit  orthogonal  vectors  (cjye^e^j  are  in  the  (%  jy  5 )  direction,  respectively,  but  the  coordinate 

system  is  only  locally  orthogonal;  this  coordinate  system  is  detailed  in  Chapter  1.  Although  the  pitch  is 
assumed  to  be  that  of  the  nose-tail  line,  it  is  not  necessary  in  the  following  development  and  will  not  be 
used  that  way  in  later  chapters. 

The  blade  shape  is  given  by 


%2  “  2?  I  c^ZS)  j  Kcj(oo),0<co<R 

=  ZS)  tEjtfj,  <5)  1 


(2-51) 


where  £*  is  the  camberline  function  and  Ej  is  the  thickness  function;  hence,  the  position  vector  of  a  point 
on  the  surface  is  given  by 

^  sin  \p  -  E  cos  Vp)i  ~  co  sin  (eb  +  w 


(2-52) 


cos  <PP  +  E  sin  v5,\  /  cos  <p  +  E  sin  \ 

- I  /  +  co  cos  \Qp  +  W+ - I  k 


where  Z5)  is  the  pitch  angle  of  the  blade-section  reference  line, 
(R(co)  is  the  rake,  and 
W(co)  is  the  warp. 

A  normal  is  given  by;  see  Chapter  1 


where  the  +  sign  refers  to  the  SB+  surface,  i.e.,  suction  or  back  side;  and  the  -  sign,  to  the  SB-  surface, 
i.e,,  pressure  or  face  side.  The  vector  element  of  area  is 

nds  dlj  dco  (2-54) 

These  expressions  can  be  substituted  into  Equation  (2-50)  to  give  precise  limits  on  the  integration.  Because 
of  nonlinearities,  no  conceptual  simplifications  result  unless  the  expressions  are  linearized. 

The  explicit  form  of  the  normal  allows  the  body  boundary  condition  to  be  simplified.  The  free- 
stream  velocity  can  be  resolved  to 

q0  =  ko\  cos  (fp  -  P)£j  +  tool  sin  (<pp  -  p)e2 
so  that  on  the  body,  the  boundary  condition  is 

N  "V<t>=  +  ]/lJ2  +  n2  ZS2  sin  (<pp  -  0)  -  ~  cos  (<pp  -  0)  j  (2-55) 

7  U 

where  j3  =  tan  - 

QZ5 


Equations  (249),  (2-50),  (246),  (247)  and  (2-55)  constitute  the  exact  formulation.  Equations 
(246)  and  (247)  are  nonlinear  in  0.  Although  Equation  (2-55)  is  linear  in  0,  an  integral  equation  results 
in  both  the  design  and  performance  applications.  Fortunately,  however,  the  integral  equation  can  be 
reduced  to  a  quadrature  in  appropriate  cases  when  the  solution  is  expanded  in  a  perturbation  series.  These 
perturbation  solutions  are  the  subject  of  Chapters  3  and  4. 

With  0  known,  the  pressure  can  be  computed  from  Equation  (2-14) 

p=yjtf-  U2  -(2q0+  V<t>)  'V<t> 
and  then  the  force  on  the  propeller  can  be  found  from 
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F  =  2  ff  p(-n)ds 

b=0  JJ  SB 

=  7  S  [2qQ  +  W  ’  V<}>/nds 

2  b=0  JJ  SB 


(2-56) 
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where  0*  and  0  are  the  values  of  0  on  the  suction  and  pressure  side  of  the  blade  at  the  points  of  inter¬ 
section  of  the  line  £j,  co  and  the  blade.  By  analogy  with  Equation  (2-25),  the  circulation  at  the  points  of 
intersection  of  the  line  co  and  the  blade  is  defined  as 

ryS;,  s?  =#£;.  E+,  ZS)  -  (Mj,  ET,  co)  =<f  -  <p~ 


rb'Z 

=  I  (V<p+  °d!+  -  V0-  °d!_~) 


(2-59) 


where  the  integration  starts  at  a  point  on  the  edge  of  the  lifting  surface  or  wake  and  goes  to  the  point 
c5  Similar  to  Equation  (2-37),  we  define 

0^;,  =<p+  +  0~ 


(V0+  °d!+  +  V0- 


(2-60) 


which  gives 


I  dCr  n  '  V^B  +  B  J  ^  '  ,  +  srn  ifip  ■ 


+  (2^^  •  V0£  +  SJ(j)+  •  V0+  +  V0  •  V0  )  sin  \p 


(2-61) 


The  thrust  coefficient  is 


7/2  pU2  irR2 


Z  /•  1  *  i 

I  d—  2  cos  p  y  U2  +  SI2  go2  cos  (<p  -|3)  0  T(ZS) 
nU2  JO  R  L  P 

-  rc‘IR  |  „  /8r« 

jc2/r  [  p  at7 


/ -  r  **  l0VB  dE 

+  2\U2  +  Q?  gj2  sm  (p  Icos  (ip  -&)[— - ^ 

jc2/r  |  p  \^7  at7 
30  ^  3^\  /  dEy~> 

+ ir  J  *  ^'®-2  ■  + *’»  ,r® 

*4  ?/  s / 

+  Vr5— )>d  — +cox^)  V05  »vr5(7—  +X/«(P 


* 


(2-62) 


fcys/  ( 
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Similarly  (he  moment  acting  on  the  propeller  blades  is 


M 


Z-l 

2 

b=() 


ff.  P(*X(-2))dS 
jj 


P 


i 


(2q()  +  70/  •  70 


(s_  X  n)  dS 


And  the  moment  about  the  /  axis  is 


Q  =  M 


(2q()  +  70/  •  70 


i  '(sXn)  dS 


(2-63) 


(2-64) 


This  equation  is  not  of  fundamental  use  in  the  following  development  and  is  not  further  reduced. 

Before  the  previous  formulation  can  be  considered  the  complete  solution,  two  points  assumed  in 
the  development  must  be  checked:  first,  that  the  integration  over  the  region  at  infinity,  which  was  neglected 
in  going  from  Equation  (2-18)  to  Equation  (2-19),  gives  at  most  a  constant;  second,  either  that  the  solution 
is  unique  or  that  further  conditions  to  establish  uniqueness  need  be  specified. 

To  examine  the  integration  over  the  boundary  at  infinity,  an  area  consisting  of  a  finite  cylinder  of 
radius  W  with  ends  at  x  =  ±  X  is  considered.  The  integration  over  this  area  as  X,  W  ■+  +  00  is  to  be 
evaluated.  By  uniqueness  proofs^  for  nonlifting  problems,  the  integrals  of  0  and  70  are  zero  for  that 
portion  of  the  solution  arising  from  the  integration  over  the  body.  Any  constant  arising  from  integration 
over  the  wake  also  gives  zero  on  the  cylinder  and  the  upstream  surface.  However,  the  downstream  surface 
must  be  investigated  in  detail  since  neither  0  nor  70  go  to  zero  as  jc  -+  +  °°. 

At  the  downstream  surface  the  vortex  sheet  rolls  up  as  it  does  for  wings;  see  Cummings; ^  How¬ 
ever,  the  total  circulation  remains  constant  and  the  lines  of  constant  circulation  become  asymptotic  to 
regular  helical  vortices  of  constant  radius  and  pitch.  The  integral  over  the  downstream  surface  requires 
evaluation  of  the  two  quantities 


and 


^Cummings,  D.,  “Vortex  Interactions  in  a  Propeller  Wake,”  Massachusetts  Institute  of  Technology, 
Naval  Architecture  Department  Report  68-12  (Jun  1968). 
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(2-66) 


where 


J 2  =  lim  | 

X+°o 


r=Xi+yl+zk 
s=yfj  +  z'km 


and  the  bar  through  the  integral  sign  means  the  vortex  itself  is  excluded;  0  and  (j>x  are  undefined  on  the 
vortex  surface,  except  in  the  limit. 

In  Appendix  A,  the  form  of  the  potential  for  a  vortex  distribution  on  a  regular  helix  is  examined, 
and  it  is  shown  that  both  0  and  0X  have  the  same  form  far  downstream.  Since  the  velocities  are  known  to 
exist  in  the  sense  of  a  Cauchy  principal-value,1  we  assume  the  potential  does  also.  In  particular, 

this  means  both  0  and  (px  are  bounded  in  the  fluid. 

The  region  of  integration  can  be  divided  into  two  regions;  one  bounded  by  a  circle  of  radius  A  and 
the  other  one  the  area  outside  this  circle  and  going  to  infinity.  By  taking  A  sufficiently  large,  integration 
over  the  area  exterior  to  the  circle  can  be  made  to  yield  as  small  a  value  as  desired,  since  0  and  its 
derivatives  go  to  zero.  In  the  interior  of  the  circle,  the  boundedness  of  0  and  the  power  of  X  in  the 
deonimnator  insure  that  the  limit  as  X  00  gives  a  zero  value  for  the  integral.  Thus  the  integral  over  the 
surface  bounding  the  fluid  does  not  contribute  to  the  expression  for  0  given  in  Equation  (2-19)  and  the 
statement  to  that  effect  just  previous  to  Equation  (2-19)  is  justified. 

To  examine  uniqueness  of  the  potential,  two  different  solutions  are  assumed  to  satisfy  Laplace’s 
equation  and  the  boundary  conditions.  From  these  a  third  solution  is  constructed  by  subtracting  them. 

This  third  solution  also  satisfies  Laplace’s  equation  and  has  zero  normal  velocity  on  the  fixed  body. 

However,  unless  the  two  solutions  have  the  same  circulation,  the  difference  solution  will  have  a  shed  vortex 
sheet  with  a  finite  value  of  normal  velocity  on  it.  Hence,  the  two  solutions  must  have  the  same  value  of 
circulation  in  addition  to  satisfying  the  boundary  conditions.  Identical  positions  for  the  shed  vortex  sheet 
are  assumed.  Then  the  difference  solution  has  no  shed  vortex  sheet,  and  0  0  for  \r\  -+  °°. 

Green’s  first  identity  in  the  form 

I  (<pv2(t>)  +  V0*V#/T  =  JTf  0  n  *70  dS  (2-67) 

V  V 


then  shows  that 


J]J  V0°V0c/r  -  0 
V 


(2-68) 


Since  V0*V0  is  a  positive-definite  form,  it  follows  that  V0=  0.  and  hence,  that  0  (the  difference  of  any  two 
solutions  with  identical  circulation)  can  be  at  most  a  constant  throughout  the  fluid.  Therefore,  any  two 
solutions  of  the  problem  which  have  identical  shed  vortex  sheets  differ  at  most  by  an  unessential  constant; 
thus  the  solution  is  unique. 
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CHAPTER  3 

REGULAR  PERTURBATIONS  IN  PROPELLER  THEORY 


A— GENERAL  FIRST-ORDER  SOLUTION 

In  this  chapter,  the  regular-perturbation  problem  is  examined  and  the  application  of  the  formulation 
to  both  design  and  performance  calculations  is  discussed. 

For  the  regular-perturbation  problem,  the  propeller  is  assumed  to  deviate  little  from  a  helical 
reference  surface  of  variable  radial  pitch  angle  \p^[o5).  (This  is  the  same  description  as  in  Chapter  2.)  In 
design  this  surface  might  be  the  approximate  position  of  the  shed  vortex  system,  and  in  performance  calcu¬ 
lations  it  would  be  the  geometrical  pitch  surface.  On  the  cylindrical  surface  co  =  constant,  a  left-handed 
coordinate  system  y ,  to)  is  constructed  (as  was  done  in  Chapters  1  and  2)  with  £y  along  the  inter¬ 
section  of  the  cylinder  with  the  reference  surface  so  that  £y  increases  with  increasing  jc.  The  variable  ^  *s 
on  the  surface,  normal  to  the  £y  axis,  so  that  increases  with  increasing  d.  The  cylinder  is  rolled  out  in 
Figure  5.  Unit  vectors  (e y,  e^e  are  in  the  (%j,  5)  direction,  respectively. 


Figure  5  —  Coordinate  System  for 
Describing  Blade  Shape 


In  terms  of  these  new  independent  variables,  the  position  vector  is  (Chapter  1) 

r_  =  ((R  +  £y  sin  cos  Vp)±  ~  w  sin  (d  ^  +  W 

Sjcosvp  +  bsinvp.  _  / 

- 7 1  +  00 cost  0^  +  W 


CO 


1/  cos  I Pp  +  %2  sin  'P, 


CO 


(3-1) 
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The  blade  is  assumed  to  be  described  as 


£2  =  e  E(\j,  <3) 

=  €  ±  6  eT  (%1  •  ^ 


(3-2) 

(3-3) 


where  E 


€ 


is  the  function  describing  the  blade  thickness  and  camber  for  the  appropriate  restrictions  on  the 
independent  variables  and  to, 

is  a  small  parameter  describing  the  departure  from  the  surface  ^  -  0, 
is  the  camberline  function,  and 
is  the  thickness  function. 


The  description  of  the  blade  offsets  differs  from  the  presentations  in  Chapters  1  and  2  by  explicitly  making 
them  0(e), 

The  position  vector  of  a  point  on  the  blade  surface  is  given  by 


s(£jy  co)  =  (<R+  sin  -  6  E  cos  <Pp)i  -  <5  sin  \  6  ^  +  W 
cos  <pp  +  eE  sin  ipp\ 


j/  +  co  cos  I  6  +  W 


CO 


(• 


sin  ^9 


(3-4) 


+ 


Si 


cos  i fp  +  e  E  sin  i 


CO 


Two  tangent  vectors  to  the  blade  surface  are  given  by 


ds 


9 s 


Hence  a  vector  normal  to  the  blade  surface  can  be  constructed 


(3-5) 
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cosip  n  +  e  !■:  sin  <pr  M 


oj  —  1+  —  f  y?.,  +  € —  sin  y? 
cl  Co/  dco  \ 


+  <r  E- 


0 /•' 


Confd 


Note  that  even  for  \p  =  d?  =  IV  -  0  and  I:  -  E(%j)  a  term  appears  in  the direction.  This  is  be¬ 
cause  the  cross  sections  arc  measured  on  the  surface  of  cylinders  of  different  radius.  For  oc^^thc  radial 
component  of  the  normal  vanishes  for  the  stated  conditions. 

The  boundary  condition  on  the  body,  Equation  (2-55),  gives  two  equations 

1/2  (N+'Vfr  -  jV-  )  =  -j/" J2  +  or  n 2  /  sin  (ipp~P) 


e - cos 

d£]  P 


n, - - —  dEr 

1/2  (N+  °V(jx-  +  N-  ’V<k-)  =  \f  U2  +  OJ2  SI2  e - cos  fy>  -ft/ 

r  P 


The  further  assumption  is  made  that 


(Vp-fr)  =  0(e) 


(This  assumption  is  necessary  only  in  design,  since  it  insures  that  the  camber  and  the  potential  are  0(e).) 

As  yet  no  approximations  beyond  the  inviscid,  incompressible  fluid  and  the  steady  velocities  have 
been  made.  However,  to  obtain  the  regular  perturbation  solution  only  the  terms  to  order  e  are  retained. 
Hence,  a  point  on  the  surface  is  given  by 


_s  =  jp  +  e_sy  +  0(ez) 

_  /  cox 

=  (6\  +  sin  <pp)JL~  cosin  I  0^  ^  H7 - J  j 

~  L  . I  . 


+  co  cos  (  6^  +  W+ - j  k  +  e  j  -E  cos  <Pp±_ 


(  i]  COS  \  E  sin  y> 

CO  cos  \dh  +  W  +  - j -  j 


/  cos  <Pj\  E  sin  ip  j 

-  oj  sin  dh  +  W  + - -  I - -  k  >  +  0(e2) 
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■^Actually  it  is  the  difference  in  velocity  across  the  blade,  Equation  (3-13),  which  should  be  taken  0  (e), 
which  is  necessarily  true  only  in  design.  If  -0  =£  0  (e),  as  would  be  true  in  performance  calculations, 

then  an  integral  equation  must  first  be  solved  to  find  the  pressure  distribution  corresponding  to  an 
uncambered,  thin  lifting  surface  at  an  angle  of  attack.  The  remaining  problem  corresponding  to  the  camber 
and  thickness  is  then  0(e). 
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whore  .S\>  is  the  projection  of  the  blade  outline  on  the  surface  =  <1 


N+  '(?'*())  =  (&  X)  cos  i p ^  sin  |  y  cos  0  +  z  sin  0 


/  dP  mr  Vp  _ dW  d(d 

+  [  £/  — - co  —  sin  i p.+  cos  <p  - 

\  '  d~  2  ~  d~~  p  1  d~ 

\  CO  7TC O  CO  CO 


I  d{f>p  ~dw  . 

z  cos  d  -  v  sin  6}  -  colt, -  -  co  —  sin  \p  +  cos  <pt.  — 

V 1  <j~~  P  P  d — 

/  \  CO  UC 0  uco 

=  l(x  -  <R  -  sin  vp)2  +  (y  +  CO  sin  d)2 


(z  -  cocos  Qj~  I  2 I2 


e  =  d  h  +  W  +  ■ 


%}  cos 


As  in  Chapter  2,  the  jump  in  potential  across  the  blade  surface  can  be  interpreted  as  a  local  circu¬ 
lation.  In  terms  of  the  velocity  jump  across  the  surface  it  is  given  by  Equation  (2-59) 


I  B  (%}'  ~  }  0  ^  “  0 


(V0+  «  dH+  -  V0  -  ♦<//  -} 


■  <Zj.  co; 


(q+  -  q -)  ‘d!  0  +  0  (e2) 


(3-13) 


where  the  path  starts  from  any  point  along  the  leading  edge  and 

rT^0> 

r (z0)  =  e  J  (q+  -  q-)  •  d[  Q  +  0  (e2) 
LE 

Similarly,  the  integral  in  Equations  (2-21)  and  (2-50)  reduces  to 


jljt 

4tt  JJ- 


io  -£L 


d%,  d  co  = 


2it  JJSd  + 
B0 


XI2  co“  |  ^  ^ 1  cos  fa  _g)  ? 

1  + -  |  - cos  ({pP  p/  d%t  dco  +0  (e2)  (3-14) 
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The  position  vector  of  the  shed  vortex  sheet  is  determined  by  Equations  (2-46)  and  (247),  which 
to  first-order  in  6  become 

(0)  (1)  9 

f^z0>  ® )  =  f 2  (z0’  ® )  +  € ?2  (z0’  d)  +  0  (*2) 


-  xte/z0 )  +  O  ^  "  ®TE  <z0» 


fe  i*Ao>  «  } 

+  e  1/2  l  { - —  ifre(S0)}da+0(e2) 

J  eTE[  n  n2Zn2  ~  I 


(3-15) 


f](z0,  d)  =  fj  (zq,  6)  +  efj  (zq,  6)  +  0  (e2) 


(3-16) 


0 

=  z0-+l/2Sl  j  e  da  +  0  (e‘ 

•'A 


i=io+eh+°(e  > 

u 

=  z0- oS+  (xTe/z0>  +77^  -  ®TE  (zoHL 


+  e  (1/2  SI  I  da)e~+  (1/2 


/J 


*jy 


(3-17) 


hiy\da>L  * 0,(2 1 


«V 


However,  since  T(zq)  in  Equation  (240)  is  of  order  e  the  first-order  solution  for  the  potential  of  the  shed 
vortex  sheet  requires  only  The  zero  order  normal  to  the  shed  vortex  sheet  becomes: 

+  _  U  3  / 

jV0=— x- —  =  (+z0£Q+—i)X  (e^s+  —  \  Xje 

30  dzn  u  °z0\ 


3  /  U°TE 

*-z0L+z0b  \XTE  ~T~ 


U 

+  —  e  n 
O  -V 


(3-18) 


45 


and  t lie  dot  product  in  the  infinite  integral  is 


+  u  V  _ 

J^.()  (£'$())  ~  zo  /  v7’/-;  +  q  (a  O'} '/•)  xl  +  M  f  sin  0 


,UTF  \  t/  [  i) 

- )  +  cos  a — }+  ms  0  J  zn  cos  a  — 

ft  /  ft  0  (I: 


>-•»  \n: 


sin  a  U/£lj  z^j~ —  -  I  xff- 


and  the  denominator  is 


v  (xtf  ^  ^  ^77.  /)}  “  *  +  z (r 


2  cozn  cos  (0  -  a+0rr)\ 


(3-20) 


where  the  position  vector  of  an  arbitrary  point  in  space  is  expressed  in  cylindrical  coordinates.  This  is 


done  to  facilitate  the  integration  since  past  experience 


11,13,14  : 


indicates  that  in  some  cases  the  integration 


can  be  done  analytically  for  such  a  representation. 

Putting  together  the  previously  described  pieces,  one  finds  the  total  first-order  potential  for  the 
disturbance  due  to  the  blades,  omitting  the  e 


7- v  it 

4n  b=0  JJ  S 


d^jtl  Cx) 


(q+  -  q-j  •  di  (j  N+  .  fl_y  _  2y/ y2+n2  %2  -~cos  fay  (i)\r-s0\- 


(3-21) 


/  2-t  rR  r 

—  2  J  Vfzo>dzo\ 

47T  h=n  J  n  J  < 


b=0  J  0 


■°°  r]lo  })k  \  (ry 

- X  -  I  •  — - ;  da 

°Y/A 9a  '9z0/  jr-j^’ 


This  equation  is  the  first-order  solution  for  the  regular-perturbation  problem.  In  the  next  two  sections,  the 
use  of  this  equation  in  design  and  performance  predictions  is  discussed,  and  in  the  last  section  the  second- 
order  solution  is  outlined. 

The  force  on  the  propeller  blade,  Equation  (2-56),  is  to  first-order  in  e 


(q0  •  vrB)Hod$l  d€y 


(3-22) 
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and  the  thrust  produced  by  all  the  blades  is 


Z-7  rr\  * - 

T  =  -  2  1  •  £  =  epZ  cos  y  yU2+Q?  go2  cqs  (\p  - P)  V(gS)  cl  go 
h=0  Jn  y 


+R 


(3-23) 


The  torque  can  be  computed  in  a  similar  manner  but  it  is  not  used  in  the  following  sections. 

A  second-order  theory  follows  by  straightforward  application  of  the  previously  described  procedure 
when  thickness  and  camber  are  expanded  in  perturbation  series.  The  second-order  effects  include  con¬ 
sideration  of  the  first-order  position  of  the  shed  vortex  position,  Equation  (3-17).  Second-order  design 
theory  is  discussed  in  Section  D. 

B-DESIGN  PROBLEM 

For  the  design  problem  considered  here  the  total  thrust  produced  by  the  Z-blades  is  specified  at 
one  operating  point,  as  well  as  the  following  geometrical  data:  ■  , 

1.  Blade  outline 

2.  Chordwise  thickness  function  and  radial  thickness  ratio 

3.  Number  of  blades 

4.  Position  of  blade-reference  line 

The  designer  also  specifies,  to  within  a  scale  factor,  the  chordwise  difference  in  pressure,  i.e.,  the 
pressure  on  the  upper  surface  minus  the  pressure  on  the  lower  surface,  which  reduces  to 


-  +  P 

p  -p  =j 


(q+P  -  (q-P  J 


Since  q  =  qg+  V<p ,  the  pressure  difference  is 

j  q0  •  V(<h  -  <P+)  -  i/2f(q<p+p  -  (q<t>p] 


p+  =  -p 


=  +p 


3.0  •  irB*i/2vitB  •  wB  J 

Also,  since  qg  has  a  zero-order  component  only  along£y,  the  significant  term  of  VlTg  is  the^ 
component.  The  velocity  difference  vTg  can  be  resolved  into  the  three  orthogonal  components  (y,  o )f 

VI" g  =  yej  +  pNq+  ofNgXej) 

From  the  boundary  condition,  Equation  (3-7),  the  component  of  p  is  known  to  first  order 

- - dET 

p  =  2  eyU2  +  oj2  Q?  — —  cos  1'P„-&)1(Nq'Nq)  +  0  (e2) 

of/. 


However,  to  0(e) ,  no  use  is  made  of  this  knowledge. 


i 
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since  V \'B  -  /YcA  7  *  (YcA  s///  =  ky„|  mv  /Vy  /J^y  *  j^ 


avn  thc  pressure  distribution  is 


p  p 


U~  +  SI2  ZS2  (sy(V  cos  +  (He2)) 


Hence  to  first  order,  the  chordwise  component  of  the  velocity  difference  is  related  to  the  pressure  difference 
across  the  blade  by  the  equation 


P  -  P 


p  Vi! 2  +  SI2  oo2  cos  ($p-($) 


(3-24) 


For  the  first-order  problem,  no  information  need  be  specified  about  the  other  components  of  the  velocity 

,  * 
difference. 


*The  component  0  can  be  found  from  V  since  is  known  in  terms  of  7^  stated  previously.  An 

alternative  is  to  note  that  when  a  tangential  discontinuity  in  velocity  exists  across  a  surface,  one  can  define 
the  surface  as  a  vortex  sheet  with  a  vorticity  distribution  given  by 

Eo 

A ’  nxfq* -q~)‘ - X  hi,  *  0  fKXe,)l 

~  ~  Vg I 


■  °H  *— Wo* Si> 


Since  j\must  satisfy 


V  •  A  =  0 


the  component  o  can  be  found 'from  the  partial  differential  equation  governing  the  components,  The 
equivalence  of  both  these  approaches  is  demonstrated  in  Chapter  5. 
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The  local  circulation,  which  is  the  Quantity  needed  in  the  calculations,  is  given  by  Equation  (2-59) 
or  to  first-order  by  Equation  (3-13).  Since  the  chordwise  component  of  velocity  difference. is  known,  d £ 
the  tangent  vector  in  Equation  (2-59),  can  be  taken  as 


o  9-  /  &  \ 

4=ii  d%j  =  —  dij  =^7  + e 


Hence  Equation  (2-59)  becomes 

4i[ 


/  9  et  \ 

+  V(4>+  +  <P-)  •  (e  — S.2J  I  d$l 

f  {  Tiy  +  MjV ^ofNgXej)  °( e; 


c2  m 


+  V\pB  °  e_2 


+  e  rr —  e? 


*,  r 


rh  (  3 ec  a et\ 

=  J  y  +  €~^-(L1+  °Z2  °@0X2l))+  V^B  '  £2e^-Jd^l 


When  the  quantities  are  expanded  in  a  perturbation  series,  one  finds 


rb  /  a  a;,  / 

-  e  J  7^  c?£7  +  J  ^  7^  +  +  dV £2 


'(EoXejA  +  V  -e2  —J  d$j  + 


dp(Cos2^)  _ 

-£2m@0X£l)ail-z}‘  Z 

tfco  2  71  CO 


cj- — sin  ip v  +  cos  <p>  - 

rfco  y  d  00 


and  £2  *  ~  v  *s  ^rom  Equation  (3-26).  For  the  first-order  problem,  only  the  leading  term  stated 

previously  is  needed 

r /'  *  e  | ;  jt<>  d(, 
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When  Equation  (3-24)  is  integrated  over  the  chord,  the  total  bound  circulation  distribution,  still 
with  the  unknown  scale  factor,  is  obtained.  This  unknown  factor  can  be  determined  from  the  equation  for 
thrust,  liquation  (3-23).  Hence  everything  needed  to  calculate  the  first-order  perturbation  potential  is  known. 

What  is  not  yet  known,  however,  is  the  magnitude  of  the  camber  and  the  blade  pitch  which  will 
produce  the  specified  thrust.  To  obtain  these  values,  the  boundary  condition,  Equation  (3-6)  is  used.  To 
first  order,  this  equation  is 


dA 


Jfo 

=  sin  (Vp-P)  -  1/2 - — — — — — 

Vu2  ~+  zs2  n2 


(3-25) 


where  \jj  is  a  principal-value  integral  given  by  Equation  (2-18). 

For  conceptual  consistency  the  blade  trailing  edge  must  lie  on  the  reference  surface  so  that  the 
trailing  vortices  spring  from  the  trailing  edge.  Thus 


This  camberline  will  consist  of  a  shape  measured  from  the  nose -tail  line  and,  in  general,  an  angle  of 
attack.  This  angle  is  added  to  the  pitch  angle  of  the  reference  surface  to  give  the  total  blade  pitch.  This 
incremental  angle  is  called  the  ideal  angle  of  attack. 

The  knowledge  of  the  camberline  completes  the  first  approximation  for  the  design  problem.  The 
regular  perturbation  solution  involves  an  integration  to  find  the  slope  of  the  camberline  and  then  a 
straightforward  integration  of  this  slope  to  find  the  shape.  However,  a  new  term  appears  in  Equation  (3-25) 
which  should  be  explicitly  mentioned  before  proceeding  with  the  analysis. 

From  the  equation  for  0  in  Equation  (2-18),  one  sees  that  Equation  (3-21)  gives  1/2  0  when  the 
principal  value  for  points  on  the  blade  reference  surface  is  taken.  Hence,  if  the  gradient  of  0  in 
Equation  (3-21)  gives  the  velocity 


V0  -  ue j  +  V£2  + 


(3-26) 


and  if  for  r  +  Yq 

1  _  _  _ 

—  7  i//g  =  0,  oi)e}  +  l >($j,  0,  u>)  e2  +  wf%},  0,  oi)e~ 

then  Equation  (3-25)  becomes 
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J6?  ~dW  ( dP  cn%2  'Pp\  _ 

+  —cos  vp  +  u—sin  ypp  -  Zj  |  — - —  I  w(£j,  0,  w > 

dZ3  do5  \dcj  2ttco  / 

In  the  literature,  intuitive  arguments  are  used  to  derive  only  the  first  term  on  the  right-hand  side  of 
Equation  (3-27).  The  formal  procedure  presented  here  indicates  that  another  term  also  contributes  to  the 
camberline  shape.  For  the  term  containing  a  multiple  of  the  distance  from  the  reference  line,  the  effect 
will  be  most  pronounced  for  points  near  the  blade  edge. 

In  the  discussion  section  of  Chapter  5,  an  estimate  of  the  contribution  from  this  term  is  made.  In 
addition,  other  aspects  of  the  design  problem  which  simplify  the  integrals  are  discussed. 

C— PERFORMANCE  CALCULATIONS 

To  determine  the  performance  of  a  propeller  over  a  given  /-range,  it  is  assumed  that  the  position 
vector  of  a  point  on  the  surface  is  given,  that  is,  the  geometry  is  completely  specified.  The  unknowns  of 

interest  are  the  thrust  and  torque  as  a  function  of  J  and  possibly  the  pressure  distribution.  The  pressure 

e  3  Et 

distribution  near  the  leading  edge  will  not  generally  be  accurate  because  -  |  is  in  general  infinite  and 

a£;  LE 

hence  not  of  order  e  as  assumed  in  the  derivation.  In  the  design  problem,  this  point  was  not  critical  but  in 
performance  calculations  it  will  generally  give  results  for  the  pressure  distribution  which  are  not  uniformly 
valid  at  the  leading  edge.  However,  it  is  integrable  and  hence  the  total  thrust  can  be  found. 

In  calculating  the  propeller  performance,  the  circulation  is  unknown,  and  the  zero-order  position 
of  the  shed  vortex  sheet  is  known.  Modes  for  the  chordwise  distribution  of  velocity  given  in  Equation  (3-24) 
can  be  assumed  which  satisfy  the  Kutta  condition  and  which  have  variable  spanwise  coefficients.  Then,  in 
principle,  the  integral  equation  formulation  in  Equation  (3-25)  can  be  solved  for  the  spanwise  coefficients. 
Sugai^a  was  able  to  successfully  evaluate  numerically  such  a  procedure.  Murray^  used  an  iteration  pro¬ 
cedure,  starting  with  blade-element  theory  (Glauert^)  with  a  reduction  factor  found  necessary  for  con¬ 
vergence.  As  already  mentioned,  Murray’s  performance  calculations  did  not  converge  to  the  values  set  in  his 
design. 

Presumably  the  exact  solution  in  Equation  (2-50)  could  be  used  as  the  integral  equation  without 
much  increase  in  complexity  but  no  one  has  done  this  yet. 

^Glauert,  H.,  “Airplane  Propellers,”  Division  L,  Aerodynamic  Theory,  Edited  by  W.F.  Durand,  J. 

Springer,  pp.  169-360  (1935);  also  published  by  Dover  Publications,  Inc.,  New  York  (1963). 
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Rather  than  solve  an  integral  equation  it  would  be  preferable  to  obtain  a  solution  which  could  be 
constructed  as  quadratures  as  was  done  in  the  design  problem.  Such  an  expression  is  found  in  the  next 
chapter. 

D-SECOND-ORDER  DESIGN  THEORY 

The  second-order  theory  follows  by  retaining  terms  to  order  e  in  the  previous  development.  Most 
of  the  second-order  effects  are  known  from  the  first-order  solution;  however,  the  second-order  increment  in 
circulation  must  be  evaluated.  This  is  found  by  considering  the  second-order  thrust 


T  =  Tjd +  e2  V(2))  +  e2  T2 

where  T is  given  by  Equation  (3-23)  and 


(3-28) 


fR 

I  duo  I  i 

J  0  h2  _ 


V \!/  •  vr,n!  )  cos  +  2  Yu2  +  Si2  Z32 


sin  (ipp-P)  cos  <pp  e2  •  +  cos  (<pp~P)  sin  <pp 


dE/v  a  a  et 

at;  3  Hj  3 


(3-29) 


All  quantities  in  this  integral  are  known  from  the  first-order  solution.  Since  in  design  the  thrust  is  set  and 
since  T is  the  set  value,  the  second-order  terms  must  vanish: 


e2  Tj(T(2>)  +  e2  T2  =  0 


which  fixes  the  value  of  r727 


,c;(-  He111  ■ 

=  -  I  i/2— . -  — —  — 

J  c2  L  rt^+n2w: 


2  COS  (ipp-i 3) 


(3-30) 


+  tan  ip 


f 

ar^777  a  eJ1)  a  e. 

3  \i  3 tj  d Sj 


+  tan  (Vp-P)  i2  *  V  T  J1)  d 
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In  design  the  form  of  the  chordwise  circulation  is  assumed  to  remain  constant,  i.e.,  7 (%jt  to)  in 
Equation  (3-24)  is  a  fixed  chordwise  variation;  however,  the  amplitude  of  the  function  may  change.  The 
scale  is  fixed  by  the  magnitude  of  the  bound  circulation.  If  7. is  normalized  to  give  unity  when  integrated 


over  the  chord 


C1 

J  c  ^ 


dtj  =  / 


then  the  scale  factor  will  be  just  the  total  bound  circulation.  In  this  case 


(q+  -  q~)  °  ej  =(T<1)  (co)  +  eArW  (gj)J  y*($j,  T3) 


r 1 

rJV-iWfZS)'  \ 

c/co) 

where  7 *  ($j,  to)  is  the  normalized  value  of  Equation  (3-24).  Now  we  have  previously  found  that 


h  r  c)Er 

f  y<2)  + - - 


(fJ1)  +  (/*)  e2  •  (NqX  ej)) 


+  e2  •  -  dfy 


T  r1\dEc 

Y(G)  =  e  +  e2  AY^2)  +  < — (fJ^  +  e2  •  (N q\ e/) 


+  e2  •  - 1  d$j 


(3-32) 


r  1 

A r (2)  =  r (2)  ..  f  J - (jjL(l)  +  a(l)  ^  (u+x  ^)) 

Jc,  d%l 


+  e_2  *  V  ^jJ1) - [  dtj 

^ 1 


(3-33) 
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All  quantities  appearing  under  the  integral  sign  in  liquations  (2-40)  and  (2-50)  are  expanded  in 
perturbation  series  and  the  second-order  contribution  found.  The  necessary  expansions  have  been  given 
previously  and  in  Part  A.  In  these  calculations  it  is  assumed  that  the  thickness  function  is  given  and  constant 
at  order  e.  The  quantity  sought  is  the  second-order  correction  to  the  camber.  This  is  obtained  from  the 
boundary  condition  on  the  blade  set  in  Equations  (2-10)  and  (3-6).  Because  of  the  length  of  the  ex¬ 
pression,  it  is  not  given  explicitly  but  it  follows  directly  from  the  previously  described  development. 
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CHAPTER  4 

SINGULAR  PERTURBATIONS  IN  PROPELLER  THEORY 


As  described  in  the  introduction,  the  singular  perturbation  problem  examines  the  solution  as  the 
chordlcngth  goes  to  zero.  At  a  fixed  pointer  ~  in  the  How  field  relative  to  the  propeller,  the 

chordlength-to-diameter  ratio  e"  goes  to  zero.  Since  the  geometry  is  unchanged  except  for  chordlength, 
the  blade  forces  go  to  zero  in  this  limit.  Hence  the  How  disturbance  vanishes  as  the  lifting  surface 
vanishes.  The  solution  in  this  case  is  just  a  lifting  line  (as  the  chordlcngth  goes  to  zero,  thickness  effects 
give  a  second-order  dipole  distribution  rather  than  a  first-order  source  distribution  and  hence  will  appear  in 
the  higher-order  terms  of  the  solution).  This  is  the  outer  flow.  The  flow  near  the  body  is  the  inner  flow 
and  can  be  found  by  stretching  the  coordinates  near  the  body  so  that  the  flow  field  is  magnified  sufficiently 
to  recover  the  necessary  details.  In  the  inner  limit  the  stretched  variables  are  held  constant  as  the  chordlength 
goes  to  zero. 

A— FIRST-ORDER  TERMS 


Outer  Flow 


The  starting  point  is  the  outer  flow  which  is  found  by  straightforward  linearization  of  the  exact 
solution  given  by  Equation  (2-50).  Both  0  and  V  will  be  of  the  same  order  in  the  integrals  of 
Equations  (2-49)  and  (2-50).  However  the  integral  over  the  blade  area  is  of  higher-order  than  the  integral 
over  the  shed  vortex  sheet  because  the  area  of  integration  goes  to  zero  with  the  chordlength.  The  zero- 
order  integration  over  the  shed  vortex  sheet  reduces  to  the  lifting-line  results,  i.e.,  aj^eT(iQ))  +  d ^  = 
2nb 


Z 


with  the  zero-order  position  of  the  vortex  sheet. 


Formally,  the  circulation  is  expressed  in  a 


perturbation  series  of  the  form 


T(z0,  e)  =  8  ^e)  (zq)  +  82(e>  Y^  (zQ)  +  •  •  •  (4-1) 


where  e  is  the  maximum  chordlength-to-diameter  ratio  and  6  n(e)  is  an  ordered  set  of  linearly  independent 


gage  functions,  with  the  property  6  =  o(8  )  as  e  -+  ft  A  general  form  for  them  is  e 

a*. 

JW„3  -  h 


san 


(In) 


¥ 


where  anm  are  constants.  Normally  the  terms  with  af  *  =£  0  are  considered  trans- 


cendentally  small  and  are  ignored.  The  coefficients  af}  r  will  be  found  as  the  solution  evolves.  The 
position  vector  of  the  curves  of  constant  circulation  in  the  shed  vortex  sheet  is  also  expanded  in  a 
perturbation  series  "  ■ 


(4-2) 


p:(j.  0:  ej  =l0(:0.  0)  +  6,(e)lj(z0,  0)  +  S2(e)  $2(zQ,  ())+•■■• 


From  Equation  (2-46)  and  (2-47)  the  zero-order  solution  for  a  blade-reference  line  lying  along  a  radial  line 

* 

in  the  a*  =  0  plane  is 


lo  =  zo£€J(e>+ 


and  consequently  the  normal  is  given  by 


at  H  (  u  \ 


zoi+7>ie  +  °<bi> 


Hence  for  the  outer  flow,  the  lowest-order  solution  in  Equation  (3-50)  is 


<p(x,  y,  z;  e)  =  5;  <p}fr)  ±  •  •  • 

R 

4>i(r)  =“  J  (z0)  G(r,  zQ)  dzL 


where  G(r,  Zq)  =  2  I  — 

b=oJeb 


„  r  i  u  ~ 

J  x  ~—(a-Qb)  zQ-~u>sin  (d-d) 

6b  {  [*  ~  2  +  Zs2  +zo2 ' 

Z~]  f°°  |  3  U  3  \  1 

b=o  o  \°te  +  nz0befrdy 


.7  ?  5/2 

+zq  -  2  Z ozq  cos  (6- a)  j 


r  =  \/(x  -  —7 i2  +  Zo2  +  zq 2  -  2zq  Zo  cos  (d-y-6b) 


The  reason  for  this  restriction  on  the  reference  line  will  be  explained  when  the  inner  flow  is  examined. 
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This  expression  is  also  derived  in  Appendix  A  from  the  Biot-Savart  Law  applied  to  an  elemental  horseshoe 
vortex.  Isay^  presents  a  similar  equation  which  he  bases  on  the  work  of  Yamazaki.^ 

Equation  (2-50)  also  contains  another  term  which  appears  to  be  of  order  e;  the  integration  of  the 
source  terms  over  the  boundary.  The  first-order  effect  in  these  integrals  reduces  to  an  integration  over  the 
radius  and  an  integration  of  the  thickness  slope  over  the  chord.  Since  only  closed  profiles  are  considered, 
this  integration  is  zero,  and  the  contribution  is  second  order.  In  the  outer  flow  one  expects  dipoles  rather 
than  sources  which  is  confirmed  in  the  previous  statement. 

Inner  Flow 

Unfortunately,  the  potential  i$  not  uniquely  determined  since  the  circulation  is  unknown.  The 
circulation  is  determined  by  details  at  the  body  and  hence  must  be  determined  from  an  inner  flow  for  which 
the  body  boundary  condition  is  not  lost  as  the  limit  of  zero  chordlength  is  taken.  Because  the  spanwise 
scale  is  not  changed  as  f  *»  0>  a  stretching  is  needed  for  only  the  chordwise  variables.  A  suitable 
stretching  for  the  reference  blade  along  the  z  axis  is 

x  *  x/e  and  y  -  y/e  0  <z  <R  (4-7) 

where  x  and  J  are  held  constant  as  e 0,  The  z  variable  remains  unchanged. 

The  previously  described  transformation  applied  to  the  Laplace  equation  gives 

90  90  ?  92  0 

e2 — -  (4-8) 

dx2  972  9  z2 

Now  it  can  be  seen  why  rake  and  warp  have  been  excluded  from  the  analysis:  If  rake  were  included,  a 
displacement  varying  with  z  would  have  to  be  made  in  x  before  the  stretching  in  Equation  (4-7)  could  be 
performed,.  This  displacement  would  introduce  a  term  of  order  e  on  the  right-hand  side  of  Equation  (4-8). 
Rake  would  be  equivalent  to  sweep  in  wings,  which  has  been  considered  by  Thurber7  If  warp  were  in¬ 
cluded,  a  displacement  in  y  would  be  necessary.  A  displacement  in  either  x  or  y  prior  to  stretching  would 
also  lead  to  a  transformed  Laplace  equation  with  a  term  of  order  e,  unless  the  displacement  were  constant. 
The  second-order  right-hand  side  means  that  the  first  two  nontrivial  solutions  are  governed  by  a  two- 
dimensional  Laplace  equation  with  the  spanwise  variable  a  parameter.  Experience  with  wing  perturbation 
problems  indicates  that  the  first  term  governed  by  Poisson’s  equation  is  a  practical  limit  to  the  number  of 
terms  in  the  series, 

The  potential  in  terms  of  the  inner  variables  must  also  be  expanded  in  a  perturbation  series.  The 
potential  in  the  inner  region  will  be 

0  (x,  y,  z;  e)  =  <p(ex,  ey,  z;  e)  =<b(x,  y,  z;  e) 

=  e$j(x,  y,  z)  +  v2(e)$  Jx,  y,z)+ •  •'  -  (4-9) 

The  inner  potential  is  of  order  e  because  the  velocities  in  the  inner  flow  are  required  to  be  bounded  as  e  ■*  0. 
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The  inner  boundary  condition  in  Equation  (2-10) 

n  *  V  0  =  -  gQ  •  n  for  r  €  Sg 

has  to  be  expressed  in  terms  of  inner  variables.  Instead  of  the  x  and  y  defined  in  Equation  (4-9),  it  is 
convenient  to  use  variables  measured  along  and  perpendicular  to  the  chordline  in  the  x,y  plane.  The  chord¬ 
line  is  at  an  angle  f>g(z)  relative  to  the  Jaxis  and  passes  through  the  reference  line,  the  z  axis.  The  variable 
£  j  measured  along  the  chordline  is  positive  in  the  downstream  direction,  and  the  variable  £ 2  is  normal  to 
the  £y  axis,  positive  in  upstream  direction.  This  coordinate  system  has  been  examined  in  Chapter  hand  is 
shown  in  Figure  6  below. 


Figure  6  —  Geometry  in  Plane  at  Fixed  z 
The  variables  £  j  and  £^  are  related  to  x  and  y  by 


*  =  £7  Sin  Vb(z)  -  h  cos  Vb(z> 
y  =  ~  £7  cos  vtfz)  -  \2  sin  I Pgfz) 


(4-10) 


Obviously  if  a  hub  existed  on  such  a  propeller,  a  complicated  expression  in  (x,y,z)  would  be  needed  for  its 
geometrical  specification.  For  narrow  blades  considered  here,  this  problem  is  not  significant. 

The  shape  is  assumed  to  be  given  relative  to  the  reference  surface  by  a  shape  function  (Chapter  1) 


$2  =  mj,z)  (4-11) 

Substituting  this  expression  into  the  previous  equations  for  x  and  yy  and  taking  the  appropriate 
differential,  one  can  find  a  normal  (Chapter  1) 
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With  these  relations  the  transformation  of  the  right-hand  side  of  Equation  (4-14)  becomes 


Qq  *  jY  =  ]/v2  + SI2  z2  I  sin  (<PB^$)  -  ~  (Oj,  z)  cos  (<pB~  j3)J 


+  e2  Sly 


a, 


d*B  3 E  3 E 


+  E  ■ 


dz  3 z  3£y  dz  / 

where  y  is  given  by  Equation  (4-7).  The  normal  in  Equation  (4-13)  transforms  to 

3  E 


~  =  W]*1*  ±e'0] 


d^B  bE  bE  d*B\ 

+  E — - I  k 


dz  bz  b%j  dz 
The  transformation  of  Equation  (2-10),  using  inner  variables  (Oj,  O2,  z),  then  becomes 


lib®  bE  3$ 
ey9a2  b$}  b0j 


d*B  bE  bE  d^B) 

dz  bz  b%j  dz  j 


3<T 

difB 

(  3$ 

39)  \ 

—  -f 

3z 

dz 

V  b<h~ 

°2  da2/_ 

-  }fu2  +  n2  z2  jsm  (<fi B-fy  -  cos  (<p B-fy | 


d*B  bE  bE  d*B\ 

-  _  — +  E - 

dz  dz  3  tjj  dz  J 


-  e2  Sly  (a^ -  -  —  +  E 


In  this  section  only  the  first-order  terms  are  considered.  Hence  for  <I>  =  e<J>y  the  equation  for  $>j  i 


b2<t>j  b2$} 

- + - =  0 

b  o  2  b  o2 

with  the  body  boundary  condition,  on  -  E(Oj ,  z),  0  <z  <R 


a$7 

3  a2  b boj 


bE 

Vj  (sin  cCj  --j—cos  <Xj) 


(4-16) 


(4-17) 


(4-18) 


(4- 19a) 


(4- 19b) 
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where 


dj(z)  =  Vtfz)  -  (i(z) 

The  remaining  boundary  condition  far  from  the  body  has  to  be  found  by  matching  with  the  outer  flow. 
This  process  will  be  formalized  in  the  following  section,  but  intuitively  it  is  expected  that  the  potential  in 

the  outer  flow  will  match  the  inner  flow  as  the  outer  variables  }/x ^  +  y ^  +  0,  0  <z  </?  and  the  inner 

variables  yOj^+o^  0  <z  <R.  In  the  outer  flow,  as  an  observer  approaches  the  reference  blade, 
the  flow  should  look  locally  like  a  two-dimensional  line  vortex  with  higher-order  effects  caused  by  the  in¬ 
duced  velocities  from  the  trailing  vortex  system.  Hence  for  the  first-order  inner  potential,  the  far-field 
boundary  condition  is  that  the  flow  look  like  a  line  vortex.  (Subsequent  matching  will  confirm  this 
hypothesis.) 

The  inner  potential  could  in  principle  be  found  from  an  inner  expansion  of  the  exact  potential 
given  in  Chapter  2.  Germain*^  reports  such  procedures  applied  to  wings  lead  to  considerable  difficulty 
and  that  it  is  simpler  to  solve  the  posed  two-dimensional  problem  of  Equation  (4-19).  In  fact,  using 
complex  variables  and  conformal  mapping  of  a  Circular  cylinder,  one  can  state  the  solution  immediately: 


I  -iaj  y. _ 

°2’  z>  =  V1  Re  y(al  +ia2>e  +a  f(Oj+i  o2) 

T  f(0]  +  io2)  _ia\ 

+  i  j—p  In - - -  -  (Oj  +  io2)e  1  j  (4-20) 


where 


00  an(z) 

f(o}  +  i  o2)  =  Oj  +  io2+  2  - 

n=0  (Oj+i  a2f 

is  the  inverse  mapping  function 
an(z )  are  complex  numbers  depending  on  the  section  shape 
a(z)  is  the  radius  of  the  transformed  circular  cylinder 
T(z)  =  4ir  a(z)V j(z)  sin  (cLj(z)  -  OLq(z))  is  the  value  of  the  circulation  when 
the  Kutta  condition  is  satisfied  and 
(trfz)  is  the  angle  of  zero  lift  of  the  profile  section 

Germain,  P.,  “Recent  Evolution  in  Problems  and  Methods  in  Aerodynamics,”  Journal  of  the  Royal 
Aeronautical  Society,  Vol.  71,  No.  682,  pp.  673-691  (1967).  < 
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Several  techniques  exist  for  finding  the  mapping  function  for  a  given  section  geometry.  Most  of 
these  are  discussed  by  Thwaites.*^  It  is  not  necessary  at  this  time  to  examine  the  techniques  for  obtaining 
these  coefficients,  but  they  are  considered  known  in  terms  of  the  geometry.  Specifications  of  the  circu- 
lation  renders  the  two-dimensional  problem  unique.  In  a  review  of  the  planar  wing  problem,  Ogilvie,^ 
points  out  that  the  inner  expansion  can  be  taken  as  a  Laurent  series  but  doing  so  adds  no  new  information 
or  insight  to  the  direct  solution  from  conformal  mapping  given  previously. 

Although  the  solution  indicated  in  Equation  (4-20)  is  the  required  solution,  in  anticipation  of  the 
matching  discussed  in  the  next  section,  the  form  of  the  solution  at  large  distances  from  the  body  will  be 
examined.  Figure  7  shows  complex  number  Oj+i(J2 ,  given  by  fxeI0J 


Figure  7  —  Profile  Coordinates  for  Inner  Flow 


Hence 


co  =  tan 


The  expansion  is  straightforward,  except  for  the  log  term 


c  o 

Thwaites,  B.,  “Incompressible  Aerodynamics,”  Oxford  University  Press,  England  (1960). 

^Ogilvie,  T.F.,  “Singular  Perturbation  Problems  in  Ship  Hydrodynamics,”  University  of  Michigan, 
Department  of  Naval  Architecture,  Report  096  (Oct  1970). 
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f(H^)  iieiw  (  a0  al 

In  l— - 1  =  In  - +  bill  +  —  + - — -  + 

a  a  \  jueICO  ^2e2/oj 


a  V  a2 


Hence  for  large  ji ,  the  expansion  is 


QjfO],  o2,  z)  (  Qj 

- - - -  ~  -  2  a  sin  (cLj-cCq)  (gj+  ojq(z))  +  Re 

a?  .  a2el(ar^  aQ  a0~al 

+  —  e~i( tit 2(jj)  +  0  .  m  + -  j - .f - -f  •  »  • 
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a0  arl/2a0z 


+  2  i  a  sin  (di-CL))  — ; — + - — — 

lieiU)  n2e2ico 


a2-ll2aQai  + 


n3  e3i“ 


where  c orfz)  is  arbitrary. 

By  grouping  terms,  one  finds 


—•=  -  2a  sin  (dj-dg)  ^co  +  dJgfzjj  +  2  (2j  kn  (al’  z) 
where  kj(dp  z)  =  a  cos  (cCj-u)  +  Re  j-  +  2i  sin 


(4-21) 


k2(alfz)  =  Re  I—  e-Wf203*  +  a^r203)  +^~  sin  (dj-dg)  e  2l0J  (arl/2aQ2  )\ 

la 2  ■ 
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kjfdj,  z)  =Re  < 


ao  a0  al 

_  e-i(a}+3u)+ - e-^r3^) 

3  a  - 


2i 


+—  sin  (cCj-Oq)  e  3l0)  fa2-ll2a(p1+ll3a(f) ' 


etc. 

The  first  term  is  a  vortex;  the  second,  a  dipole.  Then  the  higher-order  terms  are  listed.  No  source  terms 
appear,  as  expected,  since  none  are  in  the  first-order  outer  solution. 


Matching 

The  two  expressions  for  0  are  the  leading  terms  of  asymptotic  expansions  valid  in  different  domains 
of  the  flow  field.  For  the  reference  blade,  the  outer  solution  is  valid  everywhere  in  the  flow  field,  except 

z</?,  Vx2Vy2  <CQ.  The  inner  solution  is  valid  for  z  </?,  }/x2  ~+  y2  or  }/x2  +  y2  ^Cj  W.  For 
an  overlap  domain  to  exist,  one  must  assume  that  the  domain  of  validity  of  the  outer  solution  can  be  ex¬ 
tended  to  y  x2  +  y2  >  CV)  (e)  where  lim - -►  w,  This  extension  of  the  region  of  validity  is  discussed 

0  _  Ctfe) 

e+0  v 

by  Kaplan^  and  Cole.^  In  this  assumed  overlap  region,  the  x  and  y  variables  are  taken  as 


*  =  x  gfe) 


y  =  y  gffl 


(4-22) 


where  x  and  y  are  constants  of  0(1),  and  gfe)  is  the  function  which  places  x  and  y  in  the  intermediate 
region.  To  insure  x  and  y  are  in  the  intermediate  region,  gfe)  must  satisfy 


gfe)  -  0 
e 

- -* 0 

gfe) 


as  e  -*0 


The  x  and  y  are  the  outer  variables,  and  the  inner  variables  are  obtained  from  Equation  (4-7)  as 


(4-23) 
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/  rr 

From  the  definition  of  an  asymptotic  sequence,  e.g.,  Erdelyi,  the  two  solutions  can  be  matched  to  order 

Hrfe)  if 


N  N  _  _£ 

L  Sn(e)  <t>n(xg,  yg,  z)  -  Z  %<*)  %(x  e.iy  l,  z) 

n=l  n=l 

lim  - - -  =  0 

e-0  Vrfe) 


(4-24) 


(N  =  1,2  "  •) 


Hence  the  dominant  terms  can  be  matched  if 


lim 

€->0 


bjfej  <pj(xg,  yg,  z)  -  e  <b,(x  e,y  t  T  z) 

- -  0 

(J-jfe) 


A  match  to  order  e,  i.e.,  fij(e)  =  e,  is  possible  by  selecting  5 j(e)  =  e  and  by  requiring  the  following  equality 
to  hold 


lim  <t>j(xg,  yg,  z)  =  <Sy°°  z)  (4-25) 

g+0 

The  matching  principle  given  by  Van  Dyke  expresses  the  same  requirements  as  described  previously  in  a 
more  concise  form. 

The  matching  principle  will  be  used  later  but  for  now  a  much  simpler  argument  can  be  used  to 
find  the  circulation  in  the  outer  flow.  The  argument  is  that  the  spanwise  circulation  must  be  independent 
of  the  chordwise  stretching. 

Thus 


(z)  =  4ira(z)  Vj(z)  sin  (aj(z)  -  arfz»  (4-26) 

Hence  the  outer  solution  is  completely  determined  to  the  first  order.  There  are  no  incompatibilities;  conse¬ 
quently,  the  first-order  solution  is  complete. 
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Erdelyi,  A.,  “Asymptotic  Expansions,”  Dover  Publications,  Inc.,  New  York  (1956). 
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A  composite  solution,  not  necessarily  unique  but  uniformly  valid  throughout  the  flow  field,  can  be 
directly  written  down.  To  within  an  arbitrary  constant  of  z,  which  to  second  order  is  only  a  parameter  in 
the  inner  flow,  the  first-order  composite  solution  is  the  sum  of  the  inner  solution  plus  the  outer  solution 
minus  the  common  part.  The  common  part  is  either  side  of  the  equality  in  Equation  (4-25),  and  from  the 
inner  solution  in  Equation  (4-21)  it  is  seen  to  be 

0 C'P ,  =  €  <^°°,  z) 

0  z  >R 

-  e  V j(zj  2a(z)  sin  (dj(z)  -  d^z))  corfz))  0  <z 


Hence  the  composite  solution  is 


(e$j(x}  y,  z) 

4*c  e  0/x,  y ,  zj  +  e&jfoj,  o2,  z) 


z  >/? 


0  <z  <R 


(4-27) 


y  +  2  e  V ja  sinfctj-CiQ)  (oj+ 0Jq)  ) 

Use  of  the  matching  principle  will  show  that  co^(z)  is  such  that  co  +  0Jq  =  tan^y/x.  This  is  shown  in 
Equations  (4-34)  and  (4-45).  Such  a  composite  solution  is  useful  in  fmding  information  in  the  flow  field 
such  as  field-point  velocities.  Since  this  overall  view  of  the  flow  field  is  not  required  in  the  present 
analysis,  it  is  not  considered  again;  however,  similar  expressions  valid  to  higher  order  can  be  obtained  just 
as  easily  when  the  individual  higher-order  solutions  are  known.  Composite  solutions  are  discussed  at  some 
length  by  both  Van  Dyke^  and  Cole.^ 

The  circulation  to  first  order  corresponds  to  what  is  known  as  blade-element  theory. ^  Its  use¬ 
fulness  is  restricted  to  narrow  propellers  with  light  loads.  The  second-order  terms  derived  in  Section  B  will 
modify  the  circulation  distribution  with  higher-order  effects  and  will  extend  the  usefulness  of  the  solution. 


B— SECOND-ORDER  TERMS 


Outer  Flow 

The  second-order  terms  in  the  outer  flow  will  arise  from  three  sources.  The  first  will  be  from  the 
integration  over  the  blade.  Although  this  could  be  extracted  by  a  Taylor  expansion  from  the  exact 
formulation  in  Equation  (2-50),  it  is  simpler  to  assume  that  a  distribution  of  dipoles  in  the  x  and  y  direction 
will  be  the  required  singularity  system.  Once  again,  if  the  assumed  form  is  not  correct  then  the  matching 
will  not  be  possible.  The  second  term  in  the  outer  series  will  be  from  the  second-order  term  in  the  cir¬ 
culation  expansion  in  Equation  (4-1).  This  will  appear  in  an  equation  similar  to  Equation  (4-5).  The  last 
second-order  contribution  will  appear  because  of  the  first-order  modification  to  the  position  of  the  vortex 
sheet.  Since  this  last  term  is  too  complex  to  treat  analytically,  an  approximation  will  be  made  to  simplify 
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where  the  second  term  is  left  in  the  inconsistent  form  because  of  the  approximation  to  be  made.  Since 
there  is  little  hope  of  analytically  integrating  Equation  (4-29),  one  seeks  a  reasonable  approximation.  The 
main  contribution  of  the  modified  position  of  the  shed  vorticity  will  be  from  points  close  to  the  line. 
Hence  a  Taylor  series  expansion  of  the  velocity  is  called  for;  however,  for  even  greater  simplicity,  all  but 
the  first  term  is  ignored,  that  is,  the  velocities  are  evaluated  at  the  lifting  line.  Hence 


h  ~zo 


u+euJl>(z0) 

- - - (e-eb) 

Slz0  +  e  up)  (zQ) 


U  +  ~e  uj1)  (zq) 

The  quantity  Zq - 

nz0  +  e  uj1)  (z0) 


1 

is - times  the  pitch  of  the  shed  vortex  at  the  lifting  line.  This 

2  7T 


pitch  is  often  called  the  hydrodynamic  pitch  and  it  is  denoted  by  P^Zq).  Hence 


w 

h  a—  i 'e~eb>  <«®> 

It  is  possible  to  propose  a  slightly  better  model  for  the  integral  in  Equation  (4-29)  since  the  induced 
velocity  both  at  the  lifting  line  and  the  value  far  downstream  are  known  (the  value  for  x  00  is  twice  that 
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at  the  lifting  line).  Numerical  calculations^  have  been  performed  for  which  the  information  necessary  to 
propose  a  variation  in  the  induced  velocities  as  a  function  of  x  or  d  could  be  obtained.  However,  the 
previous  assumption  is  universally  made  and  considerably  simplifies  the  calculations. 

A  similar  universal  and  significantly  simplifying  approximation  is  to  ignore  the  contraction  of  lines 
of  constant  circulation.  This  means  the  radial  velocity  in  Equation  (4-28)  is  approximated  by  zero.  Any 
other  constant  approximation  for  it  would  lead  to  negative  values  for  the  radial  position,  clearly  an 
unacceptable  occurrence.  Here  too,  alternatives  to  this  crude  approximation  would  be  the  numerical  cal¬ 
culations  or  results  of  actuator-disk  theory.  These  alternatives  produce  complications  not  essential  to  the 
continuation  of  the  present  investigation,  and  we  accept  the  approximations  as  customarily  used. 

These  two  approximations  involved  in  the  higher-order  description  of  the  shed  vortex  position  are 
defined  as  the  basis  of  moderately-loaded  propeller  theory.^ 

Hence  to  the  first-order,  the  approximate  position  vector  is 


w 

i  =  zo£Z+ - (O-Qb)i- 


(4-31) 
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'lo  +  eh 


u  W 

Since  this  approximation  merely  replaces  the  constant  —  in  Equation  (4-3)  by -  ,  no  further  insight  is 


ft 


27T 


added  to  the  problem  by  considering  this  term  and  from  now  on  it  will  be  ignored.  That  no  essential 
modification  to  the  theory  results  is  not  obvious  in  the  formulation  of  the  velocity  potential  but  can  be 
seen  in  the  calculation  of  induced  velocities  by  the  Biot-Savart  Law;  see  Appendix  B.  However,  the  cal¬ 
culation  of  induced  velocities  is  usually  done  with  expressions  derived  for  constant  P-. 

Thus  only  the  modification  of  the  blade-element  circulation  distribution  and  the  dipole  distribution 
comprise  the  direct  second-order  outer  potential.  Although  a  general  gage  function  must  assumec* 

for  the  circulation  distribution,  second-order  terms  from  the  blade  integration  in  Equation  (2-50)  are  easily 
—2 

seen  to  be  order  e  .  (Change  the  chordwise  integration  to  inner  variables  and  expand  the  expression  in  a 
Taylor  series.)  Thus  the  two-term  outer  solution  is 


*  cn 

Notable  exceptions  are  Erickson  et  al.  and  Cummings,  who  calculate  vortex  interactions  for  lifting¬ 
line  theory. 

•^Erickson,  J.E.  et  al.,  “A  Theory  for  VTOL  Propeller  Operation  in  a  Static  Condition,”  Curtiss-Wright 
Corporation,  Caldwell,  N.J.  (Oct  1965). 

jjc 

The  vorticity  vector  A_=  n  X  V  T  has  no  additional  term  from  the  e~  component  of  n,  since  e~x  V  T 
=  ft 
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4>  (r;  e)  =  etpj(r)  +  b^e)  <p12(r)  +  e2  tyr) 
where  (j>j  is  given  by  Equation  (4-5),  and 

j  rR 

*!&=—  r (2)  (z0)G(r,z0)dz0 

4lT  JO 


(4-32) 


/  f* 

*Ji)  =—  L 
4«b=0  L 


xnj(z0)  +  (y  +  zQ  sin  6b)  \i^zQ)  dz 


0 


\  x2  +  (y  +  zQ  sin  V2  +  (  z  “  z0  cos  6^  2 


3/2 


(4-33) 


where  r(2)  (zq)  is  the  next  term  in  the  circulation  expansion  in  Equation  (4-1),  and  [Lj(z)  and  O-^zq)  are 
the  dipole  strengths.  All  three  of  these  are  unknown. 

To  match  with  the  inner  flow,  one  needs  the  asymtotic  sequences  which  results  when 
Equation  (4-22)  is  substituted  into  the  two-term  expression  for  (j).  By  an  indirect  procedure  described  in 
Appendix  B,  the  expansion  of  <p  is  found  to  be 

tfpc,  gy,  z;  e)  =  Frfz;  e)  +  52  tan”1?- 

+  egfx  uJV  (z)  +  y  u<D  (z)J  (4-34) 

t  2  xfij  (z)  +  yn2  (z) 

+ - +H.O.T. 

4ltg  =2  =? 
xz  +yz 

where  uj^  (z)  and  uj^  (z)  are  the  axial  and  tangential  velocities  at  the  lifting  line  as  induced  by  the 
first-order  circulation  distribution. 

Inner  Flow 

For  the  second  term  in  the  inner  flow,  the  governing  equations  must  be  examined.  From 
Equation  (4-18),  the  boundary  condition  on  the  body  can  be  determined.  One  substitutes  the  two-term 
inner  expansion 


<P(X,  y,  z;  e)  =  6<J>7  (Op  o2;  z)  +  v2  <$>2(Op  o2;  z) 


(4-35) 


into  that  equation  and  finds  that  $2  satisfies 
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v2  b*2  2 

7  [l72- j 


(4-36) 


O2  - E(02 ;  z),  0  <z  <R 


Obviously,  it  is  most  convenient  to  take  lim  eJ ~  0  so  that  the  body  boundary  condition  becomes 

e+0 


3 $2  9£ 

do  2  3 


( o};z )• 


(4-37) 


on  o 2  -  1  *  z)>  0  <2  </? 

Similarly,  in  the  flow  field,  Equation  (4-8)  requires  $2  to  satisfy  a  two-dimensional  Laplace  equation.  As 

yet  the  boundary  condition  at  +  o  ^  +  00  is  unknown  and  must  be  determined  by  matching.  The 

simplest  problem  for  which  the  two-dimensional  flow  cannot  pass  through  the  body  is  that  associated  with 

an  airfoil  at  an  angle  of  attack  to  the  flow.  Hence  as  o  ^  +  0 ^  *  00  one  assumes  that  the  potential 
behaves  like 


*2(°i’  °2  ;  z)  ~  Oj  V^z)  cos  a^z)  +  O2  V ^z)  ty2) 


(4-38) 


where  and  ^(z)  are  unknown.  This  two-dimensional  problem  is  also  easily  solved,  using  complex 
variables  and  conformal  mapping.  With  the  same  notation  as  in  Equations  (4-20)  and  (4-21)  the  second- 
order  inner  potential  is 


-=  Oj  cos  <12  +  02  sin  CL2  -  2  a  sin  (C^  -  &q)  (to  *  oo^  (z)) 


(4-39) 


where  cotfz)  is  arbitrary,  and  the  z)  are  given  in  Equation  (4-21).  In  order  to  find  the  necessary 
form  of  this  expression  for  matching,  one  substitutes 


°]  ~  °]  s/e  and  a2  =  o2  g/e 


(440) 
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Hence,  the  matching  is  reduced  to  the  expression 


r<2>  tan  1  6  gfx  uW  (z)  +  y  uj ^  (z)j  + 


(Oj  cos  &2  +  s*n  *h) 


2a  sin  (C^-dg)  (oi+  c Oq)  +  0  • 


(444) 


where  n^e)  is  such  that  terms  not  explicitly  given  vanish  in  the  limit. 
This  matching  is  possible  by  selecting 


52  =  vje)  =  e' 


(445) 


co  +  oj  n  —  tan  1  f/x 


for  which 


x  ujl)  (z)  +  y  (z)  =  V2  [Oj  cos  0^  +  s*n 


(446) 


tf2)  (z)  =  4ira(z)  V^z)  s*n  (°2  ~ 


(447) 


To  determine  the  performance  of  the  propeller,  the  circulation  distribution  must  be  known  ex¬ 
plicitly.  That  is,  both  V^z)  and  tf^z^must  be  determined  in  Equation  (447).  To  find  these  terms  we 
use  Equations  (446)  and  (4-10); 

OjfuJ n  sin  Vg-uj1)  cos  1 Pg)  -  ojuj1*  cos  $g  +  uj1)  sin  <pB) 

=  Oj  V2  cos  $2  +  °2  ^2  s*n  **2 


V2  =  uJW  *  u^2 


(448) 
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(449) 


ujl)  cos  \p q  +  u{V  sin 

tan  $2  ~  ~  - - : - — - 

uj^  sin  -  uj1)  cos 

Thus  the  second-order  solution  is  complete.  In  the  next  section,  formulas  for  the  propeller  performance 
are  given. 

C— PERFORMANCE 

Determination  of  the  circulation  from  the  given  profile  coordinates  shows  that  the  singular- 
perturbation  problem  gives  a  quadrature  for  the  performance.  However,  the  appropriate  form  of  the  inte¬ 
gral  for  the  thrust  and  torque  has  not  been  determined  yet. 

The  expression  for  thrust  is  given  by  Equation  (2-57) 


T  =  Z 


I. 


'B 


[2  q0  +  V<t>)  •  V  0/  /  •  nds 


In  lifting-line  theory  the  expression  for  thrust  is  usually  obtained  from  the  outer  flow.  Here  we 
derive  the  appropriate  expression  from  the  inner  flow  since,  in  the  preceding  equation,  the  inner  potential 
is  the  appropriate  one  to  use.  Now 

I  -bE  , 

1  °  nds  =  i_'  y±  e_2  f  d^j  dz 


bE 

=  +  e\cos<pB  +-^-sin  <PB)  do j  dz 


Qq  =  +  &  z^  [cos  (<pB  -  P)ej  +  sin  (q>B  -  j3)  e2j  +  Q}’k 

=  V jfz)  [ cos  (tj  ej  +  sin  (tj  +  & 

V<p  =  V$(o]t  o2,  z;  ej 

■  ?  h, + <2  y + f 2  *2  m  k 


see  Chapter  1  for  the  gradient  in  inner  variables. 

Thus  to  the  first-order  in  e,  the  integrand 


(2q0+  V  0)  •  V  <t> 

involves  only  two-dimensional  quantities,  i.e.,  z  is  a  parameter.  By  inspection,  the  integrand  can  be  changed 
to 
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since  \qg 


is  a  constant  with  respect  to  the  chordwise  integration  and,  hence,  integrates  to  give  zero. 


In  two-dimensional  wing  theory,  one  defines 


k 

k 


+  (2o 


*(So 


2D 


2D 


+  do}  =  -jFN 


o  dE  2 

*  "2  ***!(,  1  ’  c 


(4-50) 


where  FN  and  Fc  are  the  forces  normal  to  and  along  the  chord.  If  one  defines 

U  +  eu  (1)  +  •  °  • 

fy  =  tarC 1 - 

J2z  -  (euf1)  +  *  *  •) 

then 

Fn  =  L  cos  (<PB  -  fy) 

Fc  =  L  sin  (<pB  -  fy) 

where  L  is  the  lift  of  the  blade  section. 

Hence  the  thrust  coefficient  becomes 
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2  Z 


nU2  R2  p 


+  sin  ypB  sin  (<pg  -  p.)]  dz 


fR 

I  eL  cos  Pj  dz 
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The  lift  is  given  by 


L  =  p  }/(U  ~+  e  uJ2)  ~+  ~  •  •j2  +  ~{Ez  -  (e  ~uJ2)  ~+  •  •  Y(z;  e) 


(4-51) 


(4-52) 


(4-53) 


Hence 
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CHAPTER  5 
DISCUSSION 


A-COMMENTS  ON  RESULTS 

The  lifting-line  analysis  is  essentially  the  same  as  that  given  in  the  literature,  except  for  notation. 
One  apparent  difference  is  that  the  two-dimensional  sections  are  defined  on  sections  cut  by  a  plane  per¬ 
pendicular  to  the  straight  blade-reference  line.  However,  a  Taylor  expansion  of  the  shape  defined  on  a  cy- 
Under  shows  that  to  second-order  they  are  the  same  since  the  difference  z  -  co  is  0(e  ).  Of  course  higher- 
order  terms  would  show  the  effect  of  the  new  section  definition. 

The  Ufting-line  analysis  was  developed  for  propellers  with  a  straight  lifting  line.  If  the  blade- 
reference  line  were  raked  or  warped,  the  outer-flow  potential  could  be  easily  found.  The  overall  solution, 
however,  would  not  be  as  given  in  Chapter  4,  because  the  second-order  problem  would  include  solution  of 
Poisson’s  rather  than  Laplace’s  equation.  The  appropriate  solution  for  such  cases  could  probably  be  ob¬ 
tained,  but  the  lifting-line  analysis  is  appropriate  only  for  propellers  with  narrow  blades  and  hence  is  not 
useful  for  marine  propellers.  For  marine  propellers,  the  most  important  analysis  is  the  lifting-surface 
formulation,  and  the  rest  of  this  discussion  will  concern  this  problem. 

In  the  preceding  chapters,  the  expressions  for  velocity  were  left  in  the  form  of  the  gradient 
operating  on  an  integral.  This  is  because  they  are  principal-value  integrals,  and,  if  the  differential  would  be 
taken  inside  some  of  the  integrals  (Equation  (2-19)),  they  would  not  exist  in  the  mathematical  sense;  see 
Tychonov  and  Samarskii  The  singularities  which  result  when  the  operator  is  taken  inside  the  integral  are 
usually  called  Hadamard  singularities,  and  as  shown  in  Equation  (2-19),  the  finite  part  of  the  integral  is 
defined  in  terms  of  limit  operations.  However,  if  the  expression  for  q  is  integrated  by  parts,  the  order  of 
the  singularity  will  be  reduced  for  the  limit  r  -*s,  and  the  integrals  will  exist  as  Cauchy  principal  values. 

(For  a  principal-value  integral  the  region  surrounding/^  must  be  excluded  from  the  integral,  and  the  limit 
with  the  maximum  dimension  of  the  excluded  region  going  to  zero  must  be  taken;  for  a  Cauchy  principal- 
value  integral,  restrictions  are  placed  on  the  shape  of  the  region  excluding/^). 

In  evaluating  the  expressions  numerically,  either  formulation  of  the  problem  is  acceptable  since  both 
expressions  involve  a  region  surrounding  the  singularity  which  requires  special  consideration.  In  fact, 
Sparenberg,^  Pien  and  Strom-Tejsen,^  and  Tsakonas  and  Jacobs^*  treat  the  higher-order  singularity,  while 
Pien,*^  Kerwin,^’^  Yamazki,^  and  Murray*^  work  with  principal-value  integrals.  Sparenberg"^ 
presents  an  analysis  for  both  points  of  view. 

B— EXPRESSIONS  FOR  VELOCITY  AS 
INTEGRALS  OF  FIRST-ORDER  SINGULARITIES 

Equation  (3-21)  for  the  perturbation  potential  consists  of  three  terms 

<t>=  (j)s  + (pT+ <j)v  (5-1) 


76 


In  this  equation  <p$  is  an  integral  over  the  surface  of  the  blades  corresponding  to  the  lifting  aspect  of  the 
propeller 


<t>s(l)=^-  V  Jf 

4t  b=0  "  Sn+ 


Lb^I’  ^  —Os  ’ 


d%j  d<Z 


and  (fry  is  also  an  integral  over  the  surface  of  the  blades;  however,  due  to  the  blade  thickness 


L  ff  +  oj2  £L2  cos  (ipp-\ 3)  bEj 

2ir  b=o  JJ  sD+  ;  ;  — 


d%j  doS 


and  <j>v  is  an  integral  over  the  shed  vortex  sheet 


Z~~l  R  ww 

<t>v=j  L  [  r (z0)dz0  f  n£,  •- 
471  b=0  J0  J  j r 


+  Ho 


where  F g(cj(Z3)f  73)  «  T(Z o) 

Nqs  is  given  in  Equation  (3-10)  and 

Nqv  is  given  in  Equation  (3-18) 

In  the  following  work,  the  expressions  will  be  considered  without  the  effects  of  rake  or  warp  since 
they  are  less  complicated  and  permit  comparison  with  results  in  the  literature.  The  formulation  with  rake 
and  warp  can  be  obtained  similarly. 

In  design,  the  detailed  geometry  of  the  propeller  is  to  be  determined.  Since  this  is  so,  the 
position  of  the  trailing  edge  can  be  chosen  to  simplify  the  integrals.  With  neither  rake  nor  warp  the  blade- 

2tt 

reference  lines  are  the  6  ^  values  (Q  ^  -  —  b),  If  the  trailing  edge  of  the  blade  is  chosen  to  lie  on  the 
helical  surface  passing  through  this  line,  then  in  Equation  (3-17) 


Io~zO~cS+  o^~  ®  b^l’  ®  ^  ®  TE(zq) 
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where  /  is  given  in  Equation  (B-10)  of  Appendix  B.  In  the  preceding  integral,  Zq  is  a  dummy  variable  for 
Z5.  As  shown  in  Appendix  B  this  form  is  the  same  as  that  found  in  the  literature,  and  hence  the  detailed 
explicit  form  used  for  the  numerical  analysis  is  not  repeated  here. 

When  the  pitch  of  the  blade-reference  surface  and  the  shed  vortex  sheet  are  the  same,  Pien^ 
found  that  significant  savings  in  computing  time  could  be  obtained  by  rearranging  the  limits  of  the  infinite 
integral  so  that  it  would  need  to  be  calculated  only  once.  Naturally  this  left  an  integral  over  a  finite 
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region  to  be  evaluated  for  each  new  position  on  the  blade.  For  points  everywhere  on  the  blade  except  the 
trailing  edge,  the  preceding  expression  is  a  straightforward  quadrature  without  singularities  in  the  integrand. 
When  Pien’s  procedure  is  used,  Cauchy  principal-value  integrals  occur,  but  apparently  these  are  easier  to 
handle  than  the  repeated  infinite  integral. 

The  preceding  expression  for  the  induced  velocity  is  usually  modified  to  include  the  moderately- 
loaded  assumption  discussed  in  Chapter  4  which  consists  of  replacing  the  expression  by  the  approximate 
one  for  +  e  $7,  given  hi  Equation  (4-31),  which  is  derived  from  lifting-line  theory.  Although  an  argu¬ 
ment  can  be  made  for  this  approximation  in  lifting-line  analysis,  no  justification  is  given  for  its  use  in 
lifting-surface  theory.  In  any  event,  only  a  portion  of  the  second-order  effects  are  considered  with  this 
modification;  see  Section  D  of  Chapter  3. 

The  expression  for  the  potential  due  to  blade  thickness  can  be  written 


j  Z-l  - R  y -  r-p w/ 

4>t(l)  = - £  Vu2  +  ZS2  cos(v  -P)du  •  _ '_ZL  d$j  (5-10) 

2l*  b=0  20  Jr-/7Xi  t-i/il 


cl< $  dEn 


yco)  lr-io 


where 


t-Sq  =  (x~%j  sin 


VP)L+  +  ^ ««  ^ 


%1  cos  v?  \\ 


+  {z-  ZS  cos  ^  6h  + 


%l  cos  V 


b  ~ 

CO 


Hence  qj~  V  0^  is  given  by 

Z-i  rR 


dEn 


i  rK  / - : 

q— sz  £  I  f  +  CO2  £1* 

2ir  b=0  J0 


_  fCj 

°  I  i _ 


cos(<p-P)d(5‘ 

J  C2  t-l/ 


d£j  (5-11) 


For  points  on  the  blade  r  ->Sq,  only  a  principal-value  integral  must  be  evaluated. 

Since  this  expression  is  the  same  as  that  found  in  the. literature,  it  is  convenient  to  make  the 
necessary  changes  to  demonstrate  the  equivalence.  The  changes  are  that  the  blade  surface  is  described  by 
co,  6  instead  of  £7,  <o,  and  that  pitch  P(Z5)  instead  of  pitch  angle  *Pp(<3)  is  used.  This  is  the 
description  used  by  Kerwin.22,2^  In  addition,  the  term  cos  (ipp—fi)  is  ignored,  and  constant  pitch  is 
assumed  since  the  blade  normal  is  taken  to  be  e ^ 

For  constant  pitch  and  neither  rake  nor  warp,  the  only  component  of  velocity  in  the  boundary 
condition  of  Equation  (3-27)  is 

v“V0°£2  (5-12) 
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P(o5)  Q-P(r)  <p,.  P(  <5)  _ 

-  *  /  \  OJ1 


co  sin  (  6b  +  6  -  <f>) 


P(r)  <$>-P(io)  (A2  2  ~7  ~  ,  3/2 

- I  +  +  cu  -  2  00  r  cos  (  6^  +  6  -  <t>) 


(5-15) 

Except  for  notation,  this  agrees  with  Kerwin’s  expression.  Kerwin  uses  2Ej  for  his  thickness  dis¬ 
tribution,  so  there  is  also  a  difference  factor  of  2.  When  the  pitch  is  constant,  Gj  is  a  function  of  (6  -  0). 

The  remaining  term  to  be  examined  is  the  integral  over  the  lifting  surface.  The  velocity  from  this 
component  of  0  is  given  by  ■  .  * 
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«S=V<1>S 


1  Z'1  ff  n  (so*m(L-lo)\ 


(5-16) 
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=—  l  W  - 

4n  b=0  "  SD+  |  r 


i*  no  • 

- r-  3 - - - (r-s0)  >  ds 

-irtl5  V-snf 


When  r  -  Sq  *♦  09  the  singularity  in  this  expression  is  second-order  rather  than  first-order  as  in  principal- 
value  integrals;  see  Equation  (2-19). 

To  integrate  by  parts,  the  gradient  of  the  kernel  in  Equation  (5-16)  must  be  expanded  and 
rearranged 

In  •  (z~s)\  r-s_  f  _r~A  \ 

V  ( -  I  =  (n  0  V  J  — ■ —  -f  n  X  I  V  „X -  1 

J  "  t-I?  '  '  ImI* ) 

But  for  an  arbitrary  vector  A, 


nX(VXA)  =  (nXV)XA  -  (n  •  V)A  -n(V  *A) 


Now,  since 


v  *  - 1\  =  0  (1-1*0) 

l  \r-s\3 


the  preceding  expression  becomes 


(n  •  (r-s)  \  r-s  r-s 

- —  WflXVX - -  =  -  foX  VX““~7 

|r-s|^  /  Ir-sl5  Ir-sl^ 


To  actually  perform  the  integration  by  parts,  we  determine  that 


(n X  V X  rB - -  =  rB fox  vx  — -  +  fo x  vs  rfi;x- 


Hence  the  velocity  is 
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(5-17) 


^xvT*iv/r 


One  of  Stokes’  theorems 


jfj  (U  x  V)  X  A  ds  =  jc  di  X  A 


where  c  is  the  curve  bounding  s ,  and  di  is  along  the  curve  in  the  usual  orientation.  Thus  the  expression 
for  velocity  is 


SB*0\  /  l-'-«h 


_  f  r  r»  x- 


Now  r5  is  zero  at  the  blade  leading  edge;  and  ,  at  the  trailing  edge.  The  arc  element  along  the 

blade  outline  is 

di 

di  - - dco 

d  co 


Hence 


=  t  do5 
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b=0  I  JJ  Sn+ 


vs  *VX 


r*  i-/ 

IYVIX — 77  *o> 
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(5-18) 
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The  line  integral  cancels  a  similar  contribution  from  the  trailing  vortex  sheet  in  Equation  (5-9),  and  the  sur¬ 
face  integral  is  the  Biot-Savart  Law  with  a  vorticity  distribution: 

A  =  nXVsrB 

=  nX(q+  -  <f) 

Since  the  conversion  is  general,  it  would  also  be  possible  to  express  velocities  from  the  nonlinear 
results  of  Chapter  2  in  terms  of  a  vorticity  distribution.  To  interpret  the  results  in  such  a  manner,  it  will  be 
necessary  to  define  an  inner  surface  of  the  blade  which  has  zero  flow  velocity.  Then  the  vorticity  will  be 
simply  nX  V0.  Since  no  further  insight  is  provided  with  this  approach,  it  is  not  pursued  further. 

The  explicit  expressions  for  the  linearized  velocity  in  the  case  of  the  general  propeller  geometry  can 
be  found  in  a  straightforward  manner.  Since  the  results  will  eventually  be  compared  with  equations  found 
in  the  literature,  the  geometrical  restrictions  appropriate  for  those  cases  will  be  made.  These  are  that  the 
propeller  is  neither  raked  nor  warped  and  that  the  pitch  is  constant.  For  constant  pitch 


n  =  £2 


and  for  neither  rake  nor  warp 


«xvr5  =  e2x 


arB  cos2  <pp  9r 


dco  co 


A=  —  sar-  —  + 


drB  Zjcos2  <pp  arfi 


3co  co  3 


Now  since 


"b  =  f  Tto  29  dTl 

J  c2(co) 


one  finds 


A- 1*0 -(4/ 

Vaco*' 


3  r^l  £jCos2<pp  \ 

—  I  ydri  + - — - 7 1  ej 


(5-19) 
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This  is  the  form  of  the  vorticity  vector  required  in  the  coordinate  systems  developed  here.  However,  in  the 
literature  it  is  more  common  to  work  with  cylindrical  polar  coordinates  (co,  0,x).  In  this  system  the 
parametric  representation  of  the  helical  blade-reference  surface  is  given  by 


x  =  co  tan  @  =~Z~® 

P  2ir 


y  =  -  co  sin  6 


z  —  <jJ  cos  6 


(5-20) 


For  the  surface  coordinate  system  u\iP  -  to,  8 ,  the  metric  coefficients  are 


\2 

]  +  or 


Hence  for  constant  pitch,  the  gradient  of  Fg  in  this  coordinate  system  is 


Now  since 


ar*  si  ^ 

nXvrB-e2X  e%  +  -  9(9 


(5-21) 


one  finds 


and  hence 


r t-  J 


eLEf^ 


nXvVB  =  e2X  7£ 


+  up  7 d'Q  j  e. 


This  is  the  same  expression  obtained  by  Kerwin,^  although  he  used  the  result  in  a  different 


(5-22) 


manner. 
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Explicit  formulas  for  calculating  the  velocity  from  the  Biot-Savart  Law  are  given  in  the  next 
section  and  hence  are  not  repeated  here. 

C—COMPARISON  WITH  OTHER  FORMULATIONS 

In  the  past,  linearized  propeller  lifting-surface  formulations  have  been  developed  from  three  view- 
points-a  postulated  vorticity  distribution  (Pien,^  Kerwin,^’^  and  Murray*^),  an  acceleration  potential 
(Sparenberg,^  Hanaoka,*^  Pien  and  Strom-Tejsen**1),  or  a  velocity  potential  (Yamazaki^'^4).  Both 
velocity  potential  and  acceleration  potential  formulations  are  based  on  the  equations  of  motion,  while  the 
vortex  distribution  formulation  is  based  on  laws  derived  from  the  governing  equations.  Of  course,  all  so¬ 
lutions  should  be  equivalent  so  the  development  is  only  a  matter  of  personal  preference. 

The  details  of  the  formulation  of  the  velocity -potential  function  are  not  given  by  Yamazaki.  How¬ 
ever,  examination  of  the  potential  reveals  it  to  be  compatible  with  the  velocity  potential  derived  from  the 
acceleration  potential.  Such  solutions  are  given  in  terms  of  the  pressure  difference  across  the  blade  (for 
constant  pitch  propellers,  this  is  equivalent  to  the  radial  component  of  bound  vorticity;  see  the  footnote 
following  Equation  (3-24))  and  an  infinite  integral  to  a  variable  point  on  the  blade. 

In  the  acceleration  potential  method  formulated  by  Sparenberg,  the  linearization  is  on  the  basis  of 
only  the  axial  component  of  Qq;  other  velocity  components  are  assumed  small  relative  to  it.  However, 
since  the  propeller  also  has  a  rotational  velocity,  it  is  not  clear  such  an  assumption  is  justified.  Hanaoka 
linearizes  about  the  total  free-stream  speed  and  converts  his  expression  for  the  linearized  acceleration  po¬ 
tential  to  one  for  the  velocity  potential.  For  steady  flow,  his  Equation  (3-8)  is  the  same  as  the  lifting 
portion  of  Equation  (3-21),  when  the  pitch  of  the  blade-reference  surface  and  the  shed  vortex  wake  are 
identical. 

The  formulation  on  the  basis  of  the  vorticity  distribution  is  now  almost  classical.  The  vorticity 
vector  is  specified  on  the  body  and  shed  vortex  sheet,  and  the  velocity  distribution  can  be  calculated  from 
the  Biot-Savart  Law.  For  a  surface  distribution  of  vorticity,  the  Biot-Savart  equation  gives  the  induced 
velocity  as 
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and  the  normal  to  this  sheet  is,  (from  Equation  (3-18)) 


CjJ  co 

N  - - -  n  - - -  (-  cos  i 3  /  +  sin  fie  a) 

cos  p  cos  p  u 

hence  the  vorticity  vector  is 

,  „  ar 

A  =nX(Q  -q  )  = - (sin  $1+ cos  p  cq) 

bo5 

where  tan  j3  *  ^/f2S J! 

Since  s  =  the  integration  over  the  shed  vortex  surface  is 


Id  7 


(5-24) 


where  /  is  given  by  Equation  (B-10)  in  Appendix  B;  this  is  the  form  given  previously  in  Section  B  of  this 
Chapter.  Thus  this  component  of  velocity  is  the  same  as  that  given  by  Pien,^  Morgan  and  Wrench,^  and 
others,  except  that  this  free-vorticity  starts  at  the  variable  trailing  edge  and  lies  on  the  surface  parallel  to 
the  undisturbed  free-stream  velocity. 

The  bound  velocity  on  the  blade  is  now  needed.  Since  most  of  the  previous  investigators  have 
derived  equations  appropriate  for  only  constant  pitch,  we  do  also.  This  means  that  the  zero-order  normal 
from  Equation  (3-10)  is 

.  r  + 

N=n  =  e2 


and 

Q+  -  f  =  IZi  +  oe~+  \xe2 

where  o  is  unknown,  and  7  is  given  by  Equation  (3-24). 

Hence  the  bound  vorticity  vector  is 


ye~+  oej 


(5-25) 


The  explicit  form  of  o  will  be  found  presently  but  now  the  velocity  induced  by  this  bound- 
vorticity  vector  is  desired: 


(oej  +  ye~)X(r-s0) 


d*l 
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comes 


where  Sq  is  given  in  Equation  (3-9),  which  for  constant  pitch  and  with  neither  rake  nor  warp  be- 
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y  sin  (6  -  <t>)(x-%]  sin  <pp)\ 
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(5-26) 
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where 
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Ir-s^l5  =  |  (x  -  tj  sin  i ppp  +  t2  +  Z52  -  2  r  co  cos  (  6  -  <p) p2 
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In  the  literature  it  is  common  to  work  with  6  = - rather  than  Since  for  Z o  =  constant 
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the  preceding  velocity  components  become 
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sin  (  Ob  +  6  -  <t>)  +  cos  <pp(x  -  <o  tan  <pp  0)  sin  (0  b+  0  -  <pj) 
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(5-27) 

Cont’d 


[r-s^  =  jfx  -  ZS  tan  ipp  dp  +  +  up  -  2  co  r  cos  ( 6b  +  6  -  0)  j  ^ 

In  the  literature  it  is  also  common  to  break  this  surface  integral  up  into  two  components 
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e~X  (r-s0) 
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(5-30) 


and  combine  4 y  with  the  integral  over  the  shed  vortex  sheet.  When  this  is  done  it  is  customary  to  take  the 
pitch  angle  of  the  shed  vortex  sheet  and  the  blade-reference  surface  as  the  |3 /  angle  computed  from  lifting¬ 
line  theory;  this  corresponds  to  the  hydrodynamic  pitch  introduced  in  Equation  (4-30). 

When  ^  and  qv  are  combined,  one  has  the  form: 
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Because  of  the  form  of  this  equation  it  is  also  customary  in  the  literature  to  speak  of  both  q  j  and  qv  as 
being  induced  by  a  trailing  vortex  system.  In  fact  some  authors  derive  their  equations,  postulating  that  the 
trailing  vortex  system  starts  on  the  blade. 

Because  of  this  confusion  regarding  the  o  component  of  bound  vorticity,  it  is  not  surprising  that 
the  literature  contains  several  expressions  for  it.  In  order  to  determine  a,  one  uses  the  property  that 
'a  distribution  of  vorticity  must  satisfy  the  equation 
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For  a  constant-pitch  helical  sheet,  with  the  vorticity  vector  given  by  the  component  7  in  the  radial 
direction  and  the  component  a  directed  in  the  surface  perpendicular  to  the  radial  direction  and  tangent  to  a 
cylinder,  one  can  determine  the  metric  coefficients,  Equation  (5-21),  and  reciprocal  unitary  vectors^  and 
thus  find  the  equation  for  divergence  on  the  blade 


where  P/2 it  =  Z5  tan  \p 
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and  on  the  shed  vortex  sheet 
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The  expression  in  the  wake  is  automatically  satisfied  since  T  =  F(oS)  only.  Now  one  must  find  the  form  of 
o  on  the  blade. 

From  the  first  expression,  one  obtains 

a(ZS,  0)  =  -  J  ~  fl/^Q  2  +  ^2  + 

eLE^  '  ' 

where  f(<3)  is  the  integration  “constant.”  The  value  of  this  constant  of  integration  is  obtained  by  requiring 
that  the  vorticity  vector  not  run  off  the  leading  edge  of  the  blade  into  the  irrotational  fluid.  That  is,  if  t_ 
is  the  tangent  vector  to  the  leading  edge 
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one  finds  that 
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where  y  (Z o,  6 )  is  the  radial  component  of  vorticity. 

Pien  and  Murray  do  calculate  the  velocities  with  an  expression  like  that  for  q  j '+  qv  in  Equation  (5- 
31).  Kerwin^  first  integrates  the  value  of  3a/0coon  the  blade  from  6  to  then  integrates  that  value 
over  the  blade  surface.  In  a  later  paper,  Kerwin  and  Leopold,^  change  the  formulas  to  the  form  of 


Equation  (5-31). 


D— EXAMINATION  OF  RADIAL  COMPONENT 
IN  BOUNDARY  CONDITION 

Although  the  formal  development  of  the  lifting-surface  and  lifting-line  equations  is  important 
academically,  it  does  not  in  itself  contribute  to  improved  propeller  designs.  However,  the  lifting-surface 
formulation  does  explicitly  consider  the  radial  velocity  component  in  the  boundary  condition  used  to 
derive  the  camber  distribution  in  Equation  (3-27).  Since  this  term  might  contribute  significantly  to  the 
engineering  aspects  of  propeller  design,  some  idea  of  its  effect  on  the  computed  camberline  must  be  found. 
In  order  to  examine  the  effect  of  this  term,  a  rather  crude  approximation  will  be  made  for  the  radial 
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velocity.  The  only  radial  velocity  components  found  in  the  literature  arc  those  calculated  by  Hough  and 

cn 

Ordway.  Their  results  are  for  circumferentially  averaged  velocity  components  due  to  a  lifting-line  with  a 
circulation  distribution' they  call  approximately  optimum;  for  the  radial  component  they  do  not  actually 
compute  values  at  the  lifting  line.  Nevertheless,  we  use  these  values  because  no  others  arc  available, 

CO 

Comparison  will  be  made  with  the  maximum  camber  ratios  calculated  by  Cheng  from  the  formu- 
lation  of  Pien.  Comparisons  of  the  averaged  axial  velocities  for  the  optimum  circulation  distribution  arc 
reasonably  close  to  those  computed  for  the  arbitrary  circulation  distribution  actually  used  in  the  calcu¬ 
lations  by  Cheng.  If  we  assume  the  radial  velocity  component  is  approximately  constant  and  approximately 
given  by  the  values  given  by  Hough  and  Ordway  at  the  point  closest  to  the  lifting  line,  then  the  camber  due 
to  the  new  term  is  parabolic  with  amplitude  for  symmetrical  blades 
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In  the  calculation  method  proposed  by  Pien,  u  i is  given  in  Equation  (4-51).  Cheng  tabulates  these 
values  as  well  as  the  chord-to-diameter  ratio.  The  following  table  lists  the  maximum  camber  computed  by 
Cheng  for  his  propeller  3916A  and  the  increment  computed  by  using  the  preceding  equation,  both  for  an  8- 
inch-diameter  propeller 
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Hough,  G.R.  and  D.E.  Ordway,  “The  Generalized  Actuator  Disk,”  Developments  in  Theoretical  and 
Applied  Mechanics,  Pergamon  Press,  Inc.,  New  York,  Vol.  II,  pp.  317-336  (1965). 
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Cheng,  H.M.,  “Hydrodynamic  Aspect  of  Propeller  Design  Based  on  Lifting-Surface  Theory'— Part  I— 
Uniform  Chordwise  Load  Distribution,”  NSRDC  Report  1802  (1964). 
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'  OJ' 

Cheng, 

AE  1 

R 

EcU*.  FiSs-  7  a"d  16 

c'max 

0.3 

0.088  in. 

0.0002  in. 

0.7 

0.034  in. 

'  —0.002  in. 

The  change  in  sign  in  AEc\max  reflects  a  change  in  sign  of  velocity  component  along  the  radius. 

For  this  propeller,  the  change  of  camber  is  negligible  for  the  Z o/R  =  0.3  radius.  This  is  because  of 
the  small  slope  of  the  pitch  curve  in  this  region.  For  the  0.7  radius,  the  change  in  camber  is  about  6  percent 
of  the  contribution  from  the  conventional  calculations  and  would  be  significant  for  propellers  of  larger 
diameter. 

Additional  information  about  the  importance  of  the  radial  component  can  be  surmised  by  con¬ 
sidering  the  effect  of  rake  or  warp.  The  effect  of  warp  is  especially  important  for  improved  cavitation  per¬ 
formance  and  reduced  alternating  thrust  when  operating  in  a  wake,  Miller^  and  Boswell. To  date,  in¬ 
vestigations  have  been  conducted  with  skewed  propellers  for  which  the  blade-reference  line  lies  in  the 
helical  surface  which  passes  through  a  straight-line  element  in  the  x  *  0  plane. ^  For  such  pro¬ 
pellers  the  normal  to  the  reference  surface  Nq  would  contribute  a  radial  term  approximately  like  that 
described  for  the  camber  effect;  however,  since  is  measured  from  x  =  0,  an  angle-of-attack  term  would 
also  occur.  Such  propellers  have  blades  extending  considerably  downstream.  This  overhang  requires  that 
the  propeller  hub  and  the  rudder  be  separated  more  than  normal.  If  instead  of  skewed  propellers,  only 
warped  propellers  were  utilized,  this  problem  would  not  occur.  Hence,  the  effect  of  warp  on  the  design 
will  be  considered  now. 

Skew  is  usually  given  as  a  warp  angle  with  the  understanding  that  the  blade-reference  line  lies  in 
the  helical  reference  surface  through  a  6 ^  line.  The  warp  angles  tabulated  by  Boswell^  are 
approximately: 


wfcS)  =  es 


s hub]  t 
~/f 


where  6 s  is  the  designated  amount  of  "skew.” 


^Miller,  M.L.,  "Experimental  Determination  of  Unsteady  Propeller  Forces,”  Seventh  Symposium  on 
Naval  Hydrodynamics,  Rome  (1968). 

^Boswell,  R.,  "Design,  Cavitation  Performance,  and  Open-Water  Performance  of  a  Series  of  Research 
Skewed  Propellers,”  NSRDC  Report  3339  (1971)- 
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That  portion  of  the  boundary  condition  of  Equation  (3-27)  to  be  examined  is 
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If  the  preceding  expression  for  warp  is  used,  and  if  the  radial  component  is  again  approximated  by 
w*(Zj),  then  warp  contributes  an  angle  of  attack  (or  pitch  correction)  to  the  camberlinc  shape.  The 
tangent  to  this  pitch  angle  is  given  by  the  preceding  expression 
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For  lack  of  other  data,  we  approximate  the  \p  values  at  the  warped  line  by  the  |3?-  values  for  a  straight 
blade-reference  line. 

By  applying  these  approximations  to  the  five-bladed  propeller  considered  previously,  the  pitch-angle 
corrections  would  be 


gT 

A  a  in  Degrees 

R 

0.3 

-0.2 

0.7 

0.5 

These  corrections  are  considered  significant. 

Because  of  the  change  in  sign,  the  preceding  corrections  might  tend  to  cancel  in  their  effect  on 
thrust.  However,  they  could  be  important  for  cavitation  performance  which  is  one  of  the  main  reasons  for 
designing  a  propeller  rather  than  selecting  one  from  a  series. 
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E— DESIGN  PROCEDURES  (HYDRODYNAMIC) 

In  current  design  techniques  by  Morgan,  Silovic,  and  Denny*^  and  Kerwin  and  Leopold,^  a  combi* 
nation  of  lifting-surface  and  lifting-line  results  is  used.  Briefly  the  mechanics  of  this  procedure  are  to  use 
the  second-order,  lifting-line  equations  for  the  outer  flow  to  find  the  induced  velocities  at  the  blade- 
reference  line,  circulation  distribution  and  thrust.  The  inner  flow  is  ignored,  and  hence  an  integral  equation 
must  be  solved.  This  is  usually  done  by  assuming  either  a  tan  P-  or  circulation  distribution  and  scaling  the 
chosen  quantity  to  eventually  give  the  required  thrust  (or  power).  The  resulting  circulation  distribution  is 
used  to  scale  the  chosen  radial  component  of  vorticity  7;  see  Equation  (3-31).  The  lifting-line  induced 
velocities  are  assumed  to  be  approximately  the  average  chordwise  results  from  a  lifting-surface  calculation. 
They  are  used  to  approximate  the  velocities  induced  by  the  shed  vortex  wake  and  are  used  to  give  the 
reference  surface  ^ Z5).  Hence  everything  is  known  to  compute  the  lifting-surface  results,  partially  taking 
into  account  second-order  effects.  This  procedure  has  worked  remarkably  well  for  the  performance  of  con- 
ventional  propellers,  Cox.  '  On  the  other  hand,  it  seems  to  be  difficult  to  design  a  bad  propeller  since 
even  the  earlier,  semi-empirical  design  procedures  such  as  HilPsDZ  produced  propellers  which  developed  the 
desired  thrust.  It  is  in  other  areas,  especially  cavitation  inception,  that  the  more  refined  analysis  embodied 
in  the  lengthy  lifting-surface  calculations  are  expected  to  have  a  distinct  advantage.  Unfortunately,  few  ex¬ 
perimental  comparisons  of  the  type  necessary  to  judge  the  adequacy  of  the  various  hydrodynamic  theories 
are  found  in  the  literature.  However,  when  cavitation  occurs  first  at  the  leading  edge  for  the  design 
condition,  the  design  procedure  has  not  been  sufficient  to  produce  the  specified  chordwise  pressure  distri¬ 
bution.  Such  is  the  case  with  a  recent  series  of  skewed  propellers.' ^  Some  differences  in  the  numerical 
results  of  various  lifting-surface  formulations  are  given  by  Johnson.  3 

Although  the  combination  of  lifting-line  and  lifting-surface  results  is  quite  practical,  it  is  not 
entirely  consistent.  Specifically,  the  expansion  for  the  circulation  distribution  in  the  lifting-line  analysis 
can  hardly  be  expected  to  apply  for  e  -*  7.  However,  for  elliptic  wings,  Van  Dyke  ’  has  shown  that  the 
lift  computed  from  the  integral  equation  is  quite  close  to  the  exact  value  over  the  entire  range  of  aspect 
ratios.  Thus  the  lifting-line  results  might  be  adequate  for  parametric  studies  since  calculations  can  be 
quickly  done.  For  final  design,  though,  the  lifting-surface  equations  should  be  used  since  they  will  give  a 
geometry  better  able  to  meet  design  conditions. 

Another  utilization  of  the  lifting-line  equations  is  to  approximate  $p(Z5)  and  the  position  vector  of 
the  shed  vortex  sheet.  The  determination  of  the  shed  vortex  sheet  for  the  lifting-surface  analysis  by  using 


f)  1 

Morgan,  W.B.  et  al.,  “Propeller  Lifting-Surface  Corrections,”  Society  of  Naval  Architects  and  Marine 
Engineers  Transactions,  Vol.  76,  pp.  309-347  (1965). 
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(he  lifting-line  induced  velocities  is  probably  not  accurate  near  the  lifting  surface  where  it  is  most  important; 
in  any  event,  the  current  utilization  is  only  an  approximate  second-order  correction  for  a  portion  of  the 
second-order  terms.  Other  second-order  terms  should  be  calculated  to  determine  their  magnitude  also. 

To  use  the  consistent  first-order  lifting  surface  calculations,  the  following  procedure  is  recommended. 
In  the  first-order  theory  the  approximate  position  of  the  shed  vortex  sheet  is  independent  of  the  induced 
velocities  and  is  known.  The  blade-section  pitch  angle  \p  is  first  approximated  by  the  |3  angle.  The  mag¬ 
nitude  of  the  circulation  distribution  needed  to  produce  the  thrust  can  be  calculated  from  Equation  (3-23). 
The  propeller  geometry  can  now  be  calculated.  The  angle  of  attack  of  the  section  can  be  used  to  define  a 
new  ,  and  the  calculations  can  be  repeated.  The  second-order  results  could  also  be  calculated  if 
necessary.  In  particular  the  second-order  thrust  in  Equation  (3-29)  could  be  calculated  from  the  first-order 
data  so  that  an  adjustment  in  the  circulation  could  be  made  to  meet  the  given  thrust  or  some  other  criteria 
used  to  iterate  for  the  final  solution. 

Obviously,  these  alternative  procedures  would  be  more  involved  than  the  presently-used  techniques 
but  unless  an  entirely  consistent  approach  to  the  design  problem  is  taken  one  does  not  know  the  degree  of 
approximation  involved  in  the  present  calculations  or  what  potential  exists  for  design  based  on  the  lifting- 
surface  formulation. 

In  the  equations  for  design  the  chordwise  component  of  velocity  difference  (or  the  spanwise  com¬ 
ponent  of  vorticity)  7  is  assumed  to  be  given  as  well  as  the  total  bound  circulation  curve  Y(Z5).  This  is 
compatible  with  information  which  a  designer  should  supply  for  a  propeller  designed  from  cavitation  con¬ 
siderations:  Two  general  criteria  for  the  avoidance  of  cavitation  are  that  no  suction  peaks  occur  on  the 
blade  and  that  the  shed  vorticity  be  distributed  as  uniformly  as  possible  in  the  spanwise  direction.  The 
first  conditions  met  by  selecting  the  functiona  form  of  7;  the  second,  by  the  functional  form  of  Y(Z5), 
since  the  gradient  of  Y  gives  the  shed  vorticity.  The  number  of  blades  and  the  blade-area  ratio  are  also 
part  of  these  considerations  since  the  average  pressure  times  the  blade  area  determines  the  load  which  must 
be  sufficient  to  produce  the  required  thrust. 

F— RECOMMENDATIONS 

Obviously  the  first  recommendation  is  to  include  the  radial  velocity  component  in  the  boundary 
condition  for  propellers  designed  with  variable  geometric  pitch,  warp  or  rake.  Although  the  inclusion  of 
this  term  is  not  expected  to  have  a  significant  effect  on  the  thrust,  it  could  lead  to  propellers  with  better 
cavitation  performance.  Since  cavitation  considerations  are  one  of  the  main  reasons  for  designing  a  pro¬ 
peller  rather  than  selecting  one  from  a  series,  and  since  a  crude  analysis  shows  the  term  to  have  a  significant 
effect,  detailed  numerical  calculations  should  be  undertaken. 

Since  heavily-loaded  propellers  are  becoming  more  common  in  practice,  the  second-order  lifting- 
surface  formulas  should  also  be  programmed  for  design  calculations.  One  of  the  interesting  aspects  of  such 
investigations  would  be  determination  of  the  position  vector  of  the  lines  of  constant  circulation  in  the  shed 
vortex  sheet.  Another  interesting  outcome  of  these  calculations  would  be  an  evaluation  of  the  magnitude  of 
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the  second-order  terms.  Presently,  it  is  assumed  that  the  shed  vortex  sheet  needs  to  be  described  better 
than  by  Qualitatively  this  is  saying  that  the  infinite  integral  of  ~  0  (e)  over  the  shed  vortex  sheet 
produces  a  zero-order  contribution.  The  proposed  calculations  would  evaluate  this  assumption. 

The  formal  analysis  described  here  should  also  be  extended  to  include  vorticity,  i.e.,  an  axisymmetric 
free-stream  velocity,  unsteady  flow,  and  supercavitating  propellers.  Since  unsteady  flow  with  vorticity  is  the 
standard  operating  condition  of  propellers  behind  a  ship,  these  two  are  important  areas  of  research. 


SUMMARY  AND  CONCLUSIONS 

The  mathematical  development  presented  in  the  preceding  chapters  has  been  based  on  the 
equations  of  motion  and  certain  mathematical  relations  which  have  permitted  an  exact  formulation  to  be 
constructed.  Application  of  formal  perturbation  procedures  has  yielded  linearized  solutions  appropriate  for 
numerical  calculations.  In  Chapter  5,  the  velocities  computed  from  the  lifting-surface  potential  are  shown 
to  equal  those  computed  by  the  Biot-Savart  Law.  These  expressions  are  then  shown  to  equal  those  found 
in  the  literature  for  the  approximate  geometrical  conditions. 

An  examination  of  the  literature  did  not  reveal  previous  explicit  consideration  of  the  radial 
velocity  component  in  the  first-order  lifting-surface  perturbation  problem,  Equation  (3-27).  Of  course, 
other  investigators  realized  that  such  a  term  occurred  for  variable  pitch  propellers  but  did  not  examine  the 
consequences  of  neglecting  it  for  arbitrary  pitch  distribution.  For  the  example  considered,  the  contribution 
of  this  term  was  as  much  as  6  percent  of  the  camber  ratio  determined  by  the  conventional  analysis. 

In  the  lifting-line  analysis,  only  propellers  without  rake  and  warp  were  considered.  For  these 
propellers,  the  two-dimensional  sections  were  defined  normal  to  the  straight  blade-reference  line.  However, 
to  the  order  of  the  approximation  considered,  they  could  equally  well  have  been  considered  as  defined  on 
cylinders  as  was  done  in  the  lifting-surface  problem. 

A  method  of  design  based  only  on  the  lifting-surface  formulation  has  been  discussed.  This  pro¬ 
cedure  makes  no  use  of  the  lifting-line  results  and  would  involve  greater  computing  effort  since  an 
iteration  is  appropriate. 

Extension  of  the  method  utilized  in  this  work  to  the  problem  of  unsteady  flow  and  shear  flow  is 
recommended.  In  addition,  numerical  results  should  be  obtained.  The  numerical  results  could  include  wake 
trajectories  for  the  shed  vortex  sheet. 
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APPENDIX  A 

POTENTIAL  FOR  LIFTING  LINE  DERIVED  FROM  THE 
BIOT-SAVART  LAW 


Traditional  propeller  lifting-line  theory  (Moriya,^  Lerbs,^  Morgan  and  Wrench,^)  evaluates  the 
velocity  induced  by  the  vortex  system  from  the  law  of  Biot-Savart.  The  vortex  system  considered  is  one  for 
which  an  individual  vortex  is  an  element  of  circulation  of  strength  ^T/c/coshed  at  radius  to;  then  by  the 
conservation  of  vorticity  this  element  remains  bound  in  the  line  and  leaves  it  at  the  axis  of  rotation.  It  is 
not  clear  from  the  literature  how  to  interpret  the  vortex  system  when  the  lifting  line  extends  from  a  finite 
radius  off  the  axis  to  the  propeller  tip.  (However,  see  Appendix  B.) 

Lifting-line  theory  for  planar  wings,  Ashley  and  Landall,^  uses  a  model  which  consists  of  an 
elemental  horseshoe  vortex  of  strength  I\  This  is  composed  of  two  free  vortices  of  equal  but  opposite 
strengh  a  distance  A z  apart,  joined  together  at  the  blade  by  a  bound  vortex  element  of  the  same  strength. 
Such  a  system  is  shown  in  Figure  8  for  a  propeller  lifting  line. 


Figure  8  —  Vortex  Element 
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A  sum  of  these  elements  distributed  along  the  radius  of  the  lifting  line  then  approximates  the 
vortex  system.  In  Ashley  and  Landahl,^  the  potential  for  elements  an  infinitesimal  distance  apart  is 
derived  by  integrating  the  potential  of  a  doublet  over  the  length  of  the  element.  Here  it  will  be  derived 
from  the  Biot-Savart  Law  applied  to  the  element.  The  path  of  integration  is  shown  in  Figure  9. 


Figure  9  —  Path  of  Integration 


Let  this  path  be  L.  The  velocity  induced  by  the  element  is 


One  of  Stokes*  theorems  for  vector  integration  is 


X  Vyx  u  ds 


(Al) 


(A2) 


*See  also  Lamb^  and  Kochin  et  al.^ 

^Lamb,  H.,  “Hydrodynamics,”  Sixth  Edition,  Dover  Publications,  Inc.,  New  York  (1964). 


100 


where  s  is  an  area  bounded  by  the  curve  c.  Thus 


(nX  VS)X 


r-s 


ds 


where  s  is  the  area  bounded  by  L.  But 


(nX  VJX 


r-s 


S  r-s\3 


Hence,  since  V  •  (r/r*)  =  0 ,  there  remains 


r-s 


\  f 


(A3) 


(A4) 


r-s 


(n  X  Vs)  X 


since  n  =  n(s)  only. 

Thus  the  induced  velocity  is 


(A5) 


Hence,  by  definition,  the  potential  is 


l  rr  £L 

JJ  s  |r_£h? 


•  n  ds 


°  n  ds 


(A6) 


(A7) 


See  Lamb^^  Equation  4,  page  212. 

The  shed  vortex  sheet  is  envisioned  as  being  composed  of  a  radial  distribution  of  these  elements  so  that  the 
distance  between  elements  goes  to  zero  and  the  width  of  an  element  goes  to  zero.  In  this  limit  the  vector 
element  of  area  becomes  the  negative  of  that  in  Equation  (248): 


n  ds  =  - 


>  3?\ 

—  X - \dzQda 

9cr  9  zqJ 
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and  the  potential  for  the  total  system  becomes 


l 


Z~l 

2 

b=0 


When  the  zero-order  position  vector  for  the  trailing  vortex  system  is  substituted  into  Equation  (A8), 
Equation  (4-5)  is  obtained. 

As  an  application  of  Equation  (A7),  the  potential  of  a  circulation  element  for  planar  wings  can  be 
derived.  For  wings,  the  approximate  location  of  the  lines  of  constant  circulation  in  the  shed  vortex  sheet 
is  on  the  x,  y-plane  at  constant  y  values.  Hence 


and 


n  ds  =  -  k  dx  Ay  +  0(  Ay^) 


r  Ay  r 

«=—  I 


z  dx, 


4ir 


FA y 


0 


(x-x0)2  +  (y-y0F  +  z2]  3>2 


+  Of  Ay2) 


411  (y-ygf  +  z2 


1  +  - 

Vx2  +  (y-y0)2  +  ?2_ 


which  is  Equation  (5-35)  in  Ashley  and  Landahl.^2  This  can  be  rearranged  to  give 


TA y 


(x2  +  fy-y2)  +  z2)1!2  +  x 


0  = 


4ir  (x2  +  (y-y0)2  +  z2)  -  x2  (x2  +  (y-y0)2  +  z2)1!2 

+  Of  Ay2) 


(A9) 


TA^ 


47T  fx2  +  (y-yg)2  +  z2] 1!2  -  x  [x2  +  fy-yg)2  +  z2] 7/2 


+  Of  Ay2) 


(A  10) 
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which  is  the  form  derived  by  von  Karman  and  Burgers  (p.  1 28),  by  Truckenbrodt,  and  by  Friedrichs 

(Equation  15.4).* 

Returning  to  the  propeller  problem,  one  obtains  Equation  (4-5)  by  substituting  the  zero-order 
position  vector  of  curves  of  constant  circulation  from  Equation  (4-3)  into  Equation  (A8).  Since  only  the 
kernel  of  Equation  (4-5)  is  of  interest,  it  can  be  considered  alone: 


Z-l  oo 


G(r,z0)=  2  I 

b=0  o  HI*  Vda  9z" 


r-i  (  ar  3? 

x—  I  da 


(All) 


and  the  zero-order  approximation  to  G  is 

Z-l  -oo 


G(r,z0)  -  2  j 


~z0 


b=0  0 


a  u  a  yy 

3;c  £Izq  doj  r 


(4-6) 


r  =  - — 7/^  +  +  Zq  -  ~2zq  oj  cos  z'  0  -  7  -  6^j 


A  similar  integral  has  been  evaluated  in  the  determination  of  induced  velocities  at  a  lifting  line.  For  the 

induced  velocities  at  the  lifting  line,  Moriya**  was  able  to  integrate  the  expression  analytically  to  obtain  an 

infinite  series  of  products  of  modified  Bessel  functions.  Although  the  first  impression  is  that  replacing  an 

1 3 

infinite  integral  by  an  infinite  series  is  not  much  of  an  improvement,  Lerbs  showed  that  the  series  could 
be  approximately  summed  and  thus  calculations  could  be  quickly  made.  An  excellent  explanation  of  the 
integration  process  for  the  induced  velocities  at  a  lifting  line  and  the  numerical  approximation  to  the  series 
is  given  by  Morgan  and  Wrench.  Included  in  that  article  is  an  improved  approximation  for  the  series 
summation  which  was  developed  by  Wrench. 

The  rather  circuitous  method  discovered  by  Moriya**  will  be  applied  to  the  integral  in  Equation  (4- 
6).  Unfortunately,  the  kernel  cannot  be  evaluated  at  an  arbitrary  point;  however,  the  form  far  downstream 
is  found,  and  it  is  demonstrated  that  velocities  from  this  expression  are  twice  those  at  the  lifting  line  as 
derived  from  the  Biot-Savart  Law,  as  would  be  expected. 

ft! 

The  expression  for  Gfr.Zg)  is  given  by  Equation  (4-6).  From  the  Lipschitz  integral,  Watson,  p.384, 
one  obtains 

^von  Karman,  T.  and  J.M.  Burgers,  “General  Aerodynamic  Theory -Perfect  Fluids,”  Vol.  II  of  Aerodynamic 
Theory,  edited  by  W.F.  Durand,  Springer;  also  reprinted  by  Dover  Publications,  Inc.  (1963). 

^Truckenbrodt,  E.,  “Das  Geschwindigkeitspotential  der  tragenden  Flache  bei  inkompressibler  Stromung,” 
Zeitschrift  fur  angewandte  Mathematik  und  Machanik,  Berlin,  Vol.  33,  pp.  165-173  (1953). 

^Watson,  G.N.,  “Theory  of  Bessel  Functions,”  Cambridge  University,  England  (1944). 
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(A  12) 


]_ 

r 


2zq  cocos  ( 0  -  7  -  Of,)  dt 


From  one  of  the  addition  theorems,  Watson, 


67 


P- 


128, 


-  2z0  cocos  ( 6  -  7  -  Qb) 


oo 

=  2  km  JJtoS)  Jm(tzQ)  cos  (  0  -  a) 

m=0 


(A  13) 


where  kfn  is  Neumann’s  factor,  (km  =  0 ,  m  =  0;  kfn  =  2,  m  >  0).  If  one  assumes  uniform  convergence  for 
the  integral,  the  order  of  integration  can  be  switched  and  one  obtains 


where 


CO  CO 


c/r,v=  s  j 


(A  14) 


m=0  0 


Z-l 

2 

b=0 


r°° 

3*  30  /  I  e  \  /  cos  m( 0  -  6 fr  ~  y)  dy  (A15) 

0 


Two  cases  occur,  one  for  x  <C  0  and  one  for  jc  >  0.  For  x  <  0  the  positive  sign  is  used.  For  x  >  0,  two 

£lv 

regions  exist,  one  for  y  and  one  when  it  is  greater.  Only  the  case  for  x  >  0  will  be  considered  since 

it  permits  the  form  of  the  potential  far  downstream  to  be  obtained.  Thus  for  x  >  0 


Z-l  ( 

SEm  lr-'z0’V  =  2  \-z0^:  + 

b=0  \ 


U  3 


J. 


3x  3  6 

Uyt 


Uyt 

-"J/TT 


cos  m  (  0  -  7  -  Ojj)  dy 


+  ext  i  e 


n 


cos  m  (d  -  7  -  0 b)  dy 


£lx 

U 


(A  16) 
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The  integrals  appearing  in  Equation  (A16)  can  be  found  by  integrating  by  parts  until  the  same  expression 
occurs  again  and  then  by  rearranging  the  equation.  The  resulting  expression  for<i^  is 


The  portion  of  the  integral  involving  the  exponential  disappears  when  the  form  far  downstream  is  examined. 
The  remaining  portion  of  the  integral  can  be  integrated  using  a  result  on  page  429  of  Watson.  /  The 
resulting  expression  far  downstream  is 

«i,j  (A>« 

where  Am  is  a  combination  of  modified  Bessel  Functions 


2£l 

IT 


z-i 

2 

b=0 


2  mk 
m=0 


m  i 


U 


u  \  ( 

- |  sin  m  I  6 

^ zo)  \ 


Ux 

77 
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ftoom\  (SlzOm 


'm  \  u  )  Km  \  u 


z()>o5 


(nz(/”\  .  /ft cow 

L\  u  I  Km  \  u 


Zq<CO 


The  expression  for  0  far  downstream  is 


A 

i-J  V(z0jGoo(r,z0jdz0 


which  can  be  integrated  by  parts  to  obtain 


%)  di)dz 


After  some  manipulation,  one  finds 


R  Z-l  °o  . 

c/r  /  ft* 

X  2  Kms,nm  \e-T 

1  0  b=0  m=0 


h)A'mdzO 


(fto om\ 

~) 


Slz^ri 


faZ>m 

\  u  )  k,n  \  U  ) 


Zq  >  CO 


Zq<c5 


The  expression  for  the  axial  velocity, — ;  at  6  +  0,  x  -*  0  can  be  easily  obtained.  After  summing  over  b,  the 

dx 

expression  is 


za/u  f  dr  ~  mZ^zo  , 

W—2 —J  Td*0  2  Amz 

0  0  m-0 


Tlais  is  twice  the  value  of  the  axial  velocity  at  the  lifting  line  as  determined  from  the  Biot-Savart  Law,  as 
expected. 


106 


An  expression  for  the  portion  of  the  potential  involving  the  exponential  integral  can  be  obtained  by 
expressing  the  Bessel  Functions  as  infinite  series.  The  integral  then  involves  a  triple  infinite  summation  of 
terms  with  an  integral  now  of  the  form 


I 


oo 


0 


e 

7w 


This  expression  can  be  integrated;  see  Gradshteyn  and  Ryzhik  p.  313,  3.356(1)  and  (2)^  [n  terms  of  sine 
and  cosine  integrals  plus  another  finite  summation.  The  resulting  expression  is  quite  cumbersome  and  is 
therefore  not  given.  No  use  was  made  of  it  in  the  analysis. 


^Gradshteyn,  I.S.  and  I.M.  Ryzhik,  “Tables  of  Integrals  Series  and  Products,”  (Translation  of  the  Fourth 
Edition,  Edited  by  A.  Jeffrey),  Academic  Press,  Inc.,  New  York  (1965). 
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APPENDIX  B 

INTERMEDIATE  EXPANSION  OF  OUTER  POTENTIAL 


Equation  (4-32)  states  the  second-order  potential  for  the  outer  flow.  This  expression  is  valid  for  a 
fixed  point  in  the  flow  field  as  the  parameter  e",  equals  chord-to-diameter  ratio,  goes  to  zero.  Thus  the  outer 
flow  consists  of  a  singularity  system  distributed  along  the  radial  lines  8^  0  <  oo </?. 

There  are  only  two  essential  systems  of  singularities:  one  from  the  circulation  distribution  and  one 
from  the  dipole  distribution.  The  second-order  contribution  to  the  circulation  will  be  included  in  the 
single  integral  involving  the  circulation.  The  two  integrals  will  be  (j)j  and  02>  where  (j)j  is  given  by 
Equation  (4-5),  and  (j>2  is  given  by  Equation  (4-33).  The  (j)j  term  will  be  considered  first. 

Since  considerable  difficulty  was  encountered  in  an  attempt  to  expand  the  potential  function  for 
the  lifting-line,  the  velocity  expression  was  expanded  and  from  it  the  expansion  for  the  potential  was  con¬ 
structed. 

In  Appendix  A,  the  velocity  from  the  Biot-Savart  Law  was  shown  to  be  derivable  from  a  potential. 
Now  it  is  convenient  to  return  to  the  Biot-Savart  Law  in  order  to  find  the  intermediate  expansion.  How¬ 
ever,  some  rearrangement  is  in  order  before  the  appropriate  form  for  expansion  is  obtained.  The  vortex 
system  is  envisioned  as  composed  of  a  distribution  of  elements  shown  in  Figure  8  along  each  of  the  lifting 
lines.  The  width  between  the  free  trailing  portions  of  a  vortex  element  is  constant  for  all  elements,  and  the 
elements  are  infinitesimally  close  together.  The  path  tiZconsists  of  the  segment  along  the  lines  and  the 
two  free-trailing  paths.  Two  adjacent  trailing  elements  then  have  a  common  path  of  integration.  Hence  the 
total  induced  velocity  for  N  of  these  elements  is 


where  TQ  =  T N+j  =  0 

j^n  is  the  common  path  between  the  adjacent  n  and  n+1  trailing  vortices. 

In  the  limit  as  Aco^  0y  and  the  distribution  of  circulation  becomes  continuous  with  T(0)  =  T(R)  =  0, 
the  expression  for  the  velocity  becomes 
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/  Z'1  rR  r-  ~  zOc-Z$(°lJ 

q=—  2  f  r(z0) - -XetfObldZf, 

•J  I  r  - 


-  4  71 


b=0  0 


I-  zO-(jd®b^ 


(B2) 


,  Z~'  (R<iri:0l 


A  /ft 

—  2  f 

47T  J 


^  Z 


d  Z 


*>=0  0 
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I 


r-  f 


°  1  |r-fl* 


Hence  the  first-order  velocity  distribution  results  when  T  ^  (zq)  +  *  an^  f  ^0  +  are  su^ 

stituted  into  this  equation.  Since 


U 


$0  =  z0e-o5(°)+nf°-  °b)i- 


(B3) 


one  finds 


d  c; 


■  OJ  V 

d  $0  =  zr -  d0  +—  d0  L 


n 


u 


(zoe-e +  n  i7d(? 


and  since 


U 


\r-y=l/{x--(a-  eb 


)\2  +  o52  +  Zq  -  2  CO  zq  cos  (  0  -  cty 


one  finds 


(r  -  l0)*  d  f0- 


j_  .f  ^  +  z0  sw  ay  -  Z0  cos  a  fx  -  —  ya  -  0 b;y  |  £ 
+  |  (Z  -  Zq  COS  ay  +  Zq  SW  CL  (x  -  ~  (C C  -  0^  j  / 


+  |zfl  cocos  /"  0  -  ay  -  Zfl  1 1 


d  a 


(B3) 


Hence 
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Z-l  R 


-  “  r  -  (jj  sin  (  a 

q=—  I  r  - 

4lT  [x2  +  (v  + 


ZSsin(6  -  0  fr) x  cos  O^j-xsin  O^k 


b=0  0 


[x2  +  (y  +  zQ  sin  Oh?  +  (z  -  Z0  cos  0bPl3l2 


e  Z-J  rR  or 

+5T  2  J  dl°lT0 

b=0  0  " 


Ho  " 

<L  -  $VXXT  d  a 

rOO  ~u  3a 

6h  [(x  - -p-  (a  -  6j2  +  oS2  +  zr 2  -  2  cSzq  cos  (0  -  a)/ 3!2 


Consider  first  the  expression  from  the  bound  vorticity.  To  find  the  intermediate  expansion  one  substitutes 
the  intermediate  variables  into  the  expression  and  rearranges  it  into  an  asymptotic  sequence.  By  inspection 
of  the  denominator,  one  can  see  that  a  Taylor  expansion  is  appropriate,  except  for  the  b  =  0  lifting  line. 


r  xgl+gy i 

r(Z0>-~2^2  ~=2 - 

(e 2  x*  +  vz)  +  (z- 


$Btex,gy,z)=  — 


[g2  (x2  +  y2)  +  (z  -  Zrsj2]3!2 


ci  Z-l  R 


el  £  z  sin  dfod  Zq 

—  X  I  r (zq)  ~  -  . 

4lT  b=l  J0  !(z0 sin  Qb>  +(z  -  z0  cos  eb'i  2 


+  o  (eg) 


By  symmetry  the  second  term  cancels,  i.e.,  sin  6j=-  sin  Q z_j,  etc.  The  first  integral  can  be  expanded 

2  2—7  —7 

by  a  binomial  series  in  each  region  bordered  by  (z  -  Zq)  =  g  ( x *  +  y*).  The  result  is  given  by  Ogilvie, 
Equation  (2-43)^ 


--  €  _  xj+yl 

qB  (g  g  y>  z)  -  ~z  r(z)  -  -  ^ . ~7  +  °( 

2ng  x2  +y2 


Temporarily  we  convert  to  outer  variables  . 


?r (z)  y  i  -  xj 

q-*(x’y’z,—-jT7  +  0(e) 


_  7  y  y  1  -  x  ] 

and  note  that  7  tan  — = — : - —  ,  hence 


x2  +y7 
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eV(z)  _,y 

qjj  (x,  y,  z)  =  -  -  V  tan  ‘ —  +  0(e) 


From  which  one  can  conclude  that,  to  within  a  function  of  z,  the  expansion  of  0  contains  a  term: 


=  =  er  ,y  7 

<t*R  (g  x,  g  v,  z)  =  -  —  tan  l— +  0(  eg *) 
°  2n  x 


due  to  the  bound  vortex. 


Now  the  second  part  of  the  expression,  the  velocity  due  to  the  trailing  vortices,  must  be 
determined.  Letting  y=  a-  0  b,  one  finds 

_  z-l  rR  dr  r°° 

2  J  ldl 

b=0  J0  0  0 

(  (z0  cocos  (0  -  y-  0b)  -  z02)  i_ 


+  /  —  (z  -  z o  cos  (7 +  V  +  z0  sin  (7  +  (x~-^  7)1  L 

U  U 

-  /—  (y  +  z0  sin  (y  +  d b)  +  zQ  cos  (y  +  0b)  (x  -  —  y)j  k 

/(*  +  +  Zq  sin  (y+  6b)j  2  +  (z  -  zQ  cos  (y  +  0h))2/Sl'2 


Since  the  denominator  is  of  positive  definite  form,  it  has  a  zero  only  for  x  =  y=  0 b  =  y  =  0,  z  =  Zq.  As 
shown  by  Moriya^  *  and  Morgan  and  Wrench,^  the  velocity  vector,  as  x  =  g  x  and  y  =  g  y  go  to  zero,  is 


found  to  be 
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na  (-,  z0.  1)  =  Z  zo  cos  (7  +  Oh)  z()~ 

U 

H,  (z.  z(),  7)  =  ~/z  -  z()  cos  (7  +  Oh)  -  7  z()  sin  (7  +  Of,)/ 

Uzo 

Hz  (z,  zQ,  7)  =  -  — -  (sin  (7+  Oj})  -  7  cos  (7+  Oh)/  z() 

u  L 

and  where  the  bar  through  the  integral  means  the  singular  point  is  excluded.  Except  for  notation,  these  arc 
identical  with  the  expressions  given  in  Equations  (11)  through  (13)  of  Morgan  and  Wrench,^  Equation  (7) 
of  Moriya*  *  and  the  equations  in  Appendix  2  of  Lerbs.^  As  shown  by  these  authors,  when  U/£l  is  con¬ 
stant,  the  infinite  integral  can  be  performed  analytically  for  the  axial  and  tangential  component  of  velocity. 

From  the  previously  described  form  of  the  velocity  vector  in  the  intermediate  region,  one  con¬ 
cludes  that  the  expansion  for  the  potential  due  to  the  trailing  vortices  has  the  terms 

-  -  -  (1)  =  (1)  fZ  (1) 

<l>T(gx,gy,z)  =  eg(xua  +yut  )  +  e\  uz  dz  (B12) 

Hence  the  expansion  for  the  potential  due  to  the  vortex  system  has  the  following  form 


?  <pj  (g  x,  g  y9  z)  =  <t>B  + 

e  T (z)  i  y 

=  F  (z;  e)  - - tan*1  — 

2tt  = 

+  eg(x  Ja!)  +  y  u{  )  +  H.O. T.  (B 1 3) 

The  function  F(z;e)  contains  other  unknown  integration  constants  which  are  functions  of  z  in  addition  to 
the  integral  of  The  construction  of  the  potential  expansion  from  the  velocity  components  is  not 
required  since  the  matching  could  be  done  with  the  velocity  components  but  it  is  simpler  to  work  with  the 
one  potential  expansion  rather  than  the  /  and  /  velocity  components. 

The  expression  for  <j>2  valid  in  the  intermediate  region  is  now  to  be  found.  For  the  expansion  of 
this  expression,  it  is  not  necessary  to  examine  velocity  components,  since  a  direct  expansion  of  the 
potential  is  possible.  Substitution  of  intermediate  variables  in  the  expression  for  0^  §^ves 


02  (s  x,  s  y< z) =  —  £ 


Z~1  rR  gx  Hj  (z0)  +  (gy  +  z0  sin  Bb)  n2  (zQ)  c 


b=0  0 


l(g  x)2  +  (gy  +  zn  sin  dt,)2  +  (z  -  zn  cos  d,,)2/^2 
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For  b  ^0,  the  expansion  is  regular,  and  a  Taylor  expansion  is  possible.  The  term  independent  of  g  in  the 
expansion  is 


zQ  sin  0b  \x2  (z0)  d  20 
[(zQ  sin  6  b)2  +  (z  -  z0  cos  6  b? f3^2 

Because  of  the  symmetry  of  the  lifting  lines,  the  summation  of  this  term  is  zero.  The  term  corresponding 
to  b=0  has  the  following  expansion;  see  Ogilvie,^  Equation  (2.43) 
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(B 1 5) 


Hence 
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(B16) 


Combining  this  with  the  expansion  of  the  vortex  terms,  one  finds  that  the  expansion  of  the  potential  in 
the  intermediate  region  has  the  form: 


<t>(g  x,  g  y,  z)  =  Fq  (z;e)  -  e 


V(z) 

2k 


/  y 

tan  — + 


_  (1)  =  (I) 

+  eg  (x  ua  +y  u(  )  + 


~2  x  Hj(z)  +y  V-^z) 

2k g  ^2  +  y2 


+  H.O.T. 


(4-34) 


This  is  the  expression  stated  in  Equation  (4-34). 
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